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Final  Technical  Report  for  AFOSR-86-0177  Grant 

UNSTEADY  FLOW  STRUCTURE  FROM  SWEPT  EDGES 
SUBJECTED  TO  CONTROLLED  MOTION 


Program  Monitors:  Captain  H.  Helin  and  Dr.  J.  McMichael,  May  1,  1986  -  October  31, 
1990 

Principal  Investigator:  Professor  Donald  Rockwell,  Department  of  Mechanical 
Engineering  and  Mechanics,  Lehigh  University,  Bethlehem,  PA 


SUMMARY 


This  program  addresses  the  unsteady  flow  structure  and  loading  of  delta  wings 
subjected  to  controlled  pitching  motion.  Efforts  have  focussed  on  three  primary  areas: 
generation  of  computer-aided  techniques  for  quantitative  interpretation  of  flow 
structure;  development  of  new  types  of  experimental  instrumentation  and  facilities;  and 
characterization  of  the  unsteady  flow  structure  on  delta  wings. 

In  the  area  of  computer-aided  techniques  of  quantitative  flow  visualization  of 
vortex  structure,  emphasis  has  been  on  generation  of  new  methods  of  two  and  three- 
dimensional  image  construction.  Such  images  are  produced  by  tracking  hydrogen 
bubble  timelines  in  three-dimensional  space  or  by  obtaining  a  large  number  of  particle 
images  in  a  given  plane  or  contiguous  planes  of  the  three-dimensional  flow. 

In  the  area  of  experimental  instrumentation  and  facilities,  novel  types  of 
integrated,  active  control  systems  have  been  developed  for  simultaneously  forcing  delta 
wings  and  quantitatively  determining  the  instantaneous  flow  structure.  Development  of 
new  quantitative  visualization  techniques  has  involved  flying  hydrogen  bubble  wires, 
direct  particle  tracking,  and  particle  image  velocimetry  employing  laser  scanning  and 
pulsed-laser  techniques. 


Characterization  of  the  unsteady  flow  structure  on  delta  wings  subjected  to 
various  forms  of  pitching  motion  addresses  the  substantial  time  shift  (or  phase  lag) 
between  the  response  of  location  of  vortex  breakdown  and  the  wing  motion,  even  at 
very  low  values  of  reduced  pitching  rate.  Moreover,  the  instantaneous  structure  of  the 
cross-section  of  the  leading-edge  vortices  during  the  pitching  motion  has  been  addressed. 
Preliminary  attempts  to  control  the  vortex  structure  using  local  and  global  control 
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COMPUTER-AIDED  TECHNIQUES  FOR  QUANTITATIVE 
INTERPRETATION  OF  FLOW  STRUCTURE 


COMPUTER-AIDED  RECONSTRUCTION  OF  THREE-DIMENSIONAL 
FLOWS  BY  LAGRANGIAN  TRACKING  OF  TIMELINES  (Lawson,  Magness) 

OBJECTIVES 

yj  Construct  global,  isometric  views  of  flow  structure  on  an  oscillating  delta 
wing  as  a  function  of  time-dependent  angle  of  attack. 

y/  Technique:  generate,  via  the  pulsed  wire  method,  timelines  of  bubbles. 
Track  these  bubble  lines  in  two  orthogonal  views  to  allow  construction  of 
instantaneous  three-dimensional  views. 


PRINCIPAL  FINDINGS/PROGRESS 

sj  A  new  CAD  technique  has  been  developed  for  determining  nree- 
dimensional  curves  from  two-dimensional  line  projections.  This  capability 
involves  techniques  such  as  intersection  of  projected  surfaces  generated 
from  two-dimensional  lines.  The  resultant  three-dimensional  lines  then  can 
be  connected  using  a  bicubic  spline  technique  to  form  a  three-dimensional 
surface. 

yj  Use  of  this  technique  has  revealed  that: 

•  At  given  angle  of  attack,  the  rate  of  development  of  the  vortical 
structure  on  the  wing  and  the  rate  of  entrainment  of  the  surrounding 
irrotational  fluid  are  highly  dependent  on  the  reduced  frequency  of 
oscillation. 

•  At  higher  reduced  frequency,  vortical  distortions  of  different  form 
occur  along  the  surface  of  the  wing  at  a  given  instantaneous  angle  of 
attack;  the  timeline  surfaces  show  regions  of  both  positive  and 
negative  curvature  over  the  extent  of  the  wing. 

•  At  low  reduced  frequency,  it  appears  that  the  three-dimensional  shape 
of  the  vortical  surface  tends  toward  conical  similarity;  at  a  threshold 
(higher)  value  of  reduced  frequency,  there  is  departure  from  this  type 
of  similarity  along  with  severe  chordwise  distortions  of  the  vortical 
surface. 
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II. 


COMPUTER-AIDED  ANALYSIS  OF  THREE-DIMENSIONAL  EVOLUTION 
OF  VORTEX  FLOW  ON  OSCILLATING  DELTA  WING:  APPLICATION 
OF  A  NEW  LASER-INDUCED  REFLECTION  TECHNIQUE  (Magness, 
Utsch,  Curtis ,  Lawson) 

OBJECTIVES 

V  Using  a  laser-scanning  technique  in  conjunction  with  the  generation  of 
multiple  bubble  surfaces,  determine  the  three-dimensional,  time-dependent 
evolution  of  the  flow  structure  on  oscillating,  three-dimensional  wings. 

V  Proposed  techniques: 

•  Generate  laser  sheets  to  define  the  desired  cross-section  of  the 
unsteady  vortical  motion; 

•  Determine,  in  the  plane  of  the  laser  sheet,  the  distortion  of  the 
multiple  bubble  surfaces; 

•  Track,  in  non-real-time  at  a  defined  angle  of  attack,  the  distortion  of  a 
defined  bubble  surface  extending  from  the  apex  of  the  wing  to  its 
trailing-end.  In  doing  so,  use  the  instantaneous  position  of  the  wing  as 
a  phase  reference. 

•  Construct  the  three-dimensional  surfaces  representing  the 
instantaneous  flow  distortion  at  a  given  angle  of  attack  using  a  newly 
developed  CAD  technique. 


PRINCIPAL  FINDINGS/PROGRESS 

V  An  experimental  technique  has  been  developed  to  generate  multiple 
bubble-timeline  surfaces.  A  three-dimensional  adjustable  laser-scanning 
system  has  also  been  perfected. 

yj  A  CAD  method  has  been  implemented  to  form  three-dimensional  surfaces 
from  the  individual  planes  of  the  flow  visualized  by  laser  sheets. 


III.  DEVELOPMENT  OF  NEW  METHODS  OF  THREE-DIMENSIONAL  IMAGE 

CONSTRUCTION  FROM  PLANES  OF  PARTICLE  IMAGES  (Robinson) 

OBJECTIVES 

yj  Construct  three-dimensional  images  from  velocity  fields  in  arbitrarily- 
spaced  amd  -oriented  plames  in  an  unsteady  flow  field. 

V  Apply  the  technique  to  unsteady  flow  paist  a  delta  wing  using  a  phase- 
referencing  technique  to  construct  the  three-dimensional  vortex  structure 
at  successive  instants  of  time. 
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PRINCIPAL  FINDINGS/PROGRESS 

This  technique  involves  two  principal  phases:  development  of  accurate 
interpolation  techniques  for  the  velocity  fields  within  the  planes  of  acquired 
data;  and  determination  of  the  out-of-plane  velocity  field  and  the  associated 
three-dimensional  streamline  pattern  fields  within  the  planes  of  acquired  data; 
and  determination  of  the  out-of-plane  velocity  field  and  the  associated  three- 
dimensional  streamline  pattern. 


sj  For  interpolation  of  the  velocity  field  within  each  plane,  two  techniques 
have  been  developed:  an  adaptive  Gaussian  window  (AGW)  method;  and 
a  spline,  thin  shell  (STS)  technique.  Both  of  these  techniques  have  been 
tested  for  specified  flows,  and  have  been  shown  to  give  results  that  axe  in 
good  agreement. 

sj  The  technique  for  determining  the  three- dimensioned  field  is  as  follows: 

•  Applying  the  continuity  equation  at  each  node  of  three-dimensional 
grid,  the  derivatives  in  the  x  and  y  directions  of  the  in-plane  (x,y 
plane)  velocity  components  (u,v)  are  fitted  by  a  finite  cosine  series. 
The  integrated  form  of  this  continuity  equation  is  approximated  by  a 
modified  finite  sine  series. 

•  Applying  inverse  discrete  cosine  and  sine  transforms,  using  a  Fast- 
Fourier  Transform  technique,  the  values  of  the  out-of-plane  velocity 
components  at  each  node  are  computed. 

•  With  the  three-dimensional  velocity  field  at  hand,  a  tri-cubic  spline 
interpolation  is  employed  to  obtain  the  three-dimensional  streamline 
pattern. 


IV.  INVERSE  APPROACHES  FOR  ASSESSING  THREE-DIMENSIONAL 

CONSTRUCTION  USING  THEORETICAL  SOLUTIONS  (Robinson) 

OBJECTIVES 

V  Using  known  theoretical  solutions  of  three-dimensional  flows,  simulate 
particle  tracking  occurring  in  laboratory  flows,  and  reconstruct  the  known 
three-dimensional  velocity  field  from  the  velocity  vectors  corresponding  to 
the  simulated  particle  images. 

sj  Feedback  inadequacies  of  the  reconstruction  technique  to  the 
computational  scheme  in  order  to  optimize  the  accuracy  and  speed  of  the 
image  construction  process. 
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PRINCIPAL  FINDINGS/PROGRESS 

y  Two  classes  of  theoretical  solutions  have  been  considered:  a  columnar 
vortex  which  contains  the  basic  elements  of  a  leading-edge  vortex  on  a 
delta  wing;  and  a  Hill’s  spherical  vortex  having  highly  distributed  vorticity 
and  representing  many  of  the  features  of  the  bubble  mode  of  vortex 
breakdown  of  a  leading-edge  vortex. 

v/  For  the  case  of  the  columnar  vortex,  the  streamlines  are  not  closed,  and 
the  flow  structure  is  invariant  with  the  x  direction.  As  a  consequence,  it  is 
possible  to  reconstruct  the  principal  features  of  this  vortex  with  minimum 
information.  Near  the  central  portion  of  the  vortex,  where  the  wavelength 
of  the  spiral  motion  is  minimum,  it  is  possible  to  obtain  a  faithful 
representation  of  the  streamline  pattern  for  a  distance  between  data  planes 
corresponding  a  half  wavelength  of  the  spiral  of  the  vortex. 

y  For  the  case  of  Hill’s  vortex,  the  streamlines  are  closed  and  the  flow 
structure  is  variant  with  streamwise  distance.  In  this  case,  it  is  necessary 
to  employ  a  sufficiently  large  number  of  planes  (approximately  twenty 
evenly  distributed  along  the  axis  of  the  vortex)  in  order  to  accurately 
reproduce  the  inner  portions  of  the  closed  streamline  pattern. 


V.  THREE-DIMENSIONAL  IMAGE  CONSTRUCTION  OF  LEADING-EDGE 

VORTICES  FROM  PARTICLE  IMAGE  DATA  (Robinson,  Magness) 

OBJECTIVES 

y  Determine  the  velocity  distribution  over  defined  planes  of  arbitrary  spacing 
and  location  along  the  leading-edge  of  a  delta  wing.  Develop  a  technique 
for  determining  the  separation  line  from  the  leading-edge  of  the  wing  from 
the  raw  particle  velocity  data. 

y  Determine  the  three-dimensional  streamline  pattern  corresponding  to  a 
leading-edge  vortex  on  a  stationary  delta  wing. 

y  Using  phase- referencing  techniques,  extend  the  foregoing  procedures  to  the 
unsteady  leading-edge  vortex  on  a  pitching  delta  wing. 

PRINCIPAL  FINDINGS/PROGRESS 

y  A  method  has  been  developed  for  determining  the  velocity  field  over  the 
cross-section  of  the  leading-edge  vortex.  The  separation  streamline  from 
the  edge  of  the  wing  is  evaluated  by  application  of  either  the  adaptive 
Gaussian  window  (AGW)  or  spline  thin  shell  (STS)  interpolation  technique 
to  the  raw  particle  image  data.  A  critical  assessment  of  the  accuracy  of 
the  two  techniques  shows  that  they  give  very  close  agreement  for  the 
separation  streamline.  The  separation  streamline  is  interpreted  as  a  line 
cutting  across  the  three-dimensional  streamsurface  of  leading-edge 
vortex,  coincident  with  the  plane  of  the  original  particle  image  data.  In 
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the  event  that  vortex  breakdown  occurs  within  the  core  of  the  leading-edge 
vortex,  interpolation  is  not  valid  within  the  breakdown  region.  The 
boundary  of  the  (noncircular)  domain  of  breakdown  is  determined  by  the 
interface  between  the  random  and  organized  velocity  vector  orientations  of 
the  original  particle  image  data.  This  technique  has  been  applied  to  both  a 
leading-edge  vortex  from  a  stationary  wing  and  its  corresponding 
counterpart  from  a  wing  undergoing  pitching  motion. 

sj  The  velocity  fields  in  seven  planes  covering  a  chordwise  domain  of 
approximately  two-thirds  the  chord  of  a  delta  wing  have  been  employed  to 
determine  the  three-dimensional  streamline  pattern.  In  order  to 
accomplish  this,  the  velocity  field  in  a  plane  angled  with  respect  to  the 
freestream  is  used  to  initiate  the  out-of-plane  construction  process.  Use  of 
the  three-dimensional  image  construction  method  allows  determination  of 
the  chordwise  streamline  pattern. 


VI.  TOPOLOGICAL  METHODS  FOR  INTERPRETATION  OF  FLOW 
STRUCTURE  ON  A  PITCHING  DELTA  WING 


OBJECTIVES 

yj  Develop  a  streamline  construction  technique,  to  be  applied  to 
instantaneous  velocity  fields,  that  allows  interpretation  of  the  primary 
characteristics  of  the  instantaneous  flow  structure  on  the  basis  of  critical 
point  theory. 

n/  With  the  critical  point  theory  concepts  at  hand,  and  use  of  a  previously- 
developed  three-dimensional  construction  technique,  develop  methods  for 
constructing  the  three-dimensional  and  streamline  and  streamsurface 
patterns  along  the  chord  of  a  pitching  delta  wing. 


PRINCIPAL  FINDINGS/PROGRESS 

U  An  interactive  process  has  been  developed  for  determining  the  saddle 
points  and  foci  over  am  entire  plan  of  data  at  a  three-dimensional  flow. 
This  method  adlows  construction  of  the  instantameous  streamline  pattern  or 
any  cross-section  and  has  particular  merit  for  interactively  determining  the 
locations  of  saddle  points  amd  constructing  the  intersecting  streamline 
patterns  in  the  vicinity  of  them.  This  software  allows  ready  identification 
of  the  foci  (vortex  centers)  and  whether  they  exhibit  stable  or  unstable 
characteristics. 

s/  A  method  has  been  developed  for  construction  of  the  three-dimensional 
streamline  patterns  amd  streamsurfaces  at  any  instant  (i.e.  given  angle  of 
attatck)  along  the  entire  chord  of  the  delta  wing.  This  method  takes 
advantage  of  the  previously  developed  techniques  involving  spectral  and 
finite  difference  methods  to  determine  the  out-of-plame  velocity  component, 
amd  has  a  number  of  new,  additionad  features  that  allow  for  proper 
smoothing  of  streamline  patterns  in  three-dimensional  space. 
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VII.  INTERA CTIVE  FLOW  VISUALIZA  TION- THEORETICAL  SIMULA  TION 

(Gursul) 


OBJECTIVES 

V/  Determine  the  relationship  between  visualized  flow  patterns  of  unsteady 
shear  flows  and  the  underlying  vorticity  fields. 


s/  Formulate  criteria  for  pattern  recognition  schemes  in  two-  and  three- 
dimensional  unsteady  flows  and  characterize  fallacies  in  interpreting 
visualization  of  these  flows. 


PRINCIPAL  FINDINGS/PROGRESS 

•J  The  fundamental  case  of  an  unsteady  free  sheax  flow  in  the  form  of  a  row 
of  simple  vortices  having  the  same  orientation  has  been  considered  using 
Stuart's  exact  nonlinear  solution.  This  formulation  allows  specification  of 
arbitrary  concentration  of  vorticity  for  a  given  value  of  circulation.  The 
evolution  of  streaklines  along  the  critical  layer  has  been  considered  for 
severed  values  of  the  vorticity  concentration  parameter. 

v/  An  increased  concentration  of  vorticity  actually  results  in  a  larger  cross¬ 
stream  extent  of  the  visualized  streakline  pattern  for  a  given  value  of 
circulation. 

sj  Essentially  identical  streakline  patterns  can  be  obtained  for  drastically 
different  combinations  of  vorticity  concentration  parameter  and  circulation. 
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EXPERIMENTAL  INSTRUMENTATION  AND  FACILITIES: 
NEW  TECHNIQUES 


VIII.  INTEGRATED  ACTIVE  CONTROL/QUANTITATIVE  FLOW 
VISUALIZATION/FORCE  MEASUREMENT  (Magness,  Frednksson) 

OBJECTIVES 

V  Develop  a  computer-controlled  system  that  allows  simultaneous  control  of: 
the  pitching  motion  of  a  delta  wing;  flow  visualization  using  laser-scanning 
and  pulsed-laser  illumination  of  particles;  and  measurement  of  the 
unsteady  forces  on  the  oscillating  delta  wing. 

yj  Develop  a  phase-referencing  technique  that  will  allow  acquisition  of 
velocity  and  force  information  that  cannot  be  secured  simultaneously.  This 
phase  referencing  is  to  be  specified  by  the  digital  format  of  the  motor 
controller  system. 


PRINCIPAL  FINDINGS/PROGRESS 

V  A  new,  large-scale  apparatus  for  pitching  delta  wings  has  been  designed 
and  constructed.  It  allows,  in  principle,  arbitrary  pitching  motion  and 
variation  of  the  pitching  axis.  It  is  controlled  by  a  Compumotor  controller 
system,  which  is  driven  by  the  central  microcomputer.  Corresponding 
software  in  “C”  language  has  been  developed  in  order  to  allow  forcing  of 
the  pitching  wing  in  various  ramp-type  maneuvers. 

V  A  laser  scanning  system  for  illumination  of  particles  allows  quantitative 
visualization.  The  camera  recording  the  visualization  is  fired  by  the 
microcomputer-controlled  system,  allowing  images  to  be  acquired  at 
desired  phases  of  the  pitching  motion  of  the  wing. 

yj  A  high  sensitivity  force  measurement  system  has  been  developed  using  high 
sensitivity  strain  gauges.  Force  signals  are  in  the  form  of  modulated 
carrier  signals  in  order  to  improve  accuracy. 


IX.  DIRECT  PARTICLE  TRACKING  VIA  LASER  SCANNING  (Magness,  Kuo) 

OBJECTIVES 

V  Develop  a  technique  for  illuminating  suspended  particles.  Requirements 
are:  an  intensity  of  illumination  an  order  of  magnitude  higher  than  that 
attainable  with  conventional  chopper  (rotating  disk)  arrangements;  and  a 
variable  frequency  of  illumination. 

yj  Obtain  a  sufficiently  high  concentration  of  illuminated  particle  images  and 
values  of  displacement  of  particle  images  that  will  allow  accurate 
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characterization  of  the  flow  field  by  direct  determination  of  displacements 
of  successive  particle  images. 


PRINCIPAL  FINDINGS/PROGRESS 

v/  A  galvanometer-driven  rotating  mirror  system,  externally  controlled  by  a 
function  generator,  undergoes  successive  ramp- type  oscillations  at  rates  as 
high  as  500  c/s.  The  intensity  and  duration  of  the  scanning  rate  is 
sufficiently  high  to  allow  generation  of  as  many  as  three  to  five  multiple 
exposures  of  a  typical  particle  image  during  displacement  on  the  film 
negative.  These  multiple  exposures  allow  verification  of  the  in-plane 
motion  of  each  particle. 

sj  As  many  as  three  hundred  to  five  hundred  particle  trajectories  are 
typically  acquired  over  a  given  cross-section  of  the  flow  using  Pliolite 
particles  having  characteristic  diameters  of  the  order  of  50  to  75  microns. 
These  successive  displacements  of  a  given  particle  are  digitized  using  a 
Tektronix  digitizing  tablet  from  very  large  prints  of  the  particle  images. 

s/  Determination  of  the  randomly  located  velocity  vectors,  followed  by 
application  of  either  an  adaptive  Gaussian  window  (AGW)  or  spline  thin 
shell  (STS)  interpolation  scheme  provides  the  instantaneous  velocity  field 
over  a  cross-section  of  the  flow. 

v/  The  unsteady  flow  structure  on  the  delta  wing  can  be  selectively  sampled 
using  computer-generated  firing  pulses  transmitted  to  the  camera. 

J  For  both  stationary  and  unsteady  flows,  it  is  possible  to  construct  the 
instantaneous  streamline  pattern.  Moreover,  the  separation  streamlines 
can  be  determined  by  using  a  tracking  algorithm. 


PARTICLE  IMAGE  VELOCIMETRY  VIA  LASER  SCANNING  AND 

PULSED  LASER  TECHNIQUES  (Magntss,  Kuo,  and  Konak) 

OBJECTIVES 

sj  Aim  for  high  spatial  resolution,  instantaneous  velocity  fields  over  a  defined 
cross-section  of  the  flow  field.  Use  small  particles  (two  to  ten  micron 
diameter),  in  conjunction  with  illumination  of  high-powered  YAG  lasers 
having  a  pulse  width  of  the  order  of  nanoseconds,  in  order  to  ensure  that 
the  particles  do  not  drift  on  the  film  negative  during  successive 
illumination  pulses. 

s/  Establish  a  means  of  interrogating  the  film  negative  showing  the  successive 
locations  of  the  particle  images  due  to  multiple  exposure  of  the  film. 
Employ  a  Young’s  Fringe  Technique  involving  an  optical  Fourier 
Transform  Technique,  in  order  to  determine  the  magnitude  and  direction 
of  the  velocity  at  a  given  location  of  the  film  negative. 


v/  Develop  data  evaluation  techniques  that  will  allow  evaluation, 
transformation  and  display  of  the  large  data  sets  (approximately  four  to 
five  thousand  velocity  vectors  per  cross-section  of  the  flow),  including 
methods  for  determining  the  instantaneous  streamline  patterns  and 
contours  of  constant  vorticity. 


PRINCIPAL  FINDINGS/PROGRESS 

s/  A  new  water  channel  system,  made  entirely  of  plastic  so  that  foreign 
particles  will  not  contaminate  the  marker  particles,  has  nearly  been 
completed.  It  has  fine  filtration  systems  for  properly  conditioning  the 
water. 

sj  The  first  goal  has  been  to  determine  a  method  for  effectively  combining  the 
beams  from  two  YAG  lasers  into  a  single  colinear  beam.  Design  of  this 
aspect  has  been  completed,  including  use  of  a  beamsplitter  cube  especially 
fabricated  for  high-powered  laser  applications. 

s/  An  optical  method,  unique  to  our  laboratory,  for  forming  a  laser  sheet  from 
a  cylindrical  beam  cross-section,  and  translating  it  to  desired  locations  in 
the  flow,  has  been  designed.  It  involves  use  of  beam  steering  optics  in 
conjunction  with  spherical  and  cylindrical  lenses,  all  of  which  must  be 
specially  selected  to  account  for  the  high  laser  power  of  the  order  of  one 
Joule/cm  . 
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FLOW  STRUCTURE  AND  LOADING;  PHYSICAL  FINDINGS 


XI.  FLOW  STR UCTURE  ON  OSCILLA  TING  WING  SEGMENTS  (Utsch) 
OBJECTIVES 

yj  Determine  basic  classes  of  three-dimensional  vortex  generation  from 
interacting  wing  segments  with  application  to:  wing-flap  configurations; 
and  adjustable  wings.  Investigate  new  concepts  of  lift  generation 
associated  with  controlled  motion  of  wings  and  wing  segments. 

•J  Examine  effects  of  relative  frequency,  phase,  and  amplitude  of  oscillating 
wing  segments  on  three-dimensional  flow  structure.  This  characterization 
of  the  flow  structure  will  include:  trajectories  of  the  generated  three- 
dimensional  vortices;  characteristic  scale  and  circulation  of  vortical 
structures;  and  definition  of  vortex  interference  phenomena. 


PRINCIPAL  FINDINGS/PROGRESS 

yj  Vortices  can  coexist  on  both  the  interior  and  exterior  surfaces  of  an 
oscillating  wing  segment. 

y/  At  a  given  value  of  streamwise  (chordwise)  location,  the  spanwise  position 
of  the  interior  or  exterior  vortex  remains  nearly  invariant.  This  invariance 
persists  despite  large  excursions  in  the  angle  between  the  oscillating  wing 
segments,  including  the  limiting  case  where  the  wing  segments  come  into 
contact  during  the  oscillation  cycle. 

y  At  a  threshold  value  of  dimensionless  oscillation  frequency,  the  interior 
region  of  the  wing  segments  takes  on  a  fully  turbulent  character.  At  this 
same  frequency,  large-scale  vortices  are  generated  at  closure  of  the  wing 
segment  system. 


XII.  FLOW  STRUCTURE  ON  PITCHING  DELTA  WING  BY  USE  OF  LASER¬ 
SCANNING  TECHNIQUES  (Magness) 


OBJECTIVES 

yj  Determine  the  character  of  the  vortical  flow  structure  as  a  function  of 
mean  angle  of  attack,  reduced  frequency,  and  sweep  angle  of  the  delta 
wing. 

yj  Define  the  hysteresis  of  the  three-dimensional  vortex  formation  as  a 
function  of  streamwise  (chordwise)  distance  from  the  apex  of  the  airfoil  by 
tracking  the  streamwise  development  of  the  vortex  from  the  apex  of  the 
wing  as  a  function  of  pitch  angle  a(t). 
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\J  Investigate  the  extent  to  which  the  vortex  formation  process  can  be 
suppressed  for  certain  types  of  excitation;  consider  whether  this  so-called 
suppression  occurs  over  the  entire  wing  or  over  a  limited  range  of 
streamwise  (chordwise)  distance  along  the  wing. 


PRINCIPAL  FINDINGS/PROGRESS 

>/  Hysteresis  of  the  vortex  formation  relative  to  the  wing  motion  extends  to 
relatively  low  values  of  reduced  frequency.  Moreover,  this  hysteresis 
process  emanates  from  the  apex  region  of  the  wing  and  is  evident  at  all 
locations  along  the  wing  from  the  apex  to  the  trailing-edge. 

sj  At  sufficiently  high  reduced  frequency,  and  over  a  certain  interval  of 
chordwise  distance  along  the  wing,  there  is  no  vortex  on  either  the  upper  or 
lower  side  of  the  wing,  despite  the  fact  that  such  vortex  formation  does 
occur  on  the  corresponding  stationary  wing. 

yj  At  sufficiently  high  reduced  frequency,  and  over  a  certain  interval  or 
chordwise  distance  along  the  wing,  there  is  no  vortex  on  either  the  upper  or 
lower  side  of  the  wing,  despite  the  fact  that  such  vortex  formation  does 
occur  on  corresponding  stationary  wing. 

y/  Strongly  distorted  (suppressed)  vortex  patterns  are  evident  as  the  trailing- 
edge  is  approached  when  the  wing  is  excited  at  sufficiently  high  reduced 
frequency.  This  distortion  is  due  to  the  integrated  effect  of  the  upstream 
development  of  the  flow. 


XIII.  FLOW  STRUCTURE  ON  A  PITCHING  DELTA  WING  BY  LAGRANGIAN 
TRACKING  TECHNIQUES  EMPLOYING  A  FLYING  WIRE  (Atta) 

OBJECTIVES 

yj  Determine  the  distortion  of  the  flow  structure  entering  the  plane  of  the 
wing  at  its  apex;  in  doing  so,  employ  a  wire  oriented  in  the  plane  of  the 
wing  and  located  close  to  the  apex  of  the  wing.  This  wire  is  to  “fly”  at  the 
same  velocity  as  the  apex  of  the  wing. 

yj  Relate  the  local  flow  structure  at  the  apex  of  the  wing  to  that  along  the 
entire  wing,  extending  into  the  near-wake  region.  Characterize  the  basic 
types  of  three-dimensional  vortex  formation  in  the  near-wake  of  the  delta 
wing. 

y/  Interpret  the  hysteresis  of  the  initial  stages  of  vortex  formation  near  the 
apex  to  the  manner  in  which  flow  is  swept  above  and  below  the  apex  of  the 
wing  during  the  oscillation  cycle. 

y/  Characterize  the  above- described  aspects  of  the  flow  structure  in  terms  of 
mean  angle  of  attack  and  reduced  frequency. 
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PRINCIPAL  FINDINGS/PROGRESS 

\J  Due  to  dynamic  effects,  the  flow  entering  the  plane  of  the  wing  at  the  apex 
can  pass  below,  in  contrast  to  above,  the  wing  surface;  this  effect  occurs 
even  at  relatively  high  positive  values  of  angle  of  attack  a. 

\J  Flow  entering  the  plane  of  the  wing  and  passing  by  the  apex  shows 
substantial  hysteresis,  even  at  the  lower  values  of  reduced  frequency 
examined  herein.  These  hysteresis  patterns  can  be  characterized  by  an 
effective  sweep  angle  of  the  timeline  pattern  in  the  region  of  the  apex. 
This  effective  sweep  angle  of  the  flow  structure,  relative  to  the  geometrical 
sweep  angle  of  the  wing,  shows  substantial  hysteresis. 

yj  At  a  threshold  (relatively  high)  value  of  reduced  frequency,  at  any  instant 
during  the  wing  motion,  the  flow  passes  above  the  wing  over  a  certain 
portion  of  the  wing,  and  below  the  wing  over  another  portion  of  the  wing. 

V  The  foregoing  types  of  flow  distortion  at  the  apex  and  along  the  entire 
surface  of  the  wing  lead  to  several  different  classes  of  three-dimensional 
vortex  formation  from  the  trailing-edge  of  the  wing. 


XIV.  STRUCTURE  OF  VORTEX  CORE  DEVELOPMENT  AND 
BREAKDOWN  IN  A  PITCHING  DELTA  WING  (Atta) 

OBJECTIVES 

v/  Determine  the  mechanism  of  formation  of  the  vortex  core  on  an  oscillating 
delta  wing  as  a  function  of  angle  of  attack  and  reduced  frequency. 

V  Identify  the  classes  of  instability  and  breakdown  of  the  vortex  core  relative 
to  the  vortex  core  on  a  static  wing.  Quantify  these  observations  by 
determining  the  axial  velocity  of  the  centerline  of  the  vortex  core. 

V  Define  the  hysteresis  of  vortex  core  breakdown,  relative  to  the  motion  of 
the  wing. 

PRINCIPAL  FINDINGS/PROGRESS 

yj  At  a  threshold  value  of  reduced  frequency,  there  is  no  vortex  core  along  the 
surface  of  the  wing  over  a  significant  portion  of  the  oscillation  cycle. 

V  The  onset  of  vortex  core  development  is  distinctly  different  at  low  reduced 
frequency  and  high  reduced  frequency.  At  low  reduced  frequency,  there  is 
simultaneous  development  of  the  entire  column  of  the  vortex  core  along  its 
entire  stream  wise  extent.  At  higher  reduced  frequency,  there  is  rapid 
ejection  of  the  leading-edge  of  the  core  from  the  apex  region  of  the  wing. 

yj  Irrespective  of  the  mechanism  of  development  of  the  vortex  core,  the  axial 
velocity  of  the  core  remains  constant  during  the  development  process. 
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y  The  eventual  breakdown  of  the  core  of  the  vortex  can  be  grouped  into  four 
basic  classes,  ranging  from  a  spiral-like  mode  similar  to  that  appearing  on  a 
static  wing  to  a  rapid,  column- type  collapse  of  the  entire  core. 

y  A  secondary  vortex,  along  the  leading-edge  of  the  wing,  can  develop 
simultaneously  with  the  primary  vortex.  Its  core  displays  the  same 
principal  features  of  development  and  breakdown  as  that  of  the  primary 
vortex. 


XV.  RESPONSE  OF  VORTEX  BREAKDOWN  TO  TYPE  OF  PITCHING 

MANEUVER  (Magnesa) 

OBJECTIVES 

y  Implement  active  control  techniques  that  allow  maneuver  of  the  delta  wing 
in  the  pure  pitching  mode  at  a  desired  pitching  rate  for  several  types  of 
maneuver:  simple  pitch-up;  simple  pitch-down;  and  continuous  pitch-up 
and  pitch-down  motions. 

y  Observe,  with  the  use  of  hydrogen  bubble  and  dye  injection  techniques,  the 
location  of  the  occurrence  of  vortex  breakdown  with  respect  to  the 
instantaneous  position  of  the  wing  during  the  course  of  a  given  maneuver, 
in  order  to  provide  a  basis  for  studying  the  detailed  structure  of  the  vortex 
while  the  wing  is  in  motion. 


PRINCIPAL  FINDINGS/PROGRESS 

y  Dye  injection  and  hydrogen  bubble  techniques  have  been  developed  to 
allow  movement  of  the  location  of  marker  injection  with  the  motion  of  the 
wing.  Notable  is  the  technique  for  manufacturing  the  segmented  hydrogen 
bubble  wire  and  its  installation  along  the  leading-edge  of  the  wing.  Both 
high  speed  video  and  35  mm  camera  techniques  ha,re  been  employed  to 
characterize  the  vortex  breakdown  during  the  maneuver. 

y  Some  of  the  principal  features  of  the  movement  of  the  vortex  breakdown 
location  with  respect  to  the  wing  motion  are  as  follows: 

•  Irrespective  of  whether  one  considers  pitch-up  or  pitch-down  motion, 
the  axial  location  has  its  lowest  value  when  the  location  of  vortex 
breakdown  is  at  its  furthest  upstream  location,  coinciding  with  the 
existence  of  large-scale  vortex  breakdown  over  a  substantial  portion  of 
the  wing  surface.  Higher  velocities  tend  to  occur  when  the  breakdown 
is  on  the  downstream  portion  of  the  wing. 

•  When  the  continuous  pitch-up-down  motion  of  the  wing  is  invoked, 
the  breakdown  location  remains  downstream  of  the  midchord  and  low 
values  of  the  axial  velocity  of  the  breakdown  location  are  avoided. 
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XVI.  RELATION  OF  INSTANTANEOUS  FORCE  TO  TYPE  OF  PITCHING 

MANEUVER  (Fredriksson) 

OBJECTIVES 

yj  Design  and  implement  a  high  sensitivity  force  measurement  system  for 
characterizing  the  unsteady  Bft  and  moment  on  a  delta  wing  undergoing 
unsteady  pitching  motion  in  water. 

V  Determine  the  variation  of  lift  as  a  function  of  instantaneous  angle  of 
attack  for  various  values  of  pitching  rate  and  initial  and  terminal  angles  of 
attack  during  ramp-type  maneuvers.  If  possible,  extend  these  same 
characterizations  to  unsteady  moment  on  the  wing. 

x/  Compare  the  variations  of  lift  and  vortex  breakdown  position  with  angle  of 
attack  to  determine  the  correlation  between  them. 


PRINCIPAL  FINDINGS/PROGRESS 

s/  A  specially-designed  force  measurement  system  for  use  in  water  has  been 
completed  and  implemented.  The  strain  gauge  system  is  integrated  with  a 
signal  processing  system  that  employs  a  modulated-carrier  concept  in 
contrast  to  traditional  DC  measurement  of  forces.  Data  acquisition  and 
reduction  software  have  been  developed  in  order  to  allow  characterization 
of  the  instantaneous  forces  as  a  function  of  time  or  angle  of  attack. 
Moreover,  method  for  ensemble-averaging  successive  force  measurements 
has  been  developed. 

V  The  instantaneous  lift  vs.  angle  of  attack  has  been  determined  for  several 
values  of  pitching  rate.  Characteristics  of  the  dynamic  lift  such  as 
overshoot  of  the  static  characteristic  and  hysteresis  of  the  lift  vs.  angle  of 
attack  have  been  determined.  Increasing  the  value  of  pitching  rate 
substantially  increases  the  area  of  the  hysteresis  loop  via  maintain  low 
values  of  lift  during  the  pitch-down  portion  of  the  pitching  motion. 
However,  the  overshoot  above  the  static  lift  characteristic  does  not  seem  to 
be  significantly  influenced  by  pitching  rate. 

y/  Correlation  of  the  instantaneous  location  of  vortex  breakdown  with  the 
instantaneous  lift  on  the  wing  suggests,  contrary  to  the  assumed  behavior, 
that  there  is  not  a  direct  relationship  between  movement  of  the  vortex 
breakdown  position  along  the  surface  of  the  wing  and  changes  in  lift.  This 
aspect  is  currently  under  investigation. 
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XVII.  INSTANTANEOUS  FLOW  STRUCTURE  ACROSS  LEADING-EDGE 
VORTEX  IN  RELATION  TO  TYPE  OF  PITCHING  MANEUVER  (Magness) 


OBJECTIVES 

\J  Determine  the  detailed  flow  structure  at  selected  cross-sections  on  a 
stationary  delta  wing,  including  the  detailed  velocity  distribution  across  the 
specified  cross-section  as  well  as  the  separation  line  emanating  from  the 
leading-edge  of  the  wing. 

y/  Obtain  similar  information  for  the  case  of  a  wing  undergoing  unsteady 
pitching  motion,  and  directly  compare  with  the  corresponding  stationary 
case,  in  order  to  determine  the  underlying  physical  basis  for  the  time  shift 
of  the  response  of  the  vortex  structure  with  respect  to  the  motion  of  the 
wing.  Use  this  information  to  aid  in  interpretation  of  the  nature  of  the 
chordwise  movement  of  the  vortex  breakdown  location  during  a  typical 
maneuver. 


PRINCIPAL  FINDINGS/PROGRESS 

\J  A  laser  scanning  system,  which  allows  movement  of  the  laser  sheet  to 
arbitrary  chordwise  location  and  its  orientation  with  respect  to  the 
freestream,  has  been  implemented  in  conjunction  with  the  active  controller 
system  with  the  wing  motion.  During  this  first  phase  of  the  investigation, 
efforts  have  focussed  on  the  characterization  of  the  structure  of  the  vortex 
at  the  midchord  of  the  wing,  coincident  with  the  location  of  the  pitching 
axis.  Direct  comparisons  of  the  detailed  structure  of  the  leading-edge 
vortex  have  been  made  for  the  pitch-down  (ramp-type)  motion  and 
compared  directly  with  the  structure  of  the  vortex  existing  on  a  stationary 
wing  at  the  same  angle  of  attack. 

yj  Preliminary  observations  of  the  instantaneous  structure  of  the  vortex  cross- 
section  on  the  pitching  wing,  in  comparison  with  the  stationary  wing,  are: 

•  In  absence  of  any  breakdown  within  its  core,  the  detailed  structure  of 
the  vortex  cross-section  relaxes  to  its  steady  state  value  within 
approximately  one  convective  timescale  C/Uqo  after  cessation  of  the 
wing  motion. 

•  The  existence  of  vortex  breakdown  within  the  core  of  the  vortex 
precludes  significant  relaxation  of  the  detailed  structure  of  the  vortex 
cross-section  within  time  scales  of  the  order  of  C/U<x>- 


16 


XVIII.  RESPONSE  OF  LEADING-EDGE  VORTEX  TO  TRANSIENT  SUCTION 
(Parmenter) 


OBJECTIVES 

sj  Determine  the  effectiveness  of  localized  suction,  simulating  a  point  sink, 
located  well  downstream  of  the  onset  of  vortex  breakdown.  Determine  the 
minimum  suction  amplitude  for  which  the  vortex  core  can  be  restabilized. 

\J  Characterize  the  transient  response  of  the  restabilization  and 
destabilization  of  the  vortex  core  due  to  abrupt  onset  and  abrupt  cessation 
of  the  suction  respectively. 

\J  Determine  the  nature  of  the  unsteady  restabilization  of  the  vortex  core 
from  the  state  of  full  breakdown  to  no  breakdown. 


PRINCIPAL  FINDINGS/PROGRESS 

yj  Experiments  have  focussed  on  the  process  of  restabilization  of  the  vortex 
core  over  a  distance  of  about  one-half  chord  of  the  wing.  Vortex 
breakdown  in  its  natural  state  occurs  at  midchord  and  the  tip  of  the 
suction  probe  is  located  approximately  at  the  trailing-edge  of  the  wing.  In 
this  case,  it  is  possible  to  stabilize  the  vortex  core  with  a  suction  coefficient 
C^  as  low  as  0.035.  Substantially  lower  values  of  suction  coefficient  should 
be  attainable  when  the  location  of  the  probe  tip  is  closer  to  the  location  of 
naturally-occurring  vortex  breakdown. 

yj  When  there  is  abrupt  cessation  of  the  suction  after  the  vortex  core  has 
been  stabilized,  there  occurs  substantial  hysteresis  of  the  vortex  breakdown 
location  Xi  versus  time  t,  relative  to  that  occurring  for  the  abrupt  onset  of 
suction.  The  magnitude  and  direction  of  this  hysteresis  depends  upon  the 
amplitude  of  the  suction  coefficient  C^. 

yj  When  the  suction  coefficient  reaches  a  sufficiently  high  amplitude  for 
the  (transient)  onset  of  suction,  the  total  time  elapsed  for  restabilization  of 
the  vortex  core  is  invariant  with  amplitude  of  the  coefficient  Cn.  On  the 
other  hand,  when  the  amplitude  is  sufficiently  small,  the  overall 
response  time  increases  rapidly  with  increasing  C^. 


XIX.  RESPONSE  OF  LEADING-EDGE  VORTEX  TO  TRANSIENT  BLOWING 

(Kuo,  Gu) 

OBJECTIVES 

yj  Determine  the  response  of  the  vortex  breakdown  to  localized  blowing 
applied  at  critical  locations  on  the  leading-edge  of  the  delta  wing;  at  the 
apex;  along,  and  in  line  with  the  leading-edge  of  the  wake;  and  along  the 
axis  of  the  vortex. 
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yj  Since  the  possibility  for  restabilization,  i.e.  downstream  movement  of  the 
vortex  breakdown  position,  for  the  case  of  blowing  applied  tangentially 
along  a  rounded  leading-edge. 

yj  For  the  forgoing  types  of  blowing,  apply  transient  blowing,  in  order  to 
characterize  the  overall  response  time  until  restabilization  of  the  vortex 
core  is  attained. 

sj  Determine  the  time-dependent  cross-sectional  structure  of  the  vortex 
during  the  process  of  the  transient  blowing  as  well  as  during  the  relaxation 
of  the  vortex  core  after  the  cessation  of  blowing. 


PRINCIPAL  FINDINGS/PROGRESS 

yj  A  large  3cale  half  delta  wing  has  been  designed  and  constructed.  It  allows 
arbitrary  angle  of  attack,  while  maintaining  the  location  and  orientation  of 
the  localized  blowing  at  a  fixed  location  and  orientation  in  space.  For  the 
case  of  distributed  blowing  in  the  tangential  direction  along  the  entire 
leading-edge  of  the  wing,  a  special  leading-edge,  variable  slit  arrangement 
has  been  designed.  For  all  types  of  blowing,  a  transient  function  generator 
has  been  developed.  It  involves  a  piston-like  arrangement,  driven  by  a 
stepping  motor,  which,  in  turn,  is  linked  to  the  laboratory  microcomputer. 
This  function  generator  allows  various  forms  of  time-dependent  blowing  to 
the  applied. 

yj  Although  all  of  the  forgoing  types  of  blowing  described  in  the  objectives  are 
effective,  and  can  produce,  to  varying  degrees,  downstream  movement  at 
the  location  of  vortex  breakdown,  the  most  effective  technique  seems  to  be 
the  application  of  tangential  blowing  from  the  rounded  leading-edge  of  the 
delta  wing.  Particularly  interesting  is  the  fact  that  when  this  tangential 
blowing  is  applied  with  an  alternating  so-called  blowing/suction  cycling,  at 
a  dimensionless  frequency  of  the  order  of  a  convective  time  scale,  effective 
restabilization  of  the  vortex  core  can  be  attained  over  a  distance  of  the 
order  of  one-half  the  chord  of  the  wing. 

yj  Corresponding  particle  tracking  of  the  cross-sectional  development  of  the 
vortex  chord  during  the  process  of  alternate  blowing  and  suction  is 
underway.  At  a  given  cross-section  of  the  leading-edge  vortex,  located 
downstream  of  the  original  breakdown  position,  a  complete  recovery  of  the 
non-broken  down  leading-edge  vortex  is  obtainable  by  the  application  of 
alternating  blowing/suction. 
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XX.  PUBLICATIONS 

The  major  thrust  of  this  research  program  has  been  determination  of  the  flow 
structure  on  swept  edges  subjected  to  controlled  motion.  Emphasis  has  been  on 
techniques  of  quantitatively  visualizing  the  flow,  while  simultaneously  subjecting  the 
wing  to  various  types  of  active  control.  Using  these  techniques,  it  is  possible  to  directly 
relate  the  flow  structure  to  the  excitation  conditions.  Several  basic  types  of  control 
have  been  implemented.  These  control  approaches  encompass: 

(a)  low  frequency,  sinusoidal  motion  of  swept  wings,  focusing  on  study  of  the  phase 
shift  between  the  evolution  of  the  leading-edge  vortex  and  the  motion  of  the  wing; 

(b)  high  frequency  perturbations  of  the  leading-edge  of  the  wing,  leading  to  alteration 
of  the  structure  of  the  separating  shear  layer  and  modification  of  the  time-mean 
structure  and  location  of  vortex  breakdown; 

(c)  large-amplitude  motion  of  segmented,  delta-wing  flaps  allowing  generation  of 
substantial  circulation  and  lift  at  zero  mean  angle  of  attack; 

(d)  ramp-type  motion  of  delta  wings,  allowing  description  of  the  very  long  time  scales 
of  the  vortex  responses;  and 

(e)  application  of  unsteady  source-like  transients  to  fully  stalled  wings,  providing 
substantial  extensions  of  the  length  of  the  vortex  core  prior  to  occurrence  of 
breakdown. 

Parallel  to  these  investigations,  formulation  and  instrumentation  of  new 
techniques  of  quantitatively  interpreting  the  flow  structure  with  use  of  computer-aided 
methods  have  been  addressed.  Concepts  of  pattern  recognition  and  three-dimensional 
construction  have  been  pursued. 

The  following  list  of  publications  provides  synopses  of  progress  to  date. 

Hysteresis  of  Vortex  Development  and  Breakdown  on  an  Oscillating  Wing n  by  R. 
Atta  and  D.  Rockwell,  AIAA  Journal.  Vol.  25,  No.  11,  1987,  pp.  1512-1518. 

A  delta  wing  subjected  to  sinusoidal  pitching  motion  exhibits  vortex 
formation  from  the  leading-edge  having  several  distinguishing  features  relative  to 
the  case  of  vortex  formation  from  a  stationary  wing.  There  is  substantial  phase 
shift  of  the  onset  of  vortex  breakdown  relative  to  the  instantaneous  position  of  the 
wing  such  that  maximum  breakdown  length  may  occur  near  or  at  the  maximum 
angle  of  attack.  Moreover,  it  is  possible  for  the  core  of  the  vortex  to  disappear  over 
a  portion  of  the  oscillation  cycle.  These  observations  occur  at  moderate  and  high 
values  of  reduced  frequency;  further  studies  are  addressing  the  vortex  behavior  at 
very  low  reduced  frequencies. 

sj  “Ensemble- Averaging  and  Correlation  Techniques  for  Flow  Visualization  Images ” 

by  P.  J.  M.  Kerstens  and  D.  Rockwell,  Experiments  ia  Fluids.  Vol.  6,  1988,  pp. 
409-419. 

Flow  visualization  using  marking  techniques  such  as  timelines  provides  a 
basis  for  quantitative  analysis  of  macroscale  features  of  unsteady  flows  by  global 
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ensemble- averaging  and  correlation  techniques.  In  the  visual-ensemble-averaging 
technique  described  herein,  the  timeline  positions  are  tracked  and  averaged  in 
successive  images.  The  phase  reference  for  the  averaging  process  can  take  the  form 
of  an  analog  pressure,  velocity,  or  displacement  signal,  or  a  recurring  coherent 
portion  of  the  image.  Global  correlations  of  the  timeline  patterns  are  obtained 
using  the  same  timelines  defined  for  the  ensemble-averaging  process.  A  new  type 
of  visual  correlation  function,  giving  the  correlation  between  two  timelines  in  a 
given  image  or  successive  images,  is  proposed. 

\J  “ Estimation  of  Velocity  Eigenfunction  and  Vorticity  Distributions  from  the 

Timeline  Visualization  Technique ”  by  D.  Lusseyran  and  D.  Rockwell,  Experiments 
in  Fluids.  Vol.  6,  1988,  pp.  228-286. 

For  the  case  of  quasi-periodic  flow,  it  is  demonstrated  that  use  of  the 
hydrogen  bubble  timeline  method  leads  to  reasonable  estimates  of  the  eigenfunction 
of  the  streamwise  velocity  fluctuation.  Both  amplitude  and  phase  distributions 
across  an  unstable  wake  flow  are  well-approximated.  It  is  shown  that  the  vorticity 
extrema,  as  well  as  the  degree  of  concentration  of  vorticity,  are  in  good  agreement 
with  those  calculated  from  linear  stability  theory. 

A  critical  assessment  is  given  of  the  possible  uncertainties  associated  with 
this  technique;  the  existence  of  a  finite,  but  unknown  cross-stream  velocity 
component;  bubble  rise  due  to  buoyancy  effects;  wake  defect  created  downstream  of 
the  bubble  wire;  and  resolution  of  the  digitized  image.  Furthermore,  the 
uncertainty  in  the  streamwise  velocity,  arising  from  existence  of  a  finite  cross¬ 
stream  velocity  component,  is  actually  less  than  that  corresponding  to  a  single¬ 
element  hot  film  probe  over  certain  regimes  of  operation. 

yj  “ Flow  Visualization  Via  Laser-Induced  Reflection  from  Bubble  Sheets ”  by  C. 
Magness,  T.  Utsch,  and  D.  Rockwell,  AIAA  Journal.  Vol.  28,  No.  7,  1990,  pp. 
1199-1200. 

The  structure  of  steady  and  unsteady  aerodynamic  flows  can  be 
characterized  by  generating  fluid  markers  at  desired  locations  in  the  flow  and 
tracking  them  in  three-dimensional  space.  The  hydrogen  bubble  technique  allows 
localized  injection  of  both  continuous  and  interruptive  fluid  markers.  The 
technique  described  herein  employs  multiple  sheets  of  hydrogen  bubbles  in 
conjunction  with  laser  sheet  illumination.  Arbitrary  cross-sections  of  the  unsteady 
flow  past  an  oscillating  delta  wing  can  be  characterized.  Using  a  phase-referencing 
technique,  it  is  possible  to  relate  the  visualization  of  various  cross-sections  along 
the  wing  at  a  given  value  of  instantaneous  angle  of  attack. 

yj  “Leading-Edge  Vortices  on  a  Pitching  Delta  Wing”  by  R.  Atta  and  D.  Rockwell, 
AIAAlSMMl,  Vol  28,  No.  6,  1990,  pp.  995-1004. 

Flow  past  a  stationary  delta  wing  of  moderate  sweep  angle  gives  rise  to 
leading-edge  vortices  having  non-circular  (quasi-elliptical)  cross-section,  with  the 
inclination  of  the  vortex  cross-section  changing  the  streamwise  distance  along  the 
wing.  Representations  of  vortex  swirl  angle  are  compared  with  angle  of  attack  and 
change  of  position  in  vortex  breakdown.  Vortex  breakdown  consistently  occurs  in 
an  organized  spiral  modp;  the  corresponding  maximum  values  of  swirl  angle 
approximate  those  of  vortex  breakdown  in  axisymmetric  internal  flows. 

When  the  wing  undergoes  pitching  motion,  there  occur  two  basic  types  of 
vortex  development.  At  low  frequencies,  the  vortex  core  develops  in  the  upstream 
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direction  towards  the  apex;  at  high  frequencies,  there  is  ejection  of  the  leading-edge 
of  the  vortex  core  from  the  apex  in  the  downstream  direction.  At  sufficiently  high 
frequencies,  the  classes  of  vortex  breakdown  are  distinctly  different  from  those 
occurring  in  steady  flow. 

The  evolution  in  time  of  the  leading-edge  vortices  shows  pronounced 
dynamic  hysteresis:  the  flow  structure  of  the  vortex  is  different  during  the  up-  and 
downstroke  motion  of  the  wing,  even  for  very  low  values  of  reduced  frequency.  The 
direction  of  the  dynamic  hysteresis  loop  (clockwise  or  counterclockwise)  can  be 
related  directly  to  the  type  of  vortex  development. 

During  unsteady  development  of  the  vortex  core,  the  axial  velocity  on  the 
centerline  of  the  vortex  is  constant  with  distance  from  the  apex  of  the  wing.  As 
the  location  of  vortex  breakdown  is  approached,  there  is  abrupt  deceleration  over  a 
distance  of  the  order  of  ten  percent  of  the  chord  of  the  wing.  The  distributions  of 
velocity  in  this  region  suggest  a  universal  form  for  a  wide  range  of  angle  of  attack 
in  reduced  frequency. 

\J  “Transient  Response  of  Leading-Edge  Vortices  to  Localized  Suction ”  by  K. 
Parmenter  and  D.  Rockwell,  AIAA  Journal  Vol.  28,  No.  6,  1990,  pp.  1181-1182. 

This  investigation  addresses  the  restabilization  of  a  leading-edge  vortex  on  a 
delta  wing  by  transient  suction.  Suction  is  applied  through  a  probe  well 
downstream  of  the  onset  of  vortex  breakdown.  Of  particular  interest  is  the  overall 
time  required  for  stabilization  of  the  vortex  and  hysteresis  of  the  vortex  response. 
The  focus  of  this  investigation  is  the  response  of  the  vortex  in  relation  to  the 
abrupt  onset  and  cessation  of  suction.  The  following  aspects  are  addressed: 
visualization  of  the  unsteady  evolution  of  the  vortex  core;  the  nature  of  the 
hysteresis  of  the  core  development;  and  the  overall  response  time  from  the  onset  of 
suction  to  the  stabilization  of  the  core  as  a  function  of  suction  amplitude  and 
suction  probe  location. 

yj  “Flow  Structure  Generated  by  Oscillating  Delta  Wing  Segments ”  by  T.  Utsch  and 
D.  Rockwell,  AIAA  Journal  of  Aircraft.  Vol.  27,  No.  6,  1990,  pp.  574-576. 

Oscillation  of  two  delta  wing  segments  about  a  common  axis  by  a  concentric 
drive  train  arrangement  allows  forcing  of  each  wing  segment  in  an  arbitrary 
manner.  The  resultant  flow  structure  is  visualized  using  a  new  technique  of  laser- 
induced  reflection  from  small  hydrogen  bubbles  generated  upstream  of  the  moving 
wing  segments.  The  scanning  laser  sheet,  which  is  translated  in  the  direction  of 
flow,  allows  instantaneous  slices  of  the  flow  structure  to  be  captured.  This 
approach  leads  to  phase-locked,  instantaneous  velocities,  allowing  evaluation  of  the 
circulation  of  the  leading-edge  vortices. 

For  the  case  of  sinusoidal  motion  of  the  wing  segments,  it  is  possible  to 
generate,  at  zero  angle  of  attack,  substantial  values  of  circulation.  Oscillation  of 
the  wing  segments  at  finite  angle  of  attack  produces  drastic  hysteresis  of  the  vortex 
formation  during  the  upstroke  and  downstroke  motions  of  the  wing. 
Correspondingly,  the  circulation  during  the  upstroke  motion  is  much  higher  than 
that  during  the  downstroke. 
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PUBLICATIONS  SUBMITTED 

yj  u Control  of  Vortex  Structure  on  a  Delta  Wing  by  Small  Amplitude  Perturbations  of 
Angle-of- Attack,”  by  C.-H.  Kuo,  C.  Magness,  and  D.  Rockwell,  submitted  to 
Journal  of  Fluid  Mechanics 

Small  displacement  amplitude  of  the  leading  (one  degree)  sinusoidal 
variations  in  angle  of  attack  can  substantially  alter  the  structure  of  the  leading- 
edge  vortices  on  a  highly-swept  delta  wing.  These  perturbations  are  applied  over  a 
frequency  range  which  includes  the  inherent  instability  frequencies  of:  the  shear 
layer  separating  from  the  leading-edge;  and  the  vortex  breakdown  process.  Both 
upstream  and  downstream  movement  of  the  location  of  vortex  breakdown  can  be 
affected.  Correspondingly,  the  ratio  of  the  characteristic  swirl  and  axial  velocities 
can  change  by  as  much  as  a  factor  of  two  from  that  value  for  the  unperturbed 
wing. 

Detailed  characterization  of  the  spectral  content  of  the  unsteady  vortex 
upstream  and  downstream  of  location  of  vortex  breakdown  show  highly  organized 
spectral  components  at  the  controlled  excitation  frequency;  moreover,  as  many  as 
ten  higher  order  harmonics  can  occur.  This  highly  coherent,  phase-locked  structure 
persists  along  the  entire  axial  extent  of  the  vortex  breakdown  region. 

The  foregoing  observations  of  the  instantaneous  and  mean  flow  are 
interpreted  with  the  aid  of  flow  visualization  of  the  instantaneous  flow  structure. 

PUBLICATIONS  IN  PREPARATION 

s/  “ Quantitative  Visualization  of  Unsteady  and  Separated  Aerodynamic  Flows ”  by  D. 
Rockwell,  R.  Atta,  R.  Curtis,  C.  Hefele,  A.  Klink,  R.  Lawson,  C.  Magness  and  T. 
Utsch  (submitted  to  Experiments  in  Fluids) 

Tracking  of  hydrogen  bubble  markers  in  steady  and  unsteady  three- 
dimensional  flows  can  lead  to  quantitative  information  on  the  flow  structure  in 
terms  of  velocity  fields,  streamlines,  streamsurfaces,  and  timeline  surfaces.  This 
flow  structure  is  obtained  by  dual-image  and  single-image  techniques,  employing 
diffuse  stroboscopic  lighting  and  a  scanning  laser  sheet  respectively.  Vortex  flows 
on  stationary  and  oscillating  delta  wing  segments  are  addressed  using  these 
approaches. 

These  investigations  make  use  of  an  integrated  active  control-quantitative 
flow  visualization  system  with:  high  resolution  stepping  motors  that  force  the 
wing/wing  segments;  a  custom-designed  bubble  generator  for  producing  hydrogen 
bubble  markers  in  terms  of  pulsed  grids  or  arrays  of  bubbles;  and  a  movable 
optics  system  that  translates  the  scanning  laser  sheet.  All  three  of  these  elements 
of  the  integrated  system  are  controlled  by  a  central  computer  in  the  laboratory. 

Construction  of  three-dimensional  flow  images  is  done  with  image 
processing  techniques  and  computer-aided  methods  applied  to  slices  of  the  three- 
dimensional  flow  field.  In  turn,  these  images  can  be  displayed  in  recil  time  from 
arbitrary  perspectives  on  systems  in  the  CAD  laboratory. 
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s/  “Response  of  Vortex  Breakdown  on  a  Pitching  Delta  Wing  Subjected  to  Arbitrary 
Forms  of  Pitching  Motion  ”  by  C.  Magness  and  D.  Rockwell 

This  manuscript  will  address  the  principal  features  of  the  vortex 
breakdown  response  on  a  pitching  delta  wing.  In  particular,  the  form  of  the 
pitching  motion,  involving  simple,  pitch- up  or  pitch-down  maneuvers,  in 
comparison  with  continuous  (successive)  pitch-up  and  pitch-down  maneuvers  will 
be  addressed.  To  complement  these  overall  observations  of  the  response  of  the 
vortex  breakdown,  velocity  fields  and  separation,  streamlines  of  the  cross-sectional 
structure  of  the  vortex  will  be  presented.  The  time  lags  of  the  vortex  structure 
for  a  pitching  wing,  relative  to  that  of  a  corresponding  stationary  wing,  will  be 
addressed. 

y  “Unsteady  Topology  in  the  Leading-Edge  Vortex”  by  C.  Magness  and  D.  Rockwell 

Using  critical  point  theory,  the  instantaneous  topology  is  defined  for  a  delta 
wing  undergoing  pitching  motion.  This  approach  allows  identification  of  the 
saddle  points  and  foci  of  the  instantaneous  streamline  pattern  at  given  cross- 
sections  of  the  flow  at  a  defined  angle  of  attack.  Although  analysis  is  still 
underway,  it  appears  that  the  unsteady  topology  is  substantially  different  from 
what  one  would  expect  on  the  basis  of  quasi-steady  interpretation. 

sj  “ Unsteady  Crossflow  on  a  Delta  Wing  Using  Particle  Image  Velocimetry”  by  C. 
Magness  and  D.  Rockwell 

A  technique  of  laser  scanning  is  employed  to  illuminate  micron-sized 
particles;  photographs  of  these  particles  provide  the  instantaneous  velocity 
distribution  over  the  entire  plane  of  the  flow.  Using  this  approach,  it  is  possible 
to  characterize  the  instantaneous  flow  structure  over  the  entire  crossflow  plane. 
Although  data  reduction  is  still  underway,  it  appears  as  though  it  will  be  possible 
to  characterize  the  well  known  hysteresis  effect  during  pitch-up  and  pitch-down 
motion  in  terms  of  degree  of  coherence  (unstalled  vs.  stalled)  of  the  leading-edge 
vortex. 


NON-JOURNAL  ARTICLES 

yj  “Control  of  Leading- Edge  Vortices  on  a  Delta  Wing”  by  C.  Magness,  0. 
Robinson,  and  D.  Rockwell,  AIAA  Paper  89-0999,  presented  at  AIAA  Shear  Flow 
Control  Conference,  March  13-16,  Tempe,  Arizona. 

The  response  of  vortex  breakdown  on  a  pitching  delta  wing  is  examined  for 
various  classes  of  ramp  motion,  including  pitch-up,  pitch-down,  continuous  pitch- 
up  and  pitch-down,  and  combinations  of  ramp  pitching  rates.  Depending  upon 
which  of  these  types  of  controlled  motion  is  imposed,  there  can  occur  varying 
degrees  of  phase  shift  between  the  onset  and  development  of  vortex  breakdown 
and  the  instantaneous  angle  of  attack  of  the  wing.  These  observations  suggest  that 
by  properly  tailoring  the  functional  form  of  the  pitching  maneuver,  the  puase  sniff 
of  the  vortex  breakdown  can  be  exploited  to  optimize  the  loading  on  Lhe  wing. 
Analogous  forcing  and  response  of  vortex  breakdown  on  a  stationary  wing  to 
localize  suction  and  blowing  is  also  addressed. 

Detailed  insight  into  the  flow  structure  is  acquired  by  direct  particle 
tracking  in  conjunction  with  computer-aided  construction  techniques.  This 
approach  allows  rapid  approximations  to  the  three-dimensioneil  nature  of  the 
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leading-edge  vortex  from  planar  (later  sheet)  slices  through  the  flow.  The  phase 
shift  between  the  instantaneous  structure  of  the  leading-edge  vortex  and  the 
instantaneous  position  of  the  wing  can  also  be  characterized  using  this  approach. 

sj  “On  Unsteady  Flow  Structure  from  Swept  Edges  Subjected  to  Controlled  Motion ”  by 
D.  Rockwell,  R.  Atta,  C.-H.  Kuo,  C.  Hefele,  C.  Magness  and  T.  Utsch. 
Proceedings  of  Workshop  II  on  Unsteady  Separated  Flow  held  at  U.  S.  Air  Force 
Academy,  Colorado  Springs,  Colorado,  July  26-80,  1987.  Appearing  as  Frank  J. 
Seiler  Research  Laboratory  Report  FJSRL-TR-80-0004,  September,  1988,  Project 
2307-F1-38  Air  Force  Systems  Command,  United  States  Air  Force,  pp.  299-312. 

Recent  developments  in  the  direction  of  new  experimental  methods  and 
insight  into  the  mean  and  unsteady  vortex  structure  are  provided  in  synoptic 
form. 

Development  of  an  integrated  active  control  and  quantitative  measurement 
system  has  aimed  towards  relating  the  instantaneous  flow  structure  to  the 
instantaneous  wing  motion.  Active  control  of  the  flow  structure  of  leading-edge 
vortices  on  a  delta  wing  and  delta  wing  segments  through  both  small  and  large 
(displacement)  amplitude  pitching  motion  is  attainable  if  the  excitation  frequency 
is  properly  selected.  Such  excitation  should  consider  the  following  characteristic 
frequencies:  the  fundamental  frequency  (or  its  subharmonics)  of  the  unstable 
shear  layer  separating  from  the  edge  of  the  wing;  the  frequency  of  vortex 
breakdown;  and  the  characteristic  frequency  of  the  separation  zone  downstream  of 
the  onset  of  vortex  breakdown. 

“ Flow  Visualization  and  Its  Interpretation”  by  D.  Rockwell,  R.  Atta,  L.  Kramer, 
R.  Lawson,  D.  Lusseyran,  C.  Magness,  D.  Sohn,  and  T.  Staubli.  Proceedings  of 
the  AGARD  Symposium  on  Aerodynamic  and  Related  Hydrodynamic  Studies 
Using  Water  Facilities,  Monterey,  California,  1986,  AGARD  CP-413,  pp.  29-1  - 
29-13. 

Unsteady  two-  and  three-dimensional  flow  structure  at  leading-  and 
trailing-edges  of  bodies  can  be  characterized  effectively  using  recently  developed 
techniques  for  acquisition  and  interpretation  of  flow  visualization.  The  techniques 
addressed  herein  include:  flow  image-surface  pressure  correlation;  three- 

dimensional  reconstruction  of  flow  structure  from  flow  images;  and  interactive 
interpretation  of  flow  images  with  theoretical  simulations.  These  techniques  can 
be  employed  in  conjunction  with:  visual  correlation  and  ensemble  averaging,  both 
in  a  given  image  and  between  images;  recognition  of  patterns  of  flow  structure 
from  images;  and  estimates  of  velocity  eigenfunctions  from  images. 

,J  “ Control  of  the  Leading-Edge  Vortices  on  a  Delta  Wing”  by  C.  Magness,  0. 
Robinson,  and  D.  Rockwell,  Proceedings  of  the  NASA/AFOSR/ARO  Workshop 
on  Physics  of  Forced  Separation,  NASA  Ames  Research  Center,  Moffett  Field, 
California,  April  17-19,  1990  (ed.  L.  W.  Carr) 

The  unsteady  flow  structure  of  leading-edge  vortices  on  a  delta  wing  has 
been  investigated  using  new  types  of  experimental  techniques,  in  order  to  provide 
insight  into  the  consequences  of  various  forms  of  active  control.  These 
investigations  involve  global  control  of  the  entire  wing  and  local  control  applied  at 
crucial  locations  on  or  adjacent  to  the  wing.  Transient  control  having  long  and 
short  time-scales,  relative  to  the  convective  time-scale  C/Uqq,  allows  substantial 
modification  of  the  unsteady  and  time-mean  flow  structure. 
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Global  control  at  long  time-scale  involves  pitching  the  wing  at  rates  an 
order  of  magnitude  lower  than  the  convective  time-scale  C/Ugo,  by  large 
amplitudes.  The  functional  form  of  the  pitching  maneuver  exerts  a  predominant 
influence  on  the  trajectory  of  the  feeding  sheet,  the  instantaneous  vorticity 
distribution  and  the  instantaneous  location  of  vortex  breakdown. 

Global  control  at  short  time-scales  of  the  order  of  the  inherent  frequency  of 
the  shear  layer  separating  the  leading-edge  and  the  natural  frequency  of  vortex 
breakdown  shows  that  “resident”  response  of  the  excited  shear  layer-vortex 
breakdown  system  is  attainable.  The  spectral  content  of  the  induced  disturbance 
is  preserved  not  only  across  the  entire  core  of  the  vortex,  but  also  along  the  axis  of 
the  vortex  into  the  region  of  vortex  breakdown.  This  unsteady  modification 
results  in  time-mean  alteration  of  the  axial  and  swirl  velocity  fields  and  the 
location  of  vortex  breakdown. 

Localized  control  at  long  and  short  time-scales  involves  application  of 
various  transient  forms  of  suction  and  blowing  using  small  probes  upstream  and 
downstream  of  the  location  of  vortex  breakdown,  as  well  as  distributed  suction 
and  blowing  along  the  leading  edge  of  the  wing  applied  in  a  direction  tangential  to 
the  feeding  sheet.  These  local  control  techniques  can  result  in  substantial 
alteration  of  the  location  of  vortex  breakdown;  in  some  cases,  it  is  possible  to 
accomplish  this  without  net  mass  addition  to  the  flow  of  the  field. 
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Abstract 

The  response  of  vortex  breakdown  on  a  pitching 
delta  wing  is  examined  for  various  classes  of  ramp 
motion  including:  pitch-up,  pitch-down,  continuous 
pitch-up  and  -down;  and  combinations  of  ramp  pitching 
rates.  Depending  upon  which  of  these  types  of 
controlled  motion  is  imposed,  there  can  occur  varying 
degrees  of  phase  shift  between  the  onset  and 
development  of  vortex  breakdown  and  the 
instantaneous  angle  of  attack  of  the  wing.  These 
observations  suggest  that  by  properly  tailoring  the 
functional  form  of  the  pitching  maneuver,  the  phase 
shift  of  the  vortex  breakdown  can  be  exploited  to 
optimize  the  loading  on  the  wing.  The  analogous 
forcing  and  response  of  vortex  breakdown  on  a 
stationary  wing  to  localized  suction  and  blowing  is  also 
addressed. 

Detailed  insight  into  the  flow  structure  is 
acquired  by  direct  particle  tracking  in  conjunction  with 
computer-aided  construction  techniques.  This  approach 
allows  rapid  approximations  to  the  three-dimensional 
nature  of  the  leading-edge  vortex  from  planar  (laser 
sheet)  slices  through  the  flow.  The  phase  shift  between 
the  instantaneous  structure  of  the  leading-edge  vortex 
and  the  instantaneous  position  of  the  wing  can  also  be 
characterized  using  this  approach. 

1.  Introduction 

Of  critical  importance  in  optimizing  the 
maneuverability  of  delta  wing  aircraft  is  knowledge  of 
the  unsteady  structure  of  the  leading-edge  vortices, 
including  the  occurrence  of  vortex  breakdown.  It  is  well 
known  that  this  vortex  breakdown  process  involves  an 
abrupt  transition  from  a  jet-like  to  a  wake-like  core, 
accompanied  by  a  marked  increase  in  the  level  of 
turbulence.  Vortex  breakdown  on  stationary  wings  has 


been  investigated  extensively,  and  is  assessed  by 
Wedemeyer1.  Woffelt2  demonstrated  that  the  phase 
lag,  or  hysteresis,  of  the  vortex  breakdown  response 
relative  to  the  wing  motion  could  produce  hysteresis  of 
the  forces  acting  upon  the  wing.  For  delta  wings 
subjected  to  sinusoidal  motion  at  high  reduced 
frequency,  Atta  and  Rockwell3  found  that  the 
maximum  length  of  the  vortex  core  prior  to  breakdown 
could  occur  near  the  maximum  angle  of  attack,  as 
opposed  to  the  minimum  angle  of  attack  expected  from 
quasi-steady  considerations.  Moreover,  at  sufficiently 
high  reduced  frequency,  the  vortex  core  exists  over  only 
a  portion  of  the  oscillation  cycle.  The  detailed  nature  of 
the  vortex  development  as  a  function  of  reduced 
frequency  is  addressed  by  Atta  and  Rockwell4.  At  low 
reduced  frequency,  the  vortex  core  develops  in  the 
upstream  direction  towards  the  apex,  whereas  at 
sufficiently  high  reduced  frequency  the  leading-edge  of 
the  vortex  core  is  ejected  in  the  downstream  direction. 
There  is  substantial  phase  shift  of  the  evolution  of  the 
vortex  core  with  respect  to  the  instantaneous  position 
of  the  wing,  extending  to  very  low  reduced  frequencies 
K  =  TfC/Uoo  =  0.025.  Delta  wings  subjected  to  ramp 
motion  in  the  pitching  mode  were  investigated  by 
Reynolds  and  Abtahi5.  They  found  that,  at  sufficiently 
high  angles  of  attack,  the  response  of  the  vortex 
breakdown  has  a  very  long  time  scale.  It  is  an  order  of 
magnitude  larger  than  the  convective  time  scale  at  low 
angles  of  attack,  where  breakdown  occurs  downstream 
of  the  trailing-edge.  This  long  time  scale  is  associated 
with  a  propagation  speed  of  the  position  of  vortex 
breakdown  that  is  much  smaller  than  the  convective 
speed  of  the  vortex  core. 

The  location  of  vortex  breakdown  on  a 
stationary  wing  can  be  subjected  to  analogous  types  of 
control  in  the  form  of  suction  (Parmenter  and 
Rockwell6)  or  blowing  (Visser  et  al7)  applied  locally 
along  the  axis  of  the  vortex  or  along  the  leading-edge. 
Of  particular  interest  in  this  investigation  are  responses 
of  the  vortex  breakdown  to  abrupt,  ramp-type  suction 
and  blowing. 
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An  objective  of  this  investigation  is  to:  examine 
the  overall  response  of  the  vortex  breakdown  to  several 
classes  of  controlled  ramp  (linear)-type  motions  in  the 
pitching  mode,  including  pitch-up  and  pitch-down 
motions,  continuous  pitch-up  and  pitch-down 
maneuvers,  and  abrupt  changes  in  the  ramp  rate  of  the 
pitching  motion;  and  to  study  the  response  to 
analogous,  transient  suction  and  blowing  of  the  vortex 
on  a  stationary  wing.  These  observations  will  be 
complemented  by  detailed  characterization  of  the 
structure  of  the  vortex  cross-section,  relative  to  that  on 
the  corresponding  stationary  wing.  Selected  features  of 
the  instantaneous  structure  of  the  vortex  cross-section 
during  vortex  breakdown  will  also  be  addressed. 

2.  Experimental  Systems  and  Techniques 

Experiments  were  carried  out  in  a  water  channel 
having  a  cross-section  610  mm  by  914  mm  and  a 
maximum  freestream  velocity  of  335  mm/s.  As 
represented  in  the  schematic  of  Figure  1,  delta  wings 
having  a  sweep  angle  A  =  75*  and  chord  lengths  C  = 
241  and  508  mm  and  thickness  t  =  12  mm  were 
employed.  Each  delta  wing  was  flat  on  the  upper 
(downstream)  side  and  beveled  at  an  angle  of  40*  on  the 
lower  (upstream)  side.  The  wing  was  mounted  via  a 
9.5  mm  diameter  sting  to  a  custom-designed  pitching 
mechanism  that  allows  arbitrary  pitching  axis  and 
functional  form  of  the  pitching  maneuver.  For  the 
results  given  herein,  the  wing  of  C  =  241  mm  was 
employed  at  a  freestream  velocity  of  51  mm/s, 
corresponding  to  a  Reynolds  number  Re  =  UooC/t'  = 

I. 2  x  104.  Experiments  conducted  over  the  Reynolds 

number  range  1.2  x  104  <  Re  <  3.6  x  104  showed 

undetectable  effects  of  Reynolds  number,  supporting  the 
concept  that  the  onset  and  development  of  vortex 
breakdown  is  predominantly  an  inviscid  process,  in 
agreement  with  the  survey  of  Atta8  over  a  much  wider 
range  of  Reynolds  number.  , 

The  wing  was  pitched  about  mid-chord  at 
dimensionless  rates  A  =  aC/2Uoo  having  values  0.0125 

<  A  <  0.15.  Angles  of  attack  a  over  the  range  5*  <  a 

<  55*  were  employed,  the  upper  limit  corresponding  to 
the  occurrence  of  vortex  breakdown  at  the  apex  of  the 
stationary  wing. 

Flow  visualization  involved  interactive  use  of 
three  techniques:  dye  visualization;  hydrogen  bubbles; 
and  solid  (Pliolite)  particles.  For  dye  injection,  the  dye 
marker  was  allowed  to  seep  from  0.8  mm  diameter  holes 
located  on  the  upper  surface  of  the  wing  at  a  distance 

II. 6  mm  along  the  leading-edge  from  the  apex,  and  1.5 
mm  inboard  from  the  leading-edge.  This  location 
insured  that  there  was  a  continuous  supply  of  marker  to 
the  central  portion  of  the  developing  vortex.  In 
employing  the  hydrogen  bubble  technique,  a  periodically 
crimped  wire  of  25  micron  diameter  was  located  along 
the  underside  of  the  leading-edge,  so  as  to  allow 
continuous  generation  of  streamlines/streaklines  at  a 
number  of  locations  along  the  edge.  Lighting  for  both 
the  dye  and  hydrogen  bubble  techniques  involved 
combinations  of  continuous  and  stroboscopic  lighting. 
Images  were  recorded  on  the  high  speed  (Videologic) 
video  system  at  120  frames  per  second,  which  was 
synchronized  with  the  stroboscopic  lighting;  or,  on  a 
Nikon  F-3,  35  mm  camera.  For  the  technique 


employing  Pliolite  particles,  having  typical  diameters  of 
500  microns,  desired  cross-sections  of  the  flow  were 
illuminated  using  a  scanning  laser  sheet.  A  two  watt 
Argon-ion  laser  operating  in  the  multi-line  mode 
impinged  upon  a  mirror  oscillating  at  frequencies  in  the 
range  of  30  to  40  c/s;  the  exact  value  of  the  scan  rate 
depended  upon  the  desired  particle  displacement  on  the 
image. 

An  overview  of  the  laser  scanning  system  is 
given  in  Figure  2.  A  system  of  mirrors  was  devised 
such  that  the  scanning  frequency  and  the  orientation  of 
the  scanning  sheet,  involving  translation  in  two 
orthogonal  directions  as  well  as  rotation,  could  be 
accomplished  without  refocussing  the  optics  or 
repositioning  any  of  the  mirrors.  The  image  was 
reflected,  by  another  system,  to  a  camera  located 
exterior  to  the  water  channel.  The  mirror  system 
generating  the  laser  scanning  sheet,  the  instantaneous 
angle  of  attack  of  the  wing,  and  the  firing  of  the  camera 
were  controlled  by  an  external  frequency  generator  and 
a  microcomputer  system. 

For  the  studies  described  herein,  most  of  the 
laser  scans  were  taken  at  the  midplane,  coincident  with 
a  line  parallel  to  and  between  lines  AB  and  DC  in 
Figure  3.  In  order  to  construct  three-dimensional 
images  of  the  instantaneous  flow  structure,  the 
additional  laser  sheet  locations  and  orientations  given  in 
Figure  3  were  employed.  This  concept  is  addressed  in 
the  subsequent  section  on  image  construction. 

3.  Global  Control:  Response  sf  Vortex 

Breakdown  to  Type  of  Pitching  Motion 

Representative  classes  of  ramp  motion  are 
illustrated  in  Figures  4a  and  4b.  In  Figure  4a,  the  case 
of  simple  ramp-up,  ramp-down,  and  ramp-up  with  an 
interrupted  delay  are  illustrated  for  representative 
values  of  dimensionless  pitching  rate  A  =  dC/2Uoo- 
Figure  4b  shows  continuous  pitch-up  and  pitch-down 
motions  for  three  typical  values  of  A  as  well  as  for  the 
hybrid  ramp  function  characterized  by  an  abrupt 
increase  in  pitch  rate  during  the  pitch-up  motion. 

For  the  types  of  maneuvers  described  in  Figure 
4,  both  dye  and  hydrogen  bubble  visualization  were 
carried  out  in  order  to  determine  the  instantaneous 
location  of  vortex  breakdown  x^(t)  in  relation  to 
instantaneous  angle  of  attack  a(t). 

The  essential  features  of  the  response  of  the 
leading-edge  vortex  to  the  prescribed  pitching  motion 
are  indicated  in  Figure  5  showing  plan  views  of  the 
leading-edge  vortex  development  during  a  pitch-up 
maneuver  from  a  —  15*  to  40*  at  a  pitch  rate  A  = 
0.025.  The  total  time  for  the  wing  maneuver  from  a  = 
15’  to  a  =  40’  is  UooAt/C  =  8.66.  The  left  photo 
corresponds  to  an  angle  of  attack  a  =  15*  immediately 
before  the  onset  of  motion;  the  middle  photo 
immediately  upon  completion  of  the  wing  motion  at  a 
—  40*;  and  the  right  photo  corresponding  to  a 
dimensionless  time  Uoot/C  =  0.63  after  the  motion  is 
completed  at  a  =  40*.  In  the  left  photo  there  is  no 
indication  of  vortex  breakdown;  it  occurs  downstream 
of  the  trailing-edge.  In  the  middle  photo,  the  location 
of  breakdown  has  moved  just  upstream  of  the  trailing- 
edge,  and  in  the  right  photo,  it  has  reached  its 
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equilibrium  position  well  upstream  of  the  trailing-edge. 
This  time  shift,  or  phase  lag,  is  central  to  the 
characteristics  of  the  vortex  response  described  in  the 
following. 

Overviews  of  the  vortex  breakdown  location 
x. /C  as  a  function  of  dimensionless  time  tUoo/C  are 
illustrated  in  Figure  6  for  simple  pitch-up,  pitch-down 
and  continuous  pitch-up-down  motions  of  the  wing. 
For  the  simple  pitch-up  and  -down  maneuvers,  the 
location  of  vortex  breakdown  is  allowed  to  relax  to  its 
equilibrium  position.  On  the  other  hand,  for  continuous 
pitch-up-down  motion,  the  downward  motion 
commences  immediately  after  the  upward  motion  ceases 
and  there  is  no  opportunity  for  relaxation  of  the  vortex 
breakdown.  The  extreme  values  of  reduced  pitch  rate 
are  considered:  aC/2Uoo  =  0.025  and  0.15.  For  the 
simple  pitch-up  motion,  the  slope  d(xb/C)/d(tUoo/C) 
is  large  during  the  initial  part  of  the  motion,  then 
substantially  decreases.  For  the  pitch-down  motion,  the 
converse  occurs.  Similar  trends  are  evident  from 
inspection  of  the  plots  of  Reynolds  and  Abtahi®  for 
simple  pitch-up  and  pitch-down  motion.  Irrespective  of 
whether  one  considers  pitch-up  or  -down  motion,  the 
axial  velocity  of  the  location  of  vortex  breakdown 
( corresponding  to  the  aforementioned  slope)  has  its 
lowest  value  when  the  location  of  vortex  breakdown  is 
furthest  upstream  coinciding  with  the  existence  of  large- 
scale  breakdown  over  a  substantial  portion  of  the  wing 
surface.  Higher  velocities  tend  to  occur  when  the 
breakdown  is  on  the  downstream  portion  of  the  wing. 
When  the  continuous  pitch-up-down  motion  of  the  wing 
is  invoked,  the  breakdown  location  remains  downstream 
of  the  midchord  and  low  values  of  the  axial  velocity  uj 
the  breakdown  location  are  avoided.  (Exceptions  occur 
when  the  breakdown  location  is  near  the  trailing-edge.) 
A  time  scale  T  may  be  defined  as  the  total  time  elapsed 
from  upstream  movement  of  the  vortex  breakdown 
location  across  the  trailing-edge  (xb/C  =  1)  during  the 
pitch-up  motion  to  its  downstream  movement  across 
the  edge  during  the  pitch-down  motion.  The  value  of 
the  time-scale  T  for  the  continuous  pitch-up-down 
motion  can  be  as  small  as  one-fifth  of  that 
corresponding  to  successive  pitch- up  and  pitch- down 
motions,  evident  for  the  case  aC/2Uoo  =  0.15  in  Figure 
6.  Although  the  plots  of  Figure  6  provide  insight  into 
the  speed  d(xb/C)/d(tU00/C)  of  the  vortex  breakdown 
location,  the  phase  shift  between  the  instantaneous 
breakdown  location  xb/C  and  the  angle  of  attack  a 
must  be  addressed.  For  example,  the  occurrence  of  the 
minimum  d(xb/C)/d(tUoo/C)  =  0  for  the  continuous 
ramp-up-down  motion  is  expected  to  lag  the  minimum 
da/dt  as  0  of  the  maneuver  schedule.  In  the  following, 
this  and  other  related  features  are  addressed. 

Figure  7  shows  the  response  curves  of  vortex 
breakdown  location  xb/C  vs.  angle  of  attack  a  for  the 
cases  of  simple  pitch-up  and  pitch-down  motion.  For 
both  the  upv.  ird  and  downward  motions,  the  vortex 
breakdown  location  is  allowed  to  relax  to  its  static 
value  before  the  wing  is  set  into  motion.  For  all  three 
values  of  reduced  pitching  rate  &C/2Voo  —  0.025,  0.05 
and  0.15,  there  is  substantial  departure  from  the  static 
charae*iristic  for  bot1’  the  pitch-up  and  pitch-down 
motions.  During  the  pitch-up  maneuver,  the 
instantaneous  value  of  xb/C  substantially  overshoots 


the  corresponding  static  values  at  angles  of  attack  a  in 
excess  of  that  value  at  which  the  breakdown  location  x^ 
crosses  the  trailing-edge  of  the  wing,  i.e.  at  a  =  30  . 
Upon  completion  of  the  pitch-up  motion  at  a  =  40', 
relaxation  of  the  locations  of  breakdown  xb/C  at 
constant  angle  of  attack  a  occurs  over  the  range  0.75  < 
xb/C  <  0.30.  For  pitch-down  motion,  there  is  large 
undershoot  of  the  values  xb/C  relative  to  the  static 
characteristic.  Relaxation  after  cessation  of  the  motion 
at  angle  of  attack  a  =  15‘occurs  only  for  the  highest 
pitch  rate;  it  extends  approximately  from  xb/C  =  0.45 
to  a  location  downstream  of  the  trailing-edge.  Taking 
an  overview  of  the  xb/C  vs.  a  characteristics  of  Figure 
7,  it  is  evident  that  the  process  of  relaxation  of  the 
vortex  breakdown  location  to  its  equilibrium  state  at 
constant  ( terminal)  angle  of  attack  occurs  even  for  low 
pitch  rates  during  pitch- up  motion;  this  relaxation  can 
extend  over  a  streamwise  distance  of  the  order  of  one- 
half  the  chord  of  the  wing.  The  same  type  of  relaxation 
occurs  during  pitch- down  motion  only  for  sufficiently 
high  pitch  rate. 

A  comparison  of  the  cases  of  simple  pitch-up, 
pitch-down,  and  continuous  pitch-up  and  -down  motion 
is  illustrated  in  Figures  8a  and  8b.  For  the  pitch-up 
portion  of  the  maneuver,  the  curves  of  xb/C  vs.  a  are, 
as  expected,  very  similar  at  a  given  value  of  pitch  rate. 
During  the  pitch-down  motion,  these  curves  have  a 
similar  shape  only  at  sufficiently  low  values  of  angle  of 
attack.  For  the  initial  portion  of  the  pitch-down 
motion,  there  is  an  initial  decrease  of  the  xb/C  vs.  a 
characteristic  of  the  continuous  pitch-up-pitch-down 
curve,  relative  to  the  simple  pitch-down  curve.  This 
initial  decrease  extends  over  a  larger  range  of  a  at  A  = 
0.15  (from  a  =  40"  to  a  =  29')  than  at  A  =  0.025.  In 
other  words,  the  effect  of  not  allowing  the  vortex 
breakdown  to  relax  to  its  static  state  during  the 
continuous  maneuver  is  to  produce  a  movement  of  the 
breakdown  location  x ^  upstream  towards  the  apex  for 
the  initial  decrease  in  angle  of  attack  a. 

The  effect  of  a  hybrid  ramp  pitch  rate,  relative 
to  the  constant  rate  described  in  the  foregoing,  is  given 
in  Figure  9.  The  hybrid  function  is  defined  in  Figure 
4b.  During  the  motion  of  the  wing,  there  is  an  abrupt 
increase  in  pitch  rate  from  A  =  0.025  to  0.05;  during 
the  pitch-down  motion,  the  rate  is  constant  at  A  = 
0.025.  The  angle  of  attack  varies  between  5'  <  a  < 
55’.  The  hybrid  fraction  was  chosen  such  that  the 
lower  ramp  rate  was  in  effect  for  <*(t)  <  acr,  when  ocr 
is  the  angle  of  attack  of  the  stationary  wing  at  which 
the  vortex  breakdown  occurs  at  the  trailing-edge  (acr 
~  30*).  This  lower  ramp  rate  produces  less  undershoot 
of  the  xb/C  vs.  a  characteristic  relative  to  the  static 
characteristic.  Threafter,  in  the  region  o(t)  >  ocr,  the 
higher  ramp  rate  produces  more  overshoot  of  the  static 
characteristic  than  the  lower  ramp  rate.  During  the 
pitch-up  motion,  it  is  possible  to  obtain  substantially 
higher  values  of  xb/C  relative  to  those  attainable  with 
either  of  the  two  values  of  pitch  rate,  A  as  0.025  and 
0.05,  that  form  the  hybrid.  This  observation  suggests 
that  employment  of  a  hybrid  pitch  rate  accounting  for 
the  critical  angle  of  attack  can  produce  significantly 
higher  values  of  overshoot  of  the  static  characteristic 
vs.  a  than  those  attainable  with  either  of  the 
continuously  applied  pitch  rates  forming  the  hybrid. 


3 


4.  Instantaneous  Structure 


rview  of  technique 


By  employing  the  computer-controlled  wing 
motion  and  camera-trigger  system,  in  conjunction  with 
the  scanning  laser  sheet,  it  was  possible  to  obtain 
instantaneous  slices  of  the  vortical  structure  during  a 
defined  maneuver  of  the  wing.  This  structure  could 
then  be  compared  with  that  on  the  wing  for  successive 
static  values  of  angle  of  attack  a.  The  concept  of 
particle  tracking  involved  multiple,  consecutive  images 
of  each  particle  per  exposure  time  of  the  35  mm  camera, 
accomplished  by  proper  selection  of  the  laser  scanning 
rate.  This  technique  allowed  rejection  of  any 
incomplete  trajectory,  i.e.  one  that  started  or  ended 
outside  of  the  laser  sheet.  The  typical  thickness  of  the 
laser  sheet  was  3  mm,  and  a  nominal  scan  rate  of  40  c/s 
was  employed.  The  exposure  time  of  the  camera  was 
nominally  0.25  seconds. 


Digitized  images  of  the  tracked  particles  were 
obtained  from  enlarged  views  of  the  flow  field  on  203 
mm  x  254  mm  photographs  by  use  of  a  Tektronix  4957 
digitizing  tablet  driven  by  Unigraphics  II.  This  process 
of  data  acquisition  and  digitization  is  currently  being 
refined  in  order  to  provide  a  finer  resolution  and 
complete  automation.  A  technique  using  Young’s  fringe 
analysis  will  be  employed.  Display  of  the  raw  digitized 
data,  interpolated  forms  of  the  data,  and  various  two- 
and  three-dimensional  reconstructions  of  the  images 
were  displayed  on  a  Hewlett  Packard  9000-300  work 
station. 


The  velocity  field  represented  by  the  raw  data  is 
generally  defined  at  random  locations  within  the  plane 
defined  by  the  laser  sheet.  Two  well-known 
interpolation  schemes  for  two-dimensional  scattered 
data  were  employed  to  compute  the  velocity  vectors 
throughout  the  planes  of  interest:  the  Adaptive 

Gaussian  Window  (AGW  algorithm)  and  the  Spline 
Thin  Shell  (STS)  interpolant.  The  STS  scheme 
reproduces  exactly  the  original  values  of  velocity 
corresponding  to  the  locations  of  the  raw  data,  avoids 
spurious  oscillations  of  the  interpolant,  and  is 
independent  of  the  reference  axis.  Moreover,  for  the 
three-dimensional  interpolation  scheme  (discussed 
subsequently),  the  STS  scheme  gives  analytical 
derivatives  for  the  third  velocity  component.  On  the 
other  hand,  a  limitation  of  this  technique  is  that  it  is 
generally  restricted  to  cases  where  the  number  of  data 
points  is  relatively  low,  of  the  order  of  two  hundred. 
The  AGW  method  is  simpler  than  STS  to  implement 
and  can  handle  easily  a  higher  number  of  data  points,  in 
our  case  over  five  hundred  points  per  plane.  The 
derivatives  required  for  the  third  component  in  the 
three-dimensional  reconstruction  can  be  stored  during 
the  interpolation  process  and  the  degree  of  smoothing 
can  be  controlled  through  selection  of  the  proper  value 
of  the  Gaussian  window  width. 

A  representative  velocity  field  obtained  directly 
from  the  particle  tracking  followed  by  interpolation  is 
given  in  the  image  of  Figure  10.  In  this  case,  there  are 
approximately  450  particle  tracks  in  the  original  data, 
obtained  from  overlay  of  four  successive  phase- 
referenced  photos  obtained  during  the  motion  of  the 


wing.  It  should  be  noted  that  phase  referencing  is 
appropriate  only  when  the  flow  structure  is  relatively 
invariant  from  cycle  to  cycle.  Otherwise  a  phase- 
ensemble-)  averaging  process  is  called  for.  In  the  case 
of  the  stationary  wing  at  high  angle  of  attack,  the 
occurrence  of  pronounced  vortex  breakdown  over 
substantial  cross-section  of  the  vortex  core  promoted  a 
high  degree  of  unsteadiness.  In  this  situation,  only 
instantaneous  realizations  of  the  flow  structure,  without 
superposing  or  averaging  data  from  successive  instants, 
are  meaningful.  As  a  consequence,  the  number  of 
particle  tracks  was  of  the  order  of  200. 

4.2  Relaxation  of  flow  structure  during  wing  motion 

It  is  of  interest  to  compare  the  relaxation  of  the 
flow  structure  to  its  steady  state  condition,  as  occurs  on 
a  stationary  wing,  immediately  after  the  cessation  of 
the  motion  of  the  wing.  To  illustrate  this  concept,  the 
representative  case  of  a  simple  pitch-down  motion  from 
a  =  40*  to  a  s  15’  is  addressed.  Figure  11  shows  the 
cross-section  of  the  vortical  flow  structure  at  the 
midchord  immediately  after  cessation  of  the  pitch-down 
motion  at  a  =  15*.  Both  raw  and  interpolated  velocity 
fields  are  illustrated.  In  each  series,  the  upper  image 
corresponds  to  that  at  t*  =  tUgg/C  =  0  after  cessation 
of  the  motion.  Comparing  this  image  with  subsequent 
ones  taken  at  time  intervals  of  At*  =  1.05  shows  that 
the  flow  structure  in  the  middle  and  bottom  images  is 
essentially  the  same.  In  other  words,  in  absence  of  any 
breakdown  within  its  core,  the  detailed  structure  of  the 
vortex  cross-section  relaxes  to  its  steady  state  value 
within  approximately  one  convective  time  scale  C/Ugg 
after  cessation  of  the  wing  motion,  an  observation 
consistent  with  the  qualitative  observation  of 
Lam  bourne  et  al9.  Moreover,  a  distinguishing  feature 
of  the  image  at  t*  =  0,  relative  to  subsequent  images,  is 
the  substantially  higher  elevation  of  the  vortex  center 
above  the  surface  of  the  wing.  Evaluation  of 
corresponding  values  of  circulation  and  vorticity 
concentration  are  currently  underway.  The  rapid 
relaxation  illustrated  in  Figure  11  is  possible  because 
there  is  no  region  of  significant  vortex  breakdown  in  the 
central  portion  of  the  vortex  core.  The  effect  of 
occurrence  of  vortex  breakdown  at  this  cross-section  of 
the  flow  on  the  relaxation  time  is  addressed  in  the 
following. 

The  evolution  of  the  flow  structure  prior  to  the 
cessation  of  the  wing  motion  is  depicted  in  Figure  12. 
Each  pair  of  images  directly  compares  the  flow 
structure  for  the  corresponding  stationary  wing  (top 
image)  with  that  of  the  instantaneous  structure  for  the 
wing  undergoing  the  pitch  maneuver  (bottom  image). 
The  upper  left  set  of  images  at  a  =  35’  is  taken  shortly 
after  the  onset  of  the  ramp-down  motion  of  a  =  40*. 
The  lower  left  set  is  at  a  =  20’  and  the  upper  right  set 
at  a  =  15*,  immediately  upon  cessation  of  the  wing 
motion.  The  striking  difference  between  the  pairs  of 
images  at  a  =  35’  and  20*  is  the  occurrence  of  a  region 
of  vortex  breakdown  over  a  substantial  extent  of  the 
vortex  cross-section  for  the  moving  wing  and  absence  of 
such  breakdown  for  the  stationary  wing.  The  boundary 
of  the  region  of  vortex  breakdown  was  obtained  from 
consideration  of  the  raw  photographs  and  the  raw 
velocity  fields.  The  interface  between  the  randomly 
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disposed  velocity  vectors  of  the  turbulent  breakdown 
region  and  the  organized  velocity  vector  orientation  in 
the  exterior  region  of  the  vortex  core  was  well-defined. 
Within  the  region  of  vortex  breakdown,  the  interpolated 
velocity  field  is  meaningless.  It  is  necessary  to  employ  a 
sufficiently  narrow  window  in  the  Adaptive  Gaussian 
Window  Scheme  in  order  to  obtain  a  sufficiently  close 
match  of  the  interpolated  velocity  field  with  the  raw 
velocity  data  in  the  domain  immediately  exterior  to  the 
vortex  breakdown  region. 

Two  features  of  the  pattern  shown  in  Figure  12 
are  noteworthy.  First  of  all,  for  the  ramp  motion  from 
a  =  35*  to  a  =  20*,  there  has  elapsed  a  total  time  of 
1.95  C/Uoo.  However,  there  is  little  change  in  the 
cross-sectional  structure  of  the  vortex  core.  Except  for 
a  slight  decrease  in  the  extent  of  the  region  of  vortex 
breakdown,  the  structure  of  the  vortex  cross-section  has 
changed  very  little  and  the  distance  of  the  apparent 
centers  of  the  regions  of  vortex  breakdown  from  the 
wing  surface  are  essentially  unaltered.  In  essence, 
existence  of  vortex  breakdown  within  the  core  of  the 
vortex  precludes  significant  relaxation  of  the  detailed 
structure  of  the  vortex  cross-section  within  time  scales 
of  the  order  of  Cf  U0 0.  Ail  of  the  foregoing  patterns  are 
substantially  different  from  their  counterparts  on  the 
stationary  wing.  No  vortex  breakdown  is  evident,  and 
the  centers  of  the  vortices  move  closer  to  the  wing 
surface  as  angle  of  attack  decreases.  The  second  point 
to  be  made  is  that  when  the  region  of  vortex  breakdown 
moves  downstream  of  the  plane  of  observation ,  as 
indicated  in  the  set  of  images  at  the  upper  right  of 
Figure  12,  it  leaves  a  vortical  flow  structure  that  is 
significantly  elevated  above  the  wing  relative  to  its 
counterpart  on  the  stationary  wing.  These  features,  in 
relation  to  the  angle  and  shape  of  the  separation 
streamline  from  the  leading-edge  of  the  wing,  the  overall 
circulation,  and  the  distribution  of  vorticity  are 
currently  being  evaluated. 

4.3  Competition  of  leading-edge  vortices 

Difficulties  in  evaluating  the  instantaneous  flow 
structure  are  particularly  apparent  on  a 
stationary  wing  at  high  angle  of  attack,  or  a  wing 
moving  at  sufficiently  low  pitch  rate.  Figure  13  shows 
the  instantaneous  flow  structure  at  two  successive  times 
on  a  stationary  wing  at  angle  of  attack  40*.  Not  only 
are  there  significant  asymmetries  of  the  regions  of 
vortex  breakdown  occurring  in  the  vortices  from  the  left 
and  right  leading- edges  of  the  wing,  hut  the  cross- 
sectional  shape  and  vertical  distance  of  the  breakdown 
region  above  the  wing  surface  vary  with  time.  For 
example,  in  the  left  image,  the  center  of  the  region  of 
vortex  breakdown  is  approximately  at  the  same  level  as 
on  the  right  side  of  the  wing.  On  the  other  hand,  it  is 
significantly  lower  than  that  on  the  right  side  in  the 
second  image.  Preliminary  considerations  indicate  that 
the  fluctuations  and  instantaneous  value  of  circulation 
are  substantial,  which  would  produce  corresponding 
fluctuations  in  the  rolling  moment  upon  the  wing. 

5.  Local  Control 

Whereas  global  control  involves  forcing  of  the 
leading-edge  vortex  along  its  entire  extent,  local  control 
takes  the  form  of  suction  or  blowing  applied  at  a 


specified  location  within  the  vortex.  In  the  following, 
preliminary  findings  are  addressed  to  transiently  applied 
suction  and  blowing. 

5.1  Suction 

Localized  suction  applied  downstream  of  the 
region  of  vortex  breakdown  can  result  in  restabilization 
of  the  vortex  core,  accompanied  by  downstream 
movement  of  the  onset  of  vortex  breakdown,  as 
indicated  in  Figure  14.  Parmenter  and  Rockwell6 
address  the  effectiveness  of  abruptly  applied  suction, 
accounting  for  the  effects  of  amplitude  of  the  suction 
coefficient,  and  location  of  the  suction  probe  on  the 
time  required  for  restabilization  of  the  vortex  core. 
Figure  15  shows  an  excerpt  from  this  investigation. 
The  angle  of  attack  was  a  =  35*  and  the  Reynolds 
number  Re  =  3.1  x  104.  In  this  case,  the  dimensionless 
suction  coefficient  =  [Vp/Uoo]2(Ap/Aw)  =  0.57,  in 
which  Uqo  =  freestream  velocity,  Vp  =  velocity  of 
suction  probe  inlet,  Ap  =  area  of  probe  inlet,  and  Aw 
s=  surface  area  of  wing.  Restabilization  of  the  vortex 
core  over  a  distance  of  about  a  half  chord  of  the  wing  is 
attainable  at  lower  values  of  suction  coefficient,  down  to 
=  0.035,  the  lower  limit  of  the  pumping  system.  It 
is  anticipated  that  much  smaller  values  of  suction 
coefficient  can  be  employed  when  the  distance  between 
the  suction  probe  inlet  and  onset  of  vortex  breakdown  is 
less  than  the  half  chord  of  the  wing. 

The  flow  visualization  of  Figure  15  shows  the 
structural  response  of  the  vortex  with  time  t*,  where  t* 
=  tUoo/C,  which  ranges  from  t*  =  0.10  (shortly  after 
the  application  of  suction  of  t*  =  0)  to  t*  =  2.76 
marking  complete  stabilization  of  the  core  flow.  This 
stabilization  exhibits  four  distinct  stages:  (a) 
downstream  movement  of  the  onset  of  vortex 
breakdown  and  an  imperceptible  change  of  the  radius  of 
the  breakdown  spiral  (compare  t*  s  0.10  with  t*  = 
1.55);  (b)  continued  downstream  movement  of  the 
spiral  pattern,  accompanied  by  a  decrease  in  the  radius 
of  the  spiral  (compare  t*  =  1.55  with  t*  =  1.95);  (c) 
rapid  decrease  in  radius  of  the  spiral,  eventually 
becoming  indiscernible  and  leaving  a  breakdown  bubble 
upstream  of  a  turbulent  wake  (compare  t*  =  1.95  with 
2.29);  and  (d)  stabilization  of  the  turbulent  breakdown 
region  as  it  is  drawn  into  the  probe  (compare  t*  =  2.29 
with  2.56). 

Further  consideration  of  the  effects  of  localized 
suction  involve  abrupt  cessation  of  the  suction  after  the 
vortex  core  has  been  stabilized.  There  results 
substantial  hysteresis,  relative  to  that  occurring  for 
abrupt  onset  of  suction.  The  magnitude  and  direction 
of  this  hysteresis  depends  on  the  amplitude  of  the 
suction  coefficient  (Parmenter  and  Rockwell6). 

5.2  Blowing 

Application  of  localized  blowing  along  the 
leading-edge  has  been  shown  to  extend  the  streamwise 
length  of  the  vortex  core  prior  to  onset  of  breakdown, 
reviewed  and  demonstrated  by  Visser  et  al7.  The 
purpose  of  this  phase  of  the  investigation  is  to  examine 
the  time-dependent  restabilization  of  the  vortex  core 
after  the  abrupt  onset  of  blowing.  The  concept  of  local 
blowing  along  the  edge  is  illustrated  in  Figure  17.  For  a 


5 


half  delta  wing  at  angle  of  attack  a  =  40*,  the  location 
of  vortex  breakdown  prior  to  the  onset  of  blowing  is  at 
x/C  =  0.40.  A  '.aotr  sheet  cuU  across  the  leading-edge 
vortex  location  x/C  =  0.52.  In  the  case  addressed  here, 
the  amplitude  of  the  blowing  coefficient  is  of  the 
order  of  unity;  much  lower  values  should  be  attainable 
by  optimizing  the  location  and  orientation  of  the 
blowing  or  nozzle.  It  should  be  kept  in  mind  that  the 
ratio  of  the  length  of  the  vortex  core  prior  to 
breakdown  to  the  diameter  of  the  core  for  jet  blowing 
was  of  the  order  of  200.  As  a  consequence,  there  is 
substantial  decay  in  the  centerline  velocity  of  the  jet, 
accompanied  by  widening  of  the  jet.  Consideration 
should  be  given  to  a  blowing  coefficient,  based  on  the 
effective  jet  velocity  and  width  at  the  location  of  vortex 
breakdown. 

The  photos  of  Figure  17  show  raw  particle 
images  of  the  vortex  breakdown  prior  to  the  onset  of 
blowing  and  the  restabilized  cross-section  of  the  vortex 
after  blowing  has  been  employed.  The  left  photo 
corresponds  to  time  t*  =  0  and  the  right  photo  to  t*  = 
1.36.  In  other  words,  the  leading-edge  vortex  is,  at  the 
streamwise  location  shown,  effectively  restabilized  in 
about  one  convective  time  scale  C/U<x>. 

6.  Three-Dimensional  Image  Construction 

6.1  Inverse  technique  using  theoretical  simulation 

A  major  objective  of  this  program  is  to 
construct  the  instantaneous  (or  phase-referenced), 
three-dimensional  velocity  field  from  minimum 
information  on  selected  planes  throughout  the  three- 
dimensional  space  of  interest.  Figure  3  shows  a  plan 
view  of  the  laser  sheet  orientation  that  allows 
reconstruction  of  the  three-dimensional  velocity  field. 
In  principle,  the  acquisition  of  data  on  planes  parallel  to 
and  between  planes  AB  and  CD  in  conjunction  with 
data  on  planes  DG  and  FB  allows  reconstruction  of  the 
three-dimensional  field.  Of  course,  there  are  many 
considerations  with  regard  to  the  uncertainties  arising 
from  data  density  within  each  of  the  planes  and  the 
spacing  and  orientation  of  the  planes  which  are  related 
to  the  local  gradients  of  the  three-dimensional  velocity 
field.  All  of  these  aspects  are  currently  under 
investigation  using  an  inverse  technique,  which  starts 
with  known  theoretical  solutions  of  three-dimensional 
flows  such  as  a  simple  columnar  vortex  and  a  toroidal 
vortex.  Streamlines  corresponding  to  Hill’s10  spherical 
vortex  are  illustrated  in  Figure  18.  This  flow  structure 
has  many  of  the  elements  existing  in  complex  leading- 
edge  flows,  including  a  well-defined  distribution  of 
vorticity.  By  slicing  the  three-dimensional  vortex  at 
selected  locations  and  orientations,  determining  the 
velocity  field  in  each  plane,  then  randomizing  the 
locations  of  the  velocity  vectors,  it  is  possible  to 
simulate  the  velocity  field  deduced  from  randomly 
placed  particle  tracks  that  would  be  generated  in  the 
laboratory  using  the  laser-scanning  technique.  This 
concept  is  illustrated  in  Figure  19,  showing  side,  top, 
and  isometric  views  of  the  streamline  pattern 
superposed  upon  the  randomized  velocity  field  in  a 
plane  oriented  at  45’  with  respect  to  the  axis  of  the 
vortex.  With  a  number  of  planes  of  the  randomized 
velocity  field  at  hand,  it  is  possible  to  reconstruct  the 


three-dimensional  image  and  determine  the  effect  of  the 
local  velocity  gradients  and  laser  plane  orientation  and 
spacing  on  the  accuracy  of  the  reconstruction.  This 
process  is  currently  underway. 

Irrespective  of  whether  laboratory-  or 
theoretically-generated  data  are  employed,  the 
technique  for  determination  of  the  out-of-plane  velocity 
component  is  the  same.  Applying  the  continuity 
equation  at  each  node  of  the  three-dimensional  grid,  the 
derivatives  in  the  x  and  y  directions  of  the  in-plane  (xy 
plane)  velocity  components  u  and  v  are  fitted  by  a  finite 
cosine  series.  The  integrated  form  of  this  continuity 
equation  is  approximated  by  a  modified  finite  sine 
series.  Applying  inverse  discrete  cosine  and  sine 
transforms,  using  a  Fast  Fourier  Transform  technique, 
the  values  of  the  out-of-plane  velocity  component  w  at 
each  node  can  be  computed.  With  the  three- 
dimensional  velocity  field  at  hand,  a  tricubic  spline 
interpolation  can  be  employed  to  obtain  the  three- 
dimensional  streamline  pattern. 

6.2  Construction  af  leading-edge  vortex  from 
experimental  data 

The  three-dimensional  image  construction 
technique,  described  in  the  foregoing,  was  applied  to  the 
leading-edge  vortex  on  the  stationary  wing  at  an  angle 
of  attack  a  =  20*,  for  which  vortex  breakdown  occurs 
downstream  of  the  trailing-edge.  It  should  be  noted 
that,  in  concept,  the  same  technique  illustrated  here  for 
a  steady  flow  field  can  easily  be  applied  to  an  unsteady 
flow  by  using  phase-referencing  or  ensemble-averaging 
techniques. 

Data  were  acquired  in  seven  planes,  oriented 
parallel  to  the  span  of  the  wing  and  passing  vertically 
through  the  wing  surface,  as  depicted  in  Figure  20. 
With  reference  to  Figure  3,  these  planes  can  be  viewed 
as  lying  parallel  to  lines  AB  and  CD.  To  construct  the 
three-dimensional  image,  it  is  necessary  to  provide 
initial  conditions  for  the  out-of-plane  component  by 
acquiring  data  in  a  plane  oriented  at  an  angle  with 
respect  to  those  of  Figure  20.  As  illustrated  in  Figure 
21,  the  plane  cuts  vertically  through  the  wing  at  an 
angle  of  42’  with  respect  to  its  midplane.  An 
approximation  of  the  three-dimensional  streamline 
pattern  is  illustrated  in  Figure  22.  A  solid  model 
representation  of  the  same  leading-edge  vortex  from  a 
different  perspective  is  given  in  Figure  23.  It  should  be 
emphasized  that  the  spacing  between  the  data  planes 
illustrated  in  Figure  20  is  not  small  relative  to  the 
streamwise  distance  over  which  one  complete  revolution 
of  a  typical  three-dimensional  streamline  occurs. 
Consequently,  the  results  displayed  in  Figures  22  and  23 
are  approximate  and  useful  only  for  global 
representations  of  the  flow  structure.  Improvements  in 
the  reconstruction  scheme  and  in  the  density  of  data 
acquisition  are  currently  underway. 

7.  Concluding  fismatka 

Global  control,  whereby  the  motion  of  the  entire 
wing  influences  the  structure  of  the  teading-edge  vortex 
along  its  entire  extent  and  local  control,  involving 
application  of  suction  or  blowing  at  a  defined  location 
within  the  vortex,  have  been  addressed  herein. 
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Emphasis  has  been  given  to  the  unsteady  response  of 
the  vortex  to  transient,  ramp-type  forcing. 

A  common  feature  of  all  of  these  types  of 
control  is  a  time  delay,  or  phase  shift,  between  the 
onset  of  the  control  and  the  eventual  response  of  this 
vortex.  Depending  upon  the  type  of  control,  the 
amplitude  of  the  control  forcing  function,  and  whether 
or  not  vortex  breakdown  occurs,  this  time  delay  may 
range  from  a  value  less  than  C/Ujo  to  a  value  several 
times  C/Uoo'  This  time  delay  has  been  characterized 
with  respect  to:  the  streamwise  development  of  the 
vortex  core,  emphasizing  the  changes  in  location  of  the 
onset  of  vortex  breakdown;  and  the  development  in  the 
cross-stream  plane  of  the  vortical  structure. 

For  global  motion  of  the  entire  wing, 
examination  of  the  instantaneous  location  of  the  vortex 
breakdown,  relative  to  the  instantaneous  position  of  the 
wing,  reveals  that: 

(a)  For  simple  pitch-up  and  pitch-down  motions  of  the 
wing,  the  speed  of  propagation  of  the  vortex 
breakdown  has  its  minimum  value  when  the 
location  of  vortex  breakdown  is  furthest  upstream 
and  it  tends  to  its  largest  value  as  the  location  of 
breakdown  moves  towards  the  trailing-edge. 
Exceptions  occur  near  the  trailing-edge  of  the 
wing,  where  the  discontinuity  of  boundary 
conditions  is  associated  with  low  propagation 
speeds. 

(b)  Variation  of  the  location  of  vortex  breakdown  with 
angle  of  attack  a  shows  large  departures  from  the 
static  characteristic  even  for  low  reduced  pitch 
rates  dC/^Uoo  =  0.025.  The  breakdown  location 
can  relax  to  its  equilibrium  position  over  a  distance 
of  as  much  as  a  half  chord  after  cessation  of  the 
wing  motion. 

(c)  Continuous  pitch-up-down  motions  of  the  wing 
preclude  low  propagation  speeds  of  the  vortex 
breakdown  (except  near  the  trailing-edge)  by  not 
allowing  the  breakdown  to  move  to  upstream 
portions  of  the  wing.  For  the  continuous  ramp-up- 
down  motion  of  the  wing,  the  total  time  scale  for 
the  breakdown  location  to  cross  the  trailing-edge 
during  the  upstroke,  then  again  during  the 
downstroke,  can  be  about  one-fifth  that  for  the 
corresponding  cases  of  simple  pitch-up  and  pitch- 
down  motions  where  full  relaxation  of  the 
breakdown  is  allowed  to  occur  between  the 
motions. 

(d)  For  continuous  pitch-up-down  motions  of  the  wing, 
the  effect  of  not  allowing  the  vortex  breakdown  to 
relax  to  its  static  state  during  the  maneuver  is  to 
produce  upstream  movement  of  the  breakdown 
location  towards  the  apex  for  initial  decreases  in 
angle  of  attack  a. 

(e)  Combinations  of  simple,  ramp-type  motions  to 
form  a  hybrid  pitching  motion  can  produce 
overshoots  of  the  static  Xl  vs  a  characteristic 
beyond  that  attainable  with  any  of  the  simple 
ramp  motions  alone. 


These  features  of  the  streamwise  movement  of 
the  vortex  breakdown  during  global  control  are 
associated  with  variations  in  the  cross-sectional 
structure  of  the  vortex.  The  time  scales  involved  in 
relaxation  of  the  vortex  cross-section  are  strongly 
dependent  upon  existence  of  breakdown  within  the  core. 

(a)  In  absence  of  breakdown,  the  detailed  structure  of 
the  vortex  rapidly  relaxes  to  its  static  state  within 
•ibout  one  convective  time  scale  C/Uoo- 

(b)  Existence  of  breakdown  within  the  core  can 
produce  long  delays  in  the  relaxation;  insignificant 
changes  were  observed  over  as  long  as  two 
convective  time  scales  C/Uqo  while  the  wing  is  in 
motion. 

For  local  control  in  the  form  of  suction  or 
blowing,  substantial  streamwise  movement  of  the  vortex 
breakdown  and  corresponding  restabilization  of  the 
cross-section  of  the  vortex  can  be  achieved  with  a  time 
delay  of  one  to  two  convective  time  scales.  However, 
minimization  of  the  energy  input  to  effect  the  alteration 
of  the  vortex  is  likely  to  result  in  longer  time  scales. 
Optimization  of  these  control  concepts  and 
establishment  of  analogies  with  the  vortex  response  to 
global  control  are  currently  underway. 
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Figure  3:  Location  and  orientation  of  planes  of 
laser  scanning  technique. 


Figure  1:  Schematic  of  delta  wing. 
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Figure  2:  Illustration  of  integrated  flow  visualization-active  control  system 
employing  computer-controlled  wing  motion  and  camera  unit  in  conjunction 
with  laser  scanning  technique. 
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Angle  of  Attack 


Figure  4a:  Representative  pitching  schedules  for 
simple  pitch-up  and  pitch-down  motion  of  wing, 
allowing  for  relaxation  of  leading-edge  vortices 
after  cessation  of  wing  motion. 


Figure  4b:  Representative  pitching  schedules  for 
continuous  pitch-up-down  motion  of  wing. 


Figure  5:  Visualization  of  development  of  leading-edge  vortices  via  the  hydrogen  bubble  technique  for  pitch-up 
motion  from  a  =  15’  to  o  =  40"  at  pitch  rate  crC/2U<x  =  0.025.  Left  photo  corresponds  to  a  =  15’ 
immediately  prior  to  the  onset  of  wing  motion;  middle  photo  taken  at  o  =  40'  immediately  after  cessation  of 
motion;  and  right  photo  represents  time  tUoo/C  =  0.63  after  cessation  of  motion. 
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Figure  6:  Time  history  of  location  of  vortex  breakdown  for  simple  pitch-up,  simple  pitch- 
down,  and  continuous  pitch-up-down  motions  of  wing.  Cases  of  relative  low  (top  graph)  and 
high  (bottom  graph)  pitching  rates  are  illustrated. 


Figure  7:  Location  of  vortex  breakdown 
as  function  of  angle  of  attack  for  three 
representative  values  of  pitching  rate  for 
simple  pitch-up  and  pitch-down  motions, 
where  flow  structure  is  allowed  to  relax  to 
static  condition  of  wing  motion. 
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Figure  8:  Location  of  vortex  breakdown 
as  function  of  angle  of  attack  for  various 
types  of  pitching  motion  at  relatively  low 
(top  graph)  and  high  (bottom  graph) 
pitching  rates. 
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Figure  9:  Location  of  vortex  breakdown 
as  function  of  angle  of  attack  for 
continuous  pitch-up-down  motions  of 
wing.  Hybrid  pitching  schedule  is  defined 
in  Figure  4b. 
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Figure  10:  Representative  velocity  field 
obtained  directly  from  particle  images 
followed  by  interpolation.  Image*  acquired 
at  midchord  and  at  a  =  15*  immediately 
after  cessation  of  pitch-down  motion  from 
a  =  40*.  Pitch  rate  6C/2Uoo  =  0.15. 


Figure  11:  Raw  and  interpolated  data  showing  relaxation  of  instantaneous 
vortical  flow  structure  at  midchord  immediately  after  cessation  of  pitch- 
down  motion  from  a  =  40*  to  a  s  15*  at  pitch  rate  &C/2Uoo  =  0.15. 
Upper  image  corresponds  to  t*  =  tUoo/C  =  0  after  cessation  of  motion; 
middle  image  to  t*  =  1.05;  lower  image  to  t*  =  2.10. 


Figure  12:  Instantaneous  vortical  flow  structure 
at  midchord  at  successive  instants  of  time  during 
pitch  down  maneuver  (bottom  image  in  each  set) 
relative  to  structure  for  stationary  wing  (top 
image  in  each  set).  Upper  left  set  of  images 
corresponds  to  a  =  35*;  lower  left  set  to  a  =  20*; 
and  upper  right  set  to  a  =  15*.  Closed  contours 
represent  boundaries  of  region  of  vortex 
breakdown.  Pitch  rate  AC/2Uoo  =  0-15. 


Figure  13:  Instantaneous  flow  structure 
at  two  successive  times  near  midcbord 
(x/C  =  0.51)  of  stationary  wing  at  a  = 
40*.  Extent  and  location  of  region  of 
vortex  breakdown  oscillates  with  time. 


Figure  14:  Schematic  illustrating  concept 
of  suction  applied  to  vortex  in  region 
downstream  of  onset  of  vortex 
breakdown. 
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Figure  15:  Visualization  of  process  of  stabilization  of  vortex  core  after  abrupt  onset  of  localized  suction  at  t*  = 
0.  Left  boundary  of  photo  is  at  x/C  =  0.55.  Vortex  breakdown  in  absence  of  suction  occurs  at  xj C  =  0.45 
Suction  probe  is  located  at  xp/C  =  0.95.  Amplitude  of  suction  coefficient  Cp  =  0.57.  Re  =  3.1  x  104. 


Figure  16:  Schematic  illustrating  down¬ 
stream  movement  of  vortex  breakdown 
due  to  blowing  along  leading-edge  of  wing 
applied  at  its  apex. 


Figure  17:  Particle  images  across  cross-section  of  leading-edge  vortex  at  laser  sheet  location  x/C  =  0.52  before 
(t*  =  tUoo/C  as  0)  and  after  (t*  =  1.36)  abrupt  onset  of  blowing  applied  along  leading-  edge  of  wing  at  its 
apex.  In  absence  of  blowing,  vortex  breakdown  occurs  at  x^/C  =  0.40.  Re  =  18.8  x  104.  (Photos  by  C.-H. 
Kuo  and  S.  Guzy). 
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Figure  18:  Representation  of  stream¬ 
lines  of  Hill’s  spherical  vortex 
employed  for  simulation  of  three- 
dimensional  construction. 


Figure  19:  Superposition  of  stream¬ 
lines  of  Hill’s  spherical  vortex  and 
randomized  velocity  field  in  plane 
cutting  at  an  angle  of  45*  through  the 
center  of  vortex. 


Figure  20:  Velocity  fields  obtained  directly  from 
particle  images  in  planes  cutting  vertically 
through  wing  which  is  stationary  and  at  angle  of 
attack  a  =  20*. 
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Figure  21:  Velocity  field  in  plane  cutting 

vertically  through  wing  at  angle  of  42*  with 
respect  to  midplane  of  wing. 


Figure  22:  Approximation  of  three-dimensional  Figure  23:  Solids  model  representation  of  leading- 

streamline  pattern  of  leading-edge  vortex  edge  vortex  of  Figure  22. 

obtained  from  information  in  images  of  Figures 
21  and  22. 
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ABSTRACT 

This  investigation  addresses  active  control 
of  the  leading-edge  vortices  formed  from  delta 
wings  and  wing  segments.  A  variety  of 
modifications  of  the  flow  structure  are 
characterized  with  quantitative,  global  flow 
visualization  in  conjunction  with  LDA 
techn i ques . 


INTRODUCTION 


The  structure  of  flow  past  a  stationary  delta 
wing  has  been  studied  extensively  over  the  past 
few  decades.  Early  investigations  established  that 
flow  separation  from  the  leading-edge  of  the  wing 
gives  rise  to  a  spiral-vortex  sheet  on  the  upper 
surface  of  the  wing;  the  consequence  of  this 
vortex  pattern  is  an  increase  in  lift,  above  that 
which  would  occur  without  leading-edge  separation 
(e.g.,  Lee';  Ornberg2;  El  le^)  .  The  study  of 
Kuechemann1* ,  which  involves  a  nonlinear  lifting- 
surface  theory,  describes  this  effect.  More 
recent  studies  (Kuechemann^;  Hoei jmakers°;  Rizzi 
et  a  1 7 ;  Woodson  and  DeJarnette8)  delineate  details 
of  the  primary-secondary  vortex  system  adjacent  to 
the  leading-edge  and  the  corresponding  separation 
and  reattachment  lines. 


The  unsteady  flow  structure  on  a  delta  wing 
subjected  to  controlled  motion  has  received  less 
attention.  The  investigation  of  Lambourne  et  al® 
demonstrated  the  concept  of  phase  shift  between 
development  of  the  vortex  pattern  and  the  stepwise 
motion  of  the  wing  in  the  plunging  mode.  For  the 
case  of  simple  harmonic  motion  in  the  pitching  mode, 
Gad-el -Hak  and  Ho'^  observed  an  analogous  phase 
shift  of  the  growth-decay  cycle  of  the  leading- 
ecge  vortices.  This  phase  shift,  or  hysteresis, 
during  the  course  of  the  wing  oscillation  cycle 
was  characterized  in  terms  of  height  or  thickness 
of  the  dye  blob  marking  the  vortex;  the  degree  of 
hysteresis  was  found  to  be  a  strong  function  of 
reduced  frequency.  Further  aspects  of  the  unsteady 
three-dimensional  structure  on  various  wings  are 
addressed  in  the  subsequent  investigations  of 
Gad-el -Hak  and  Ho">'2.  They  emphasize  the 
importance  of  mutual  induction  between  leading- 
and  trai 1 ing-edge  vortices  in  determining  the 
overall  flow  structure. 


*  Based  on  presentation  at  AFOSR  Workshop  on 
Unsteady  and  Separated  Flows,  U.  S.  Air  Force 
Academy,  Colorado  Springs,  Colorado,  July  26  - 
30,  1987. 


Central  to  the  unsteady  structure  of  the 
leading-edge  vortices  is  the  time-dependent 
development  of  the  vortex  core  along  its  stream- 
wise  extent  and  i fs  eventual  breakdown.  This 
breakdown  is  characterized  by  abrupt  transition 
from  a  jet-like  to  a  wake-1 ikq  core  of  the  vortex, 
accompanied  by  a  substantial  increase  in 
turbulence  activity.  Such  vortex  breakdown  has 
been  investigated  for  flow  past  a  variety  of 
stationary  wing  configurat 
Uedemeyer'3.  For  nonstati 
motion,  Atta  and  Rockwell' 
of  the  maximum  breakdown  position  of  the  core  near 
the  maximum  angle  of  attack,  rather  than  near  the 
minimum  value;  substantial  overshoot  of  the  core 
breakdown  length  relative  to  its  steady  counterpart; 
and  existence  of  the  vortex  core  only  over  a 
fraction  of  the  oscillation  cycle.  There  have  been 
extensive  and  insightful  studies  of  vortex  break¬ 
down  in  various  types  of  ax i symmetric  tube 
arrangements  as,  for  example,  described  by  Garg 
and  Leibovich'5,  Escudier-8  and  Sarpkaya1'.  These, 
as  well  as  other  experimental  studies  and  a  wide 
range  of  theoretical  approaches  beyond  the  scope  of 
this  introduction  are  incisively  assessed  by 
Le  ibovich'8 > '9.  In  essence,  the  investigations  of 
vortex  breakdown  addressed  therein  focussed  on 
steady  inflow  conditions.  Exceptions  include  the 
investigations  of  Lambourne2^  and  Sarpkaya2'  who 
reveal  interesting  aspects  of  transient  vortex 
development  and  overshoot  of  the  vortex  breakdown 
posi tion. 

There  are  a  number  of  unexplored  features  of 
unsteady  vortex  flows  on  oscillating  delta  wing 
and  flap  arrangements.  These  issues  include;  flow 
structure  of  the  vortex  development  and  breakdown; 
evolution  of  the  mean  and  unsteady  velocity  fields 
of  the  vortex  in  the  pre-  and  post-breakdown  regions, 
including  alteration  of  the  velocity  eigenfunctions 
and  spectral  energy  transfer;  and  the  effects  of 
mutual  interference  between  flap/wing  combinations. 
Particularly  helpful  for  investigating  these  aspects 
are  quantitative  flow  visualization  techniques 
which,  when  employed  in  conjunction  with  laser- 
Doppler  and  pressure  measurement  techniques,  can 
provide  both  local  and  global  views  of  the  evolving 
flow  structure. 


ions,  as  reviewed  by 
onary  (sinusoidal)  wing 
*  observed:  occurrence 


EXPERIMENTAL  SYSTEM  AND  TECHNIQUES 

All  experiments  were  carried  out  in  one  of  the 
three  water  channel  systems  custom-designed  for 
unsteady,  separated  flows.  The  cross-sectional 
areas  of  the  test  sections  of  these  channels  ranged 
from  1.5  ft2  to  6  ft2.  In  order  to  allow 
observation  of  the  flow  from  an  arbitrary 


perspective,  all  channels  of  the  test  section  were 
made  entirely  of  glass  or  Plexiglas,  furthermore, 
in  the  case  of  the  largest  channel,  the  entire 
system  upstream  and  downstream  of  the  test  section 
is  made  of  f*V C  and  bronze,  in  order  to  preclude 
problems  with  corrosion.  The  chord  C  length  of 
the  delta  wings  ranged  f-om  C  =  50.8  mn  to  C  = 

508  mm  and  the  corresponding  Reynolds  number 
range  extended  over  5.8  x  10^  <  U„C/v  <  3-9  x  101*. 
Values  of  Re  up  to  Re  =  4.5  x  105  are  attainable 
in  the  large-scale  water  facility,  while  still 
allowing  reasonably  low  absolute  frequencies 
(-.0.5  Hz)  to  be  employed  to  attain  high  values 
of  reduced  frequency.  This  concept  simplifies 
the  mechanical  forcing  of  the  wing,  as  well  as 
data  acquisition.  In  all  cases,  the  viewing 
surface  was  the  flat  machined  surface  of  the  wing; 
the  opposite  side  of  the  wing  was  machined  at  an 
angle  of  a  =  10°  to  20°,  depending  upon  the  wing 
under  consideration.  All  wings  were  of  Plexiglas, 
for  the  technique  of  illumination  by  laser- 
induced  reflection,  the  wing  surface  was  machined 
to  a  one  micron  finish.  In  cases  where  diffuse 
stroboscopic  lighting  was  employed,  the  wing 
was  spray-painted  flat  black  for  employment  of  the 
hydrogen  bubble  technique  and  flat  white  for  use 
of  the  dye  injection  method  of  visualization. 

The  experimental  objectives  require  employment 
of  laser-Ooppier  anemometry  techniques  in 
conjunction  with  quantitative  flow  visualization  in 
order  to  characterize  flow  structure.  Eventual 
measurements  of  the  unsteady  surface  pressure  will 
allow  establishment  of  a  direct  relationship 
between  the  instantaneous  loading  and  the 
instantaneous  flow  structure.  The  first  stage  of 
laser-Ooppier  measurements  described  herein  were 
made  with  a  single  channel  two  watt  Argon-ion 
system  with  a  beam  expander  to  optimize  signal  to 
noise  ratio.  Quantitative  flow  visualization 
involved  tracking  of  hydrogen  bubble  markers.  The 
intention  of  these  tracking  techniques  is  to 
circumvent  the  many  fallacies  associated  with 
interpretation  of  traditional  dye  and  smoke 
injection  techniques  (lusseyran  et  a  1  22;  Hama23). 
Preliminary  pressure  measurements  have  been  carried 
out  with  a  high  sensitivity  Kulite  transducer 
coated  with  a  Paralene  film  to  inhibit  corrosion; 
alternately,  a  PCS  transducer  operating  on  the 
piezo  electric  principle  was  also  employed.  Further 
details  of  these  techniques  are  described  by 
Rockwel  I  et  a  I 

A  primary  objective  of  this  program  is  to 
develop  an  integrated  active  control  system  that 
allows  forcing  of  wings  and  wing  segments  with 
arbitrary  functions  and  simultaneously  accommodates 
acquisition  of  instantaneous  quantitative  flow 
visualization,  laser-Ooppier,  and  pressure 
information.  The  first  version  of  this  system 
has  been  completed.  Forcing  is  provided  by 
Superior  Electric  stepping  motors  with  appropriate 
gear  reduction;  or  with  Compumotors  running  at 
ten  to  twenty  thousand  steps  per  revolution.  These 
motors  are  interfaced  with  a  Zenith-241  computer. 
Also  interfaced  with  the  computer  are:  the  output 
of  the  IDA  system;  the  marker  pulse  generator  and 
switch  system  for  producing  grids  of  visualization 
marker  at  a  desired  phase  of  the  wing  motion;  and 
the  microprocessor  linked  to  the  servomotor  drive 
system  that  will  eventually  translate  the  laser 
scanning  system  (In  the  direction  of  the  flow)  In 
order  to  allow  three-dimensional  sweeps  of  the 


flow  structure. 


Illumination  of  the  flow  structure  is  achieved 
by  a  two-watt  Argon-ion  laser  (with  a  multi-line 
mirror  holder)  impinging  upon  a  system  of  mirrors. 
This  system  allows  a  rotating  mirror  to  sweep  the 
laser  at  an  arbitrary  angle  with  respect  to  the 
freestream  while  the  entire  optical  train  is 
translated  in  the  freestream  direction.  Alternately, 
diffuse  illumination  can  be  employed;  it  involves 
stroboscopic  lights  having  a  power  of  90  watts,  a 
flash  duration  of  ten  microseconds,  and  a  repetition 
frequency  of  120  Hertz.  Or,  a  higher  powered  strobe 
having  an  intensity  of  8.3  Joules  at  120  cycles  per 
second  was  employed.  Recording  of  the  visualization 
images  was  achieved  with  either  one  or  two  video 
cameras  having  a  zoom  lens  capability.  Each  of  the 
two  video  cameras  could  be  connected  to  the  main¬ 
frame  of  the  video  system  and  the  corresponding 
images  displayed  simultaneously  using  the  split¬ 
screen  capability  of  the  mainframe. 


Figures  la  and  lb  show  simplified  overviews  of 
the  two  basic  types  of  image  acquisition:  a  duai- 
view  method  (Figure  la);  and  a  single-view  method 
(Figure  lb).  In  concept,  each  of  these  techniques 
provides  the  capability  for  tracking  visualization 
markers  in  three-dimensional  space. 


For  the  dual  image  method  of  Figure  la, 
stroboscopic  lighting  is  employed.  The  lights  are 
angled  with  respect  to  the  surface  of  the  bubble 
markers  in  order  to  provide  optimum  reflection  when 
viewed  through  the  cameras.  Mirror  rrangenents 
located  external  to  the  test  section,  as  well  as 
sufficiently  far  downstream  of  the  region  of 
interest  so  as  to  preclude  interference,  were 
employed  to  facilitate  viewing  the  flow  structure 
from  various  perspectives.  Various  locations  and 
orientations  of  25  micron  platinum  wires  allow 
generation  of  hydrogen  bubbles.  In  the  arrangement 
shown  in  Figure  la,  the  wire  flies  at  the  same 
velocity  as  the  apex  of  the  wing;  this  is 
accomplished  by  mounting  the  wire  between  two  probe 
supports  that  are  rigidly  attached  to  the  pitching 
axis  of  the  wing.  The  orientation  of  the  wire  is 
adjusted  to  allow  examination  of  either:  the  core 
flow,  by  ensuring  continuous  injection  of  marker 
at  or  near  the  apex  of  the  wing;  or,  the  exterior 
flow,  by  injecting  markers  from  a  wire  in  the 
plane  of  the  wing  at  the  apex  (Rockwell  et  al2l*>^5). 


For  the  single  Image  method  of  Figure  lb,  the 
principle  is  to  employ  the  localized  injection  of 
bubble  markers  of  defined  width  and  length  in 
conjunction  with  a  laser  sheet  (or  sheets)  of 
finite  thickness  in  order  to  determine  the  velocity 
field.  As  opposed  to  the  previous  technique,  only 
a  single  video  camera  is  employed  with  various 
mirror  arrangements  in  order  to  obtain  the  desired 
cross-sectional  view  of  the  flow  image.  This  image 
is  produced  by  laser-induced  reflection  from  the 
hydrogen  bubble  markers.  By  use  of  an  optically 
transparent  wing.  It  Is  possible  to  view  both  sides 
of  the  wing  simultaneously.  Bubble  markers  are 
generated  upstream  of  the  wing  from  a  single 
segmented  wire,  arrangements  of  them,  or  a  wire  grid 
as  schematically  portrayed  in  Figure  lb.  The 
basic  concept  of  this  technique  is  that  by 
generating  bubble  markers  of  defined  frequency  and 
duration  and  employing  a  laser  sheet  of  finite 
thickness,  estimates  of  the  velocity  field  can  be 
made  from  a  single  view. 
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Figure  la:  Dual  image  technique  for  tracking 
three-dimensional  flow  structure. 
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Figure  lb:  Single  Image  technique  for  determining 
two-  and  three-dimensional  flow 
structure. 


In  order  to  facilitate  quantitative 
interpretation  of  flow  visualization  images,  a 
variety  of  image  processing  techniques  have  been 
developed.  For  overviews  of  these,  the  reader  is 
referred  to  Rockwell  et  and  Smith^°.  In 

short,  quantitative  techniques  for  characterizing 
the  unsteady  process  of  bursting  in  a  turbulent 
boundary  layer  (Lu  and  Smith^)  and  vorticity  in 
an  unstable  wake  flow  (Lusseyran  and  Rockwell^”) 
have  been  employed.  In  addition,  efforts  have  been 
focussed  on  the  following  techniques:  reconstruc¬ 
tion  of  three-dimensional  flow  structure  using  a 
phase-reference  technique  (Ongoren  et  al^9)  ; 
ensemble-averaging  between  images  and  correlation 
within  or  between  images  (Kerstens  and  Rockwell 3®); 
and  Fourier  discriptor  techniques  for  pattern 
recognition  (Gumas  and  Rockwell^). 


QUANTITATIVE  FLOW  VISUALIZATION 

Flow  past  a  stationary  delta  wing  having  a 
sweep  angle  X  »  45°  and  inclined  at  an  angle  of 
attack  a  »  15°  with  respect  to  the  freestream 
produces  the  flow  structure  given  in  Figure  2a; 
the  photos  show  only  a  portion  of  the  flow  field. 

In  this  case,  the  wing  is  made  of  brass,  requiring 
modification  of  the  single  image  technique  of 
Figure  lb;  use  of  an  appropriate  mirror  arrangement 
allowed  simultaneous  views  of  the  upper  and  lower 
sides  of  the  wing.  Bubbles  were  generated 
continuously  from  the  grid  located  upstream  of  the 
wing,  and  the  laser  was  translated  to  various 
distances  x  from  the  apex  in  order  to  provide  the 
desired  images.  It  should  be  emphasized  here,  and 
in  subsequent  photos  of  this  type,  that  the 
thickness  of  the  bubble  sheets  forming  the  grid 
pattern  is  of  the  order  of  25  microns.  Use  of 
high  gain  on  the  video  system,  for  purposes  of 
illustration,  greatly  exaggerates  the  actual 
thickness.  In  this  sequence  of  photos,  the 
direction  of  circulation  about  the  wing  is  readily 
apparent.  Moreover,  it  clearly  increases  in  the 
streamwise  direction. 

By  tracking  each  of  these  nodes  of  the  bubble 
grid  in  three-dimensional  space,  and  making  use  of 
photos  in  addition  to  those  of  Figure  2a,  it  is 
possible  to  determine  the  three-dimensional 
streamline  pattern  of  the  flow.  This  is 
accomplished  by  digitizing  the  node  locations,  then 
using  CROSS  THREE-D  of  the  Unigraphics  CAO  software 
system.  The  three-dimensional  pattern  is  displayed 
on  a  HcAuto  0-100  C  color  terminal  allowing 
rotation  of  the  three-dimensional  image  in  real 
time.  Figures  2b  through  2d  show  selected  results. 
In  Figure  2b,  an  end  view  of  the  streamline  pattern 
is  illustrated.  Figure  2c  shows  an  end  view  of 
stream  surfaces,  formed  by  connecting  streamlines 
whose  upstream  positions  are  located  in  defined 
vertical  planes  prior  to  encounter  with  the  wing. 
Figure  2d  shows  the  same  stream  surfaces  from  a 
different  perspective.  Various  combinations  of 
horizontal  surfaces,  not  shown  here,  were  also 
constructed.  With  this  type  of  streamline 
information.  It  is  possible,  In  concept,  to  proceed 
with  calculation  of  the  corresponding  velocity  field 
and  circulation.  Occurrence  of  flow  separation 
along  the  upper  surface  of  the  wing  (see  photo 
corresponding  to  x/C  *  0.8  In  Figure  2a)  produces 
inaccuracies  of  the  streamline  pattern  in  and 
immediately  around  that  region.  Moreover,  there 
is  inadequate  spatial  resolution  there  that 
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Figure  2:  (a)  End  view  of  distorted  grid  of  hydrogen  bubble  lines  (actual  thickness  of  25  microns  greatly 

exaggerated  by  high  video  gain);  (b)  end  view  of  streamline  pattern  obtained  by  tracking  nodes 
of  distorted  grid  In  three-dimensional  space  and  connecting  them  with  CR0SS-3-D  CAD  technique; 
(c)  end  view  of  stream  surfaces  obtained  from  connecting  streamlines  In  initially  vertical 
plane  before  encounter  with  the  wing;  and  (d)  an  additional  view  of  (c)  from  an  elevated 
perspective. 


precludes  calculation  of  the  velocity  and  vorticity 
distributions.  This  technique  for  determining  the 
flow  structure  exterior  to  the  vortex  core  is 
currently  being  extended  to  unsteady  flows  past 
oscillating  delta  wings/ 

In  order  to  determine  the  streamline  pattern 
and  velocity  field  of  the  vortex  core,  it  is  most 
effective  to  generate  markers  locally  from  a 
single  wire  or  an  arrangement  of  wires  located 
within  the  core  itself.  The  photo  of  Figure  3a 
shows  representative  markers  from  a  single, 
segmented  bubble  wire  located  in  the  core  of  a 
vortex  on  a  delta  wing  having  a  sweep  angle 
>  =  75°,  an  angle  of  attack  a  ■  20°,  and  at  a 
streamwise  location  x/C  ■  0.87.  In  the  case  of 
steady  flow,  translation  of  the  laser  sheet  in 


the  streamwise  direction  provides  a  family  of 
marker  locations  in  three-dimensional  space  which, 
in  tur-  par.  be  connected  to  give  the  streamline 
pattern  of  the  flow.  Note  the  distortion  in  the 
marker  pattern  on  the  left-side  of  the  photo  as 
well  as  to  the  right  of  the  apparent  center  of  the 
vortex.  These  distortions  are  caused  by  the  shear 
layer  that  originally  sepsrated  from  the  edge  of 
the  wing.  Averaging  provides  the  streamline 
pattern  shown  in  Figures  3b  through  3d.  In 
Figure  3b.  the  view  is  in  the  upstream  direction 
along  the  axis  of  the  center  of  the  vortex. 

Figures  3c  and  3d  show  this  streamline  pattern  at 
different  orientations;  in  both  of  these  cases, 
both  the  streamlines  and  lines  corresponding  to 
constant  streamwise  coordinate  are  displayed. 


(c) 


Figure  3:  (a)  Elements  of  bubbles  generated  from  two  different  segmented  wires  located  in  vortex  core; 

(b)  end  view  of  streamlines  obtained  by  translating  laser  sheet  In  streamwise  direction  and 
tracking  bubble  elements  in  three-dimensional  space;  (c)  pattern  of  (b)  from  another 
perspective  including  lines  representing  constant  values  of  streamwise  coordinate;  and 
(d)  view  of  (c)  from  another  perspective. 
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Figure  it:  (a)  End  view  of  delta  wing  segments 
oscillating  about  a  common  axis  in 
harmonic  version  of  clap  and  fling 
maneuver,  with  top  photo  taken  at 
x/C  «  0.7  and  bottom  photo  x/C  “0.8 
and  a  time  delay  between  photos 
corresponding  to  the  transit  time  of  the 
freestream  between  two  stations  (bubble 
sheets  have  a  thickness  of  25  microns, 
greatly  exaggerated  here  due  to  video 
gain);  (b)  excerpt  of  upper  half  of 
velocity  field. 

Using  the  general  types  of  techniques  shown 
in  conjunction  with  Figures  2  and  3i  It  is  possible 
to  determine  directly  the  velocity  field  by  pulsing 
the  voltage  applied  to  the  bubble  wires,  such  that 


discrete  markers  are  generated.  Such  methods  can 
be  applied  to  both  steady  and  unsteady  flows,  and 
efforts  are  continuing  in  this  direction  in  our 
laboratory.  In  fact,  distributions  of  axial  and 
swirl  velocity  components  obtained  for  the 
conditions  described  in  conjunction  with  Figure  3a 
show  a  reasonable  approximation  to  the  expected 
d  i  stri but  ions. 

Determination  of  the  instantaneous  velocity 
field  in  the  cross-stream  plane  is  of  obvious 
importance  in  determining  the  circulation  for 
various  wing-flap  or  flap-flap  arrangements.  The 
photos  in  Figure  4a  correspond  '  r.  >  simple  harmonic 
version  of  the  clap  and  fling  motion  of  two  delta 
wing  segments  (or  "flaps")  oscillating  about  a 
common  axis,  corresponding  to  the  center  of  the 
corresponding  stationary  delta  wing.  The  top  photo 
was  taken  at  x/C  *  0.7  and  the  bottom  one  at  x/C  ■ 
0.8.  The  wing  segments  are  in  the  process  of 
closing,  and  the  time  shift  between  the  two  photos 
was  calculated  to  correspond  to  the  time  required 
for  a  defined  bubble  segment  to  travel  between 
x/C  »  0.7  to  0.8  in  the  freestream  region.  The 
axial  velocity  c'mwonent  at  any  location  in  the 
cross-section  can  be  estimated  by  pulsing  the 
front  of  the  bubble  grid  and  calculating  the 
"glow  time"  corresponding  to  passage  of  this  pulsed 
front  of  bubbles  through  the  laser  sheet  of  finite 
thickness.  Outside  the  vortex  core,  the  axial 
velocity  deviates  relatively  little  from  the 
freestream  value,  simplifying  the  tracking  process. 
For  the  present  purposes,  it  is  the  circulation 
about  the  vortex  core  that  is  of  primary  interest. 
The  velocity  field  corresponding  to  the  upper  half 
of  the  flow  is  shown  in  Figure  4b.  It  was  obtained 
by  averaging  the  node  displacements  about  the  plane 
of  symmetry  of  the  wing  segments,  then  carrying  out 
a  linear  interpolation.  Circulation  calculated 
using  two  different  circuits  differed  by  less  than 
31,  thereby  suggesting  that  the  circuit  for 
determining  the  circulation  was  exterior  to  the 
vortex  core.  The  dimensionless  circulation 
corresponding  to  the  pattern  of  Figure  4a  is  about 
one-half  that  generated  in  a  ramp-type  fling  of  a 
two-dimensional  wing  system  (Spedding  and 
Maxworthy32)  with  no  mean  flow  in  the  axial 
direction. 


PHASE  SHIFT  AN0  HYSTERESIS  OF  VORTEX  CORE 
DEVELOPMENT  AND  BREAKDOWN 

When  a  delta  wing  is  pitched  about  its 
trai 1 ing-edge  in  sinusoidal  motion,  examination  of 
the  flow  structure  over  the  range  of  reduced 
frequency  0.025  <  K  <  1.7,  and  mean  angle  of 
attack  5°  <  a  <  20  shows  that  there  are  two 
fundamental  types  of  vortex  development  during  the 
course  of  an  oscillation  cycle.  At  low  reduced 
frequencies,  the  core  development  is  of  the  form 
shown  in  the  left  column  of  photos  of  Figure  5, 
while  at  higher  reduced  frequency,  the  form  is  as 
shown  in  the  right  column.  Considering  first  the 
case  of  low  reduced  frequency  at  a  ■  9°+ ,  the  dye 
has  been  swept  from  the  mid-portion  of  the  wing  to 
a  region  along  Its  leading-edge.  At  a  ■  1 5°'t’ , 
onset  of  vortex  core  development  is  evident  in  the 
right  portion  of  the  photo,  and  at  a  “  I8°t  and 
19®+,  this  core  development  progresses  in  the 
upstream  direction  towards  the  apex. 
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delineated  in  forthcoming  writeups.  Among  them  is  Figure  5  shows  that  at  the  maximum  instantaneous 

the  fact  that  the  process  of  vortex  breakdown  shows  angles  of  attack  shown  therein,  the  core  maintains 

new  types  of  structure  beyond  that  of  the  tradi-  its  integrity  well  downstream  of  the  immediate 

tional  spiral  mode  observed  on  stationary  wings.  vicinity  of  the  apex.  A  particulary  dramatic 

These  breakdown  mechanisms  include  abrupt  onset  of  illustration  of  this  phase  shift  occurs  at  a  mean 

a  very  large-scale  spiral  mode;  or  a  column-type  angle  of  attack  5  »  20°  for  K  ■  0.13  (not  shown 

instability  that  can  supercede  the  classical,  here).  In  this  case,  the  maximum  streamwise  length 

small-scale  spiral  mode.  of  the  vortex  core  occurs  at  a  -  29°+.  The  issue 

arises  as  to  whether  this  sort  of  drastic  phase 

It  is  evident  from  the  photos  of  Figure  5  that  shift  can  be  effectively  employed  in  phased  motions 
there  is  suDstantial  phase  shift  between  the  of  control  surfaces  during  a  typical  maneuver, 

process  of  vortex  development  and  the  instantaneous 

position  of  the  wing.  For  the  stationary  wing.  These  phase  shifts  are  manifested  in 

the  vortex  core  is  completely  broken  down  at  the  hysteresis  loops  when  the  instantaneous  position 

apex  for  a  ■  20°.  In  contrast,  the  photo  of  of  the  vortex  breakdown  is  plotted  as  a  function 


Figure  6;  Hysteresis  loops  characterizing  Instantaneous  position  of  vortex  core  and  vortex  breakdown  as  a 
function  of  instantaneous  angle  of  attack. 
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of  instantaneous  angle  of  attack  a.  It  is  essential 
to  verify  the  character  of  such  loops  by  considering 
various  types  of  dye  injection  and  hydrogen  bubble 
flying  wire  methods  in  order  to  preclude  fallacious 
interpretation.  Figure  6  shows  that  even  at  a 
reduced  frequency  as  low  as  K  »  0.025.  there  is 
substantial  hysteresis  relative  to  the  characteris¬ 
tic  of  the  stationary  wing.  The  form  and  sense  of 
the  hysteresis  loop  is  generally  the  same  for  higher 
values  of  K  until  K  »  0.28  is  attained.  At  K  *  0.28, 
the  form  of  the  hysteresis  loop  takes  on  a  figure 
eight  form,  and  at  K  =  0.76  and  higher  values  of  K, 
the  form  and  the  sense  of  the  hysteresis  loop  change 
drastically.  Whereas  the  hysteresis  loops 
corresponding  to  lower  values  of  K  are  in  the 
clockwise  direction,  the  sense  of  higher  frequency 
loops  is  in  the  counterclockwise  direction.  These 
differences  in  form  and  sense  of  the  hysteresis 
loops  are  directly  related  to  the  manner  in  which 
the  vortex  core  develops,  as  discussed  in  conjunc¬ 
tion  with  Figure  5.  At  K  =  0.025,  it  is  of  the 
low  frequency  type,  at  K  »  0.76  of  the  high 
frequency  type,  and  at  K  =  1.28  there  is  transition 
between  these  two  extremes. 

In  addition  to  the  primary  vortex  described 
in  the  foregoing,  there  also  exists  an  unsteady 
secondary  (counterrotating)  vortex  between  the 
primary  vortex  and  the  leading-edge  of  the  wing. 
Typical  hysteresis  loops  of  the  secondary  vortex, 
in  relation  to  their  primary  counterparts,  are 
given  in  Figure  7.  Over  the  range  of  reduced 
frequencies  considered,  the  following  trends, 
exemplified  in  Figure  7,  are  evident  for  the 
secondary  vortex:  the  same  form  (shape)  of 
hysteresis  loop  as  for  the  primary  vortex;  a  smaller 
value  of  maximum  breakdown  length  than  that  of  the 
primary  vortex;  and  a  phase  shift  with  respect  to 
the  primary  vortex,  i.e.  lag  of  breakdown  of  the 
primary  vortex  relative  to  that  of  the  secondary 
vortex. 

To  characterize  the  time-dependent  evolution 
of  the  vortex  core,  characteristic  axial  and  swirl 
velocities  are  required.  Figure  8  shows  the 
variation  of  the  axial  velocity  component  u  along 
the  core  centerline  as  a  function  of  distance  x 
from  the  apex.  For  reference,  the  case  K  »  0  is 


also  shown.  Remarkable  is  the  essentially  constant 
value  of  u  during  the  initial  stage  of  the  core 
developmen?  followed  by  an  abrupt  decrease  to  zero, 
corresponding  to  the  instantaneous  stagnation 
condition  at  the  onset  of  vortex  breakdown.  Not 
shown  here  are  normalized  velocity  variations 
within  the  deceleration  region  over  a  range  of  mean 
angle  of  attack  ci  and  reduced  frequency  K.  Such 
normalization  suggests  a  universal  velocity 
variation  for  the  unsteady  core  that  is  essentially 
coincident  with  that  corresponding  to  the  stationary 
wing . 

Under  current  consideration  is  a  comparison 
of  various  functional  forms  of  the  wing  forcing, 
e.g.,  sinusoidal  vs.  ramp  forcing.  This  comparison 
will  allow  further  insight  into  the  mechanisms  of 
phase  shift  of  vortex  development  and  breakdown 
with  respect  to  the  wing  history. 


PERTURBATIONS  OF  VORTEX  CORE 

Perturbations  applied  to  a  swept  leading-edge 
can  provide  further  insight  into  the  possibility 
of  vortex  control  by  various  means  of  alteration 
of  boundary  conditions  at  the  leading -edge  in  the 
form  of  very  small  flaps,  pulsed  blowing,  etc. 

In  determining  the  appropriate  excitation  frequency 
f  ,  consideration  should  be  given  to:  the 
frequency  of  vortex  breakdown  (Leibovich*®) ;  and 
the  inherent  instability  frequency  of  the  shear 
layer  separating  from  the  edge  (Pierce^1*;  Payne  and 
Neison^S;  Gad-el-Hak  and  B1ackwelder^°) .  Here,  the 
case  of  a  delta  wing  having  a  sweep  angle  \  «  75°, 
at  an  angle  of  attack  o  *  35°  i  10°,  will  be 
addressed.  Oscillations  were  in  the  pitching  mode 
about  the  trail ing-edge;  the  amplitude  of  the  apex 
motion  was  0.0175  C  and  the  dimensionless  frequency 
of  oscillation  was  f  r  /u  «  0.027,  which  corre¬ 
sponds  to  the  natural  Breakdown  frequency  of  the 
core  on  the  stationary  wing.  These  parameters  apply 
to  Figures  9  through  14.  In  Figures  9  through  13, 
data  was  acquired  at  x/C  •  0.38,  upstream  of  vortex 
breakdown.  In  Figure  14,  the  station  x/C  »  0.883 
was  considered. 

The  mean  distributions  of  the  axial  and  swirl 
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Figure  7:  Representative  hysteresis  loops  of  secondjry  (counter-rotating)  vortex  In  comparison  with  Its 
corresponding  primary  vortex. 
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velocity,  designated  by  0  and  w  are  given  in 
Figure  9-  The  forms  of  these  distributions  are 
similar  to  those  corresponding  to  the  stationary 
wing  (f  »  0)  .  Particularly  important  is  the 
asymmetry  of  chese  velocity  distributions  with 
respect  to  the  center  of  the  vortex  core. 

Asymmetry  arises  from  mutual  induction  by  the 
vortices  from  both  of  the  leading-edges  of  the 
wing,  as  well  as  from  distortion  of  the  unsteady 
shear  layer  shed  from  the  edges  of  the  wing. 
Clearly,  theoretical  formulations  that  assume 
symmetry  of  the  velocity  and  vorticity  field  of 
the  vortex  core  must  be  interpreted  with  care. 
Calculation  of  criticality  conditions  and 
instability  criteria  for  the  vortex  core  should 
account  for  the  nonax i symmet r ic  mean  vorticity 
field.  Further  studies  of  the  mean  velocity  field 

K 
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•  0.28 
v  0.50 
a  0.76 
O  1.15 
9  1 .60 
a  i -91* 


5  -  10° 
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Figure  9:  Mean  axial  and  swirl  velocity  components 
on  a  delta  wing  (sweep  angle  A  »  75°)  at 
a  mean  angle  of  attack  5  »  35°  undergoing 
oscillations  of  magnitude  Aa  »  1°  at  a 
frequency  ferc/uc  *  0.027.  Velocities 
determined  at  x/C  ■  0.38. 


Figure  8:  Variation  of  instantaneous  axial 

velocity  along  centerline  of  vortex 
core  as  a  function  of  distance  from  the 
apex  of  wing  for  a  range  of  reduced 
frequency. 


are  focussing  on  determination  of  the  u  and  w 
distributions  for  forcing  at  fundamental,  as  well 
as  sub-  and  superharmonics  of  the  inherent 
instability  of  the  shear  layer  separating  from  the 
edge  and  the  vortex  breakdown  frequency. 

Rigorously  speaking,  one  should  consider  the 
entire  cross-section  of  the  flow,  and  not  simply 
the  velocity  distribution  through  the  center  of  the 
vortex  core,  as  suggested  in  Figure  9-  As 
illustrated  in  Figures  10  and  10b,  the  distribution 
of  the  axial  velocity  component  u  across  the  flow 
shows  further  features  of  the  asymmetry.  These 
three-dimensional  plots  extend  over  elevations 
from  the  wal I  of  z/s  ■  0. I  to  z/s  »  1.1,  in  which 
z  is  normal  to  the  wing  surface  and  s  is  the  local 
semi -span. 

The  amplitude  and  phase  distributions  of  the 
axial  fluctuation  component  u  at  the  forcing 
frequency  f  are  given  in  Figure  II.  These 
distributions  were  taken  at  the  same  elevation 
z  ■  z.  as  for  the  mean  velocities  of  Figure  9.  The 
distributions  of  amplitude  and  phase  at  a  spanwise 
location  near  the  edge  of  the  wing  are  character¬ 
istic  of  those  observed  in  unstable  shear  flows: 
two  adjacent  peaks  of  unequal  amplitude  with  a 
phase  shift  of  about  ir  across  them.  Further 
inboard,  near  the  apparent  center  of  the  vortex, 
there  occur  two  additional  peaks  with  a  gradual 
phase  shift  of  about  it  across  them.  These 
amplitude  and  phase  distributions  extending  across 
the  vortex  through  the  separating  shear  layer  form 
what  might  be  termed  a  composite  eigenfunction;  it 
embodies  both  the  eigenfunction  of  the  separating 
shear  layer  as  well  as  that  of  the  integrated 
history  of  the  vortex  of  the  core.  A  clear  picture 
of  this  composite  eigenfunction  requires  similar 
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(b) 


Figure  10:  (a)  Three-dimensional  representation 

of  the  distribution  of  the  mean  velocity 
component  u  over  the  cross-section  of 
the  flow,  (b)  Same  as  (a)  but  viewed 
from  different  perspective. 

information  for  the  swirl  component  of  velocity; 
it  is  currently  being  acquired. 

Contours  of  constant  magnitude  of  the 
fluctuating  velocity  u  over  the  flow  cross-section 
are  given  in  Figure  12.  The  location  of  the  star  * 
corresponds  approximately  to  the  apparent  center 
of  the  mean  vortex.  The  region  of  high  intensity 
near  the  wall  just  inboard  of  the  edge  is 
tentatively  attributed  to  the  secondary  (counter¬ 
rotating)  vortex  in  that  region.  Above  the  edge, 
the  elongated  regions  of  high  intensity  are 
clearly  due  to  the  separated  shear  layer.  As  the 
shear  layer  evolves  in  the  direction  of  the  mean 
swirl,  there  is  a  region  of  large  fluctuation 
amplitude,  located  above  the  vortex  center,  whose 
leading  portion  Is  distorted  In  the  direction  of 
the  swirl.  Just  below  the  apparent  center  of  the 


Figure  II:  Distribution  of  axial  component  of 

fluctuating  velocity  u  and  its  phase 
S.  through  the  center  of  the  mean 
vertex. 

vortex,  there  is  yet  another  peak  of  u.  The 
structure  suggested  by  the  contours  of  Figure  12, 
namely  that  of  an  unstable  three-dimensional  shear 
layer,  continuously  feeding  into  the  inner  portion 
of  the  vortex,  is  in  accord  with  flow  visualization 
(not  shown  here). 

Up  to  this  point  we  have  addressed  amplitude 
and  phase  variations  only  at  the  excitation 
component  f  .  However,  due  to  nonlinearity  of  the 
flow  structure,  other  spectral  components  will  be 
present.  Detailed  spectra  have  been  acquired  for 
a  range  of  excitation  conditions,  and  representa¬ 
tive  ones  are  shown  in  Figure  12.  In  Figure  13a, 
the  spectra  were  taken  at  the  location  of  the 
large  amplitude  peak  at  y/s  *  0.90  in  Figure  II  and 
in  Figure  13b  at  the  location  of  the  peak  at  y/s  » 
0.45.  As  evidenced  In  Figures  13a, b  and  in  other 
spectra  not  shown  here,  the  organized,  multiple- 
peaked  spectral  content  persists  across  the  entire 
extent  of  the  vortex  core  including  the  separating 
shear  layer.  Remarkable  is  the  highly  ordered 
nature  of  these  spectra  showing  as  many  as  nine 
well-defined  peaks  over  the  frequency  range.  Each 
of  these  peaks  corresponds  to  a  multiple  of  the 
fundamental  spectral  component,  representing  the 
forcing  frequency  f  .  The  possible  interpretation 
of  these  ordered  spectra  In  terms  of  concepts  of 
nonlinear  Interaction  is  underway.  There  are,  of 
course,  many  issues  related  to  the  spectral 
evolution  of  the  flow  in  the  streamwise  direction. 
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Figure  12:  Contours  of  constant  magnitude  of  u 
over  the  flow  cross-section. 

In  particular,  the  nature  of  the  spectral  energy 
transfer  as  the  vortex  enters  and  passes  through 
the  breakdown  region  is  of  particular  importance. 

External  forcing  can  alter  the  mean  flow 
character  and  thereby  the  vortex  breakdown 
location  in  several  interesting  fashions,  as  will 
be  addressed  in  forthcoming  writeups.  We  address 
here  one  such  observation.  The  vortex  on  the 
opposite  side  of  the  wing  (relative  to  the  one  under 
consideration)  is  stabilized  such  that  its  break¬ 
down  occurs  beyond  the  tra i 1 ing-edge.  For 
excitation  at  frequency  f  equal  to  that  of  the 
natural  breakdown  frequency  f*,  Figure  Ik  compares 
representative  mean  and  fluctuating  velocity 
distributions  without  and  with  forcing  at  the  same 
streamwise  location  x/C  *  O.883.  Without  forcing, 
there  is  clearly  vortex  breakdown:  the  mean 
velocity  distribution  u  shows  a  wake-like  profile 
and  the  velocity  goes  through  zero  to  a  negative 
value  as  the  vortex  centerline  is  approached.  The  0 
distribution  of  the  orqanized  component  at  the 
breakdown  frequency  f*  is  also  illustrated.  On 
the  other  hand,  in  the  case  of  forcing  at  the 
naturaj  frequency  of  vortex  breakdown,  i.e.  at 
fe  *  fo,  vortex  breakdown  does  not  occur.  The 
distance  that  the  breakdown  of  the  forced  vortex 
is  delayed,  relative  to  the  unforced  vortex, 
corresponds  to  about  twenty  percent  of  the  chord  C. 
The  mean  velocity  u  of  the  forced  vortex  shows 
the  same  general  form  as  that  of  Figure  9-  The 
magnitudes  of  the  spectral  components  u(fe)  and 
u(2fe)  are  comparable.  Determination  of  the 
possibility  of  a  simultaneous  or  Intermittent 
presence  of  symmetric  and  spiral  modes  of 
instability  requires  corresponding  phase 
distributions;  this  issue  is  currently  under 
consideration. 


CONCLUDING  REMARKS 

Active  control  of  the  flow  structure  of 
leading-edge  vortices  on  a  delta  wing  and  delta 
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Figure  13:  Representative  spectra  of  the  velocity 
fluctuation  u. 

wing  segments  through  both  small  and  large 
(displacement)  amplitude  pitching  motion  is 
attainable  If  the  excitation  frequency  is  properly 
selected.  Such  excitation  should  consider  the 
following  characteristic  frequencies:  the 
fundamental  frequency  (or  its  subharmonics)  of  the 
unstable  shear  layer  separating  from  the  edge  of 
the  wing;  the  frequency  of  vortex  breakdown;  and 
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Figure  14:  Comparison  of  representative  mean  and  fluctuating  velocity  distributions  of  the  unforced  and 
forced  vortices  at  the  same  streamwise  location  x/C  ■  O.883. 


the  characteristic  frequency  of  the  separation 
zone  downstream  of  the  onset  of  vortex  breakdown. 
Depending  upon  the  frequency  selected,  there  can 
occur  a  variety  of  interesting  modifications  of  the 
leading-edge  vortices.  Among  them  are  alteration 
of:  the  mean  velocity  field  of  the  vortex;  the 
mechanism  of  vortex  breakdown;  and  the  phase  shift 
of  the  instantaneous  vortex  breakdown  position 
relative  to  the  instantaneous  angle  of  attack. 

Development  of  an  integrated  active  control 
and  quantitative  measurement  system  is  aimed 
towards  relating  the  instantaneous  flow  structure 
to  the  instantaneous  wing  motion.  80th  global 
and  local  characterization,  exterior  to  and  within 
the  vortex  core,  are  insightful  and  constitute 
the  thrust  of  further  efforts  in  this  direction. 
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SUMMARY 

Unsteady  two-  and  three-dimensional  flow  structure  at  leading-  and  trai 1 Ing-edges  of  bodies  can  be 
characterized  effectively  using  recently  developed  techniques  for  acquisition  and  Interpretation  of  flow 
visualization.  The  techniques  addressed  herein  Include:  flow  tmage-surface  pressure  correlations; 
three-dimensional  reconstruction  of  flow  structure  from  flow  Images;  and  Interactive  interpretation  of 
flow  images  with  theoretical  simulations.  These  techniques  can  be  employed  In  conjunction  with:  visual 
correlation  and  ensemble-averaging,  both  within  a  given  Image  and  between  images;  recognition  of  patterns 
of  flow  structure  from  Images;  and  estimates  of  velocity  eigenfunctions  from  Images. 


1.  INTRODUCTION 

There  are  a  wide  variety  of  methods  for  visualizing  and  Interpreting  flow  phenomena.  The  overviews 
of  Enrich*  and  rierzk’rch  et  al*  assess  a  range  of  techniques  Including,  for  example:  smoke  and  dye 
injection;  laser-induced  fluorescence;  and  density  gradients  Induced  In  mixing  of  dissimilar  fluids. 

These  methods  have  provided  valuable  Insight  Into  a  number  of  features  of  the  flow  structure. 

Flaw  visualization  using  water  as  a  working  fluid  has  several  advantages.  Perhaps  the  foremost  Is 
the  ease  of  marking  and  tracking  fluid  elements  at  relatively  long  time  scales,  l.e.  low  frequencies,  in 
comparison  with  the  corresponding  scales  In  a  gaseous  medium.  If  Instantaneous,  as  opposed  to  time- 
averaged,  surface  pressure,  lift,  etc.  In  conjunction  with  the  visualized  flow  structure  Is  desired,  then 
the  time  and  length  scales  offered  by  a  water  facility  are  highly  advantageous.  Although  concurrent, 
quantitative  measurement  of  unsteady  velocity  and  pressure  may  be  more  delicate  and  laborious,  judicious 
choice  of  the  most  physically  meaningful  locations  for  measurement  of  the  flow  field,  Identified  with  the 
aid  of  flow  visualization,  can  reduce  the  amount  of  effort  In  this  direction. 

Interpretation  of  visualized  steady  flows  Is  relatively  straightforward;  however  this  Is  not  the 
case  for  unsteady  flows,  especially  for  perturbed  shear  flows.  It  Is  well  known  from  linear  stability 
thaory  that  there  exist  drastic  gradients  of  amplitude  and  phase  of  fluctuating  velocity  and  pressure 
across  an  unsteady  shear  flow.  As  a  consequence,  the  apparent,  visualized  flow  structure  may  not  be 
compatible  with  the  actual,  underlying  vortlclty  field  (Hama3;  Lusseyran  and  Rockwell1*).  Clearly, 
considerable  care  must  be  exercised  In  Interpreting  unsteady  phenomena.  By  proper  choice  of  visualization 
technique,  and  employment  of  simple  theoretical  simulations,  fallacies  of  Interpretation  can  be  precluded. 

In  the  series  of  Investigations  described  herein,  we  address  techniques  for  Implementing  and 
Interpreting  a  particular  class  of  flow  visualization  methods:  generation  and  tracking  of  lines  of  elements 
marked  at  specified  times,  l.e.  timelines.  Emphasis  is  on:  correlating  the  Instantaneous,  visualized  flow 
structure  with  Instantaneous  surface  pressure;  reconstruction  of  three-dimensional  flow  structure  from 
dual  views  of  the  flow  field;  and  Interpretation  of  the  flow  visualization  In  conjunction  with  basic 
classes  of  theoretical  simulation,  in  addition,  possibilities  are  addressed  for:  visual  correlation  and 
en samb I e- averaging  of  visualization  images;  pattern  recognition  of  flow  structure  using  Images;  and 
evaluation  of  Images  for  determination  of  velocity  eigenfunction  and  vortlclty. 


2.  EXPERIMENTAL  TECHNIQUES 

All  experiments  were  carried  out  In  one  of  three  water  channel  systems,  custom  designed  for  study  of 
unsteady,  separated  flows.  The  cross-sectional  areas  of  the  test  sections  of  these  cham  .Is  ranged  from 
1.5  ft*  to  6  ft*.  In  all  cases,  the  test  sections  were  made  entirely  of  plextglas  or  glass  to  facilitate 
multl-dlmensional  views  of  the  flow  structure.  Moreover,  all  components  of  the  largest  channel,  including 
the  test  section  contraction.  Inlet  and  outlet  tanks,  and  all  attendant  piping  are  made  of  plexiglas  or 
PVC  In  order  to  preclude  difficulties  arising  from  system  corrosion. 

Figure  I  shows  an  experimental  arrangement  that  Is  generic  to  the  variety  of  techniques  employed  In 
the  fluid  mechanics  laboratory  at  Lehigh  University.  A  t>i«(.lnur»  m!.b,  or  i..  a.  .--.'ige.nar.i  of  wires, 
generates  hydrogen  bubble  timelines.  The  bubble  lines  or  sheets  are  then  tracked,  recorded,  and 
processed  by  employing:  video  cameras  synchronized  with  the  mainframe  of  the  video  system;  a  digitizing 
system;  and  an  Image  processing  system  located  In  the  computer-aided  design  (CAD)  laboratory. 

Using  a  custom-designed  function  generator,  It  Is  possible  to  generate  hydrogen  bubble  rimelines  of 
desired  frequency  and  width.  Two  stroboscopic  (90  watt)  lights,  Interfaced  with  the  mainframe  of  the 
video  system.  Illuminate  these  bubble  timelines. 
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Figure  I:  Schematic  of  technique  for  acquiring  and  processing  flow  visualization  using  two-camera,  split¬ 
screen  video  system,  digitizer,  and  computer-aided  design  system.  Technique  shown  Is  generic  to  the  variety 
of  arrangements  employed  In  the  Lehigh  University  fluid  mechanics  laboratories. 


Two  simultaneous  images,  corresponding  to  two  arbitrary  views  of  the  flow  field,  can  be  displayed  on 
the  split-screen  video  system.  A  variety  of  mirror  arrangements  provides  views  of  the  flow  field  from  a 
number  of  perspectives.  These  mirrors  are  located  either  external  to  the  transparent  test  section,  or 
within  the  flow  field  at  a  location  sufficiently  far  downstream  of  the  region  of  Interest. 

Video  images  can  be  digitized  using  a  Colorado  video  (model  270A)  digitizer  or  a  digitizing  tablet. 
Images  obtained  from  the  digitizer  are  put  on  file  In  the  CAD  laboratory,  allowing  application  of  a  range 
of  Image  pre-processing  and  processing  techniques. 

Pre-processing  of  digitized  Images  enhances  the  definition  of  the  bubble  timeline  patterns.  These 
processing  operations  include:  a  "band-pass"  spatial  filter;  a  thresholding  operation;  a  thinning 
algorithm;  a  two-dimensional,  low-pass  filter;  and  other  related  techniques,  which  are  described  In  detail 
by  Gumas  and  Rockwell.  As  a  follow-on  to  this  image  pre-processing.  It  Is  possible  to  employ  a  curve¬ 
tracking  technique  which,  in  essence,  transforms  timelines  having  multiple  grey  levels  to  continuous 
well-defined  timelines  in  binary  form.  Gumas®  describes  the  curve-tracking  algorithm,  which  Is  based  on 
the  concept  of  linear  predictive  coding. 


3.  QUANTITATIVE  INTERPRETATION  OF  VISUALIZATION  IMAGES 

.  e 

Figure  2  shows  a  typical  timeline  visualization  photo,  corresponding  to  anti  symmetrica  I  vortex 
shedding  from  a  cylinder  oscillating  In  the  streamwlse  direction  (Ongoren  and  Rockwell?).  Also 
illustrated  is  a  representation  of  the  digitized  timelines,  providing  a  basis  for  quantitative  Inter¬ 
pretation  of  the  flow  structure.  A  timeline  is,  by  definition,  made  up  of  all  fluid  elements  originally 
marked  at  the  same  time  at  the  upstream  location  of  the  bubble  wire.  Taking  these  timelines  as 
reprasentat iva,  we  address  possibilities  for:  visual  correlation  and  ensemble-averaging  of  Images;  pattern 
recognition  using  visualization  images;  and  image  evaluation  for  velocity  eigenfunction  and  vortlclty. 
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Figure  2:  Typical  flow  visualization  photograph  representing  vortex  shedding  In  antlsymmetrlcal  mode 
arising  from  streamwlse  perturbations  of  cylinder.  Also  shown  is  digitized  representation  of  upper  half 
of  timeline  pattern  (Ongoren  and  Rockwell';  Kerstens  and  Rockwell**). 


Visual  correlation  and  ensemble-averaging  of  Images 

Although  correlations  of  velocity  and  pressure  at  a  point,  or  between  two  points.  In  a  flow  field 
have  received  considerable  attention  over  the  past  few  decades,  relatively  little  effort  has  been  devoted 
to  defining  a  more  globally-oriented  correlation  based  on  flow  visualization.  In  doing  so.  It  Is  necessary 
to  define  quantitatively  the  elemental  portions  of  a  visualization  Image.  One  possibility  Is  to  define 
the  local  tangent  of  each  of  the  Individual  timelines  as  a  function  of  position  along  the  timeline.  This 
local  tangent  Is  termed  the  tangent  angle  function  t(x,y).  If  we  know  the  function  ♦  along  all  the 
timelines  In  an  Image,  we  have.  In  effect,  defined  the  Image  on  a  physically  meaningful  basis. 

The  obvious  approach  to  defining  a  correlation  between  two  points  A-l  and  B-l  on  two  different 
timelines  A  and  B  Is  to  consider  the  product  of  the  respective  tangent  angle  functions  +A-I+B-I-  However, 
as  shown  by  Kerstens  and  Rockwall”,  this  definition  Is  Inadequate,  as  It  cannot  distinguish  certain 
features  of  two  curves  oriented  arbitrarily  with  respect  to  each  other.  Rather,  It  Is  useful  to  employ 
a  correlation  function  based  on  the  difference  between  tangent  angle  functions  ^a-I~*B-1*  Employing  such 
a  definition,  It  Is  then  possible  to  define  a  global,  visual  correlation  between  two  timelines  by 
integrating  over  the  spatial  extent  of  the  normalized  timelines.  This  approach  can  lead  to  various  types 
of  auto-  and  cross-correlations  between  defined  regions  of  the  flow  field.  Ideally,  It  would  be  desirable 
to  carry  out  these  correlations  for  segmented  timelines,  where  a  timeline  Is  broken  Into  discrete  segments 
through  use  of  a  segmented  bubble  wire;  then,  by  tracking  these  segments  and  carrying  out  the  afore¬ 
mentioned  visual  correlation,  a  comparison  can  be  made  with  the  measured  vortlclty  field. 

The  visual  average  of  several  successive  Images  follows  from  the  spatial  averages  of  respective 
timelines  In  a  number  of  successive  Images  (Kerstens  and  Rockwell^).  In  selecting  the  Images  to  be 
averaged,  a  physically  meaningful  criterion  must  be  applied.  If  each  of  the  selected  Images  corresponds 
to  a  deflnnd  phase  of  an  unsteady  event  within  the  flow,  such  as  sn  extrema  displacement  of  a  body, 
than  the  averaging  Is  of  the  ensemble-averaging  type.  In  principle,  a  variety  of  sensors.  Including 
pressure  and  velocity  transducers,  can  be  employed  as  a  phase  reference  for  the  ensemble  averages; 
moreover,  an  Interesting  possibility  Involves  use  of  a  repetitive,  coherent  portion  of  the  Image  itself 
In  order  to  trigger  the  ensemble-averaging. 

Pattern  recognition  using  Images 

In  the  event  that  the  flow  Is  non-perlodlc.  It  Is  useful  to  determine  how  often  certain  large-scale 
features  of  the  flow-structure  occur.  Consequently,  It  Is  necessary  to  formulate  a  technique  for  rapidly 
identifying  these  structures  from  flow  Images.  If  a  library  of  basic  types  of  flow  patterns  can  be 
described,  then  comparison  of  visualized  flow  patterns  with  this  library  provides  a  basis  for  characterizing 
the  frequency  of  appearance  of  elemental  types  of  flow  structure.  This  approach  represents  a  type  of 
pattern  recognition. 

A  technique  employed  for  Identification  of  machine  parts  (Persoon  and  Fu^) ,  handprint  character 
recognition  (Granlund*”) ,  and  recognition  of  aircraft  silhouettes  (Wallace  and  Mitchell”)  Involves 
definition  of  a  library  of  Fourier  descriptors  for  basic  types  of  patterns.  The  Fourier  descriptor  Is 
simply  the  Fourier  transform  of  the  tangent  angle  function  b,  •*  previously  defined.  The  advantage  of 
the  Fourier  descriptor  technique  Is  that  the  curve  description  Is  essentially  Independent  of  orientation, 
position,  and  scale  of  the  curve.  Whereas  the  traditional  use  of  the  Fourier  descriptor  technique  Is  for 
closed  curves,  flow  Images  are  typically  made  up  of  open  curves  In  the  form  of  streakllnes,  pathllnes,  or 
timelines.  Gumas  and  Rockwell^  describe  Implementation  of  the  Fourier  descriptor  technique  for  open 
curves  of  the  timeline  type. 

Image  evaluation  for  velocity  eigenfunction  and  vortlclty 

The  displacement  of  a  timeline  over  a  given  Interval  of  time  Is  related  to  the  flow  velocity. 
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Consequently,  It  Is  possible  to  relate  the  time  of  flight  of  timelines  to  the  velocity  distribution 
across  the  flow.  In  fact,  lu  and  Smith*2  have  shown  that  the  higher  order  statistics  of  a  turbulent 
boundary  layer  can  be  determined  using  this  concept;  this  approach  can  lead  to  Insight  beyond  that 
attainable  with  traditional  point  measurements. 

A  primary  advantage  of  the  timeline  technique  Is  that  It  gives  Instantaneous  Information  across  the 
entire  flow;  It  therefore  leads  to  the  Instantaneous  streamwlse  velocity  component  across  the  flow. 
Consequently,  this  method  Is  well-suited  to  characterizing  flows  for  which  there  are  substantial  gradients 
of  velocity  amplitude  and  phase  In  the  cross-stream  direction.  That  Is,  the  local  deflection  of  a  time¬ 
line  will  be  In  accord  with  the  local  fluctuation  amplitude  and  Its  relative  phase.  Lusseyran  and 
Rockwel l”  have  demonstrated  that  this  technique  can  provide  a  reasonable  estimate  of  the  eigenfunction 
of  the  streansvlse  velocity  component.  Through  application  of  the  continuity  equation,  It  leads  to 
approximation  of  contours  of  constant  vortlclty  tn  an  unstable  wake  flow. 


4.  FLOW  IMAGE  -  SURFACE  PRESSURE  CORRELATION 

In  order  to  determine  the  underlying  mechanisms  that  produce  unsteady  pressure  fluctuations  along  a 
surface,  it  Is  desirable  to  correlate  the  Instantaneous,  visualized  flow  structure  with  the  Instantaneous 
surface  pressure.  One  possibility  Is  to  record  simultaneously  the  visualized  flow  and  the  time-dependent 
surface  pressure  from  a  large  array  of  surface-mounted  transducers.  If  the  flow  structure  Is  highly 
organized  then.  In  principle,  It  Is  possible  to  employ  a  single  transducer  switched  to  successive  locations 
along  the  surface.  8y  keeping  track  of  pressure  amplitude  and  phase,  the  Instantaneous  pressure  field  can 
be  reconstructed. 

In  acquiring  these  pressure  measurements,  It  Is  necessary  to  have  an  appropriate  phase  reference, 
which  can  take  the  form  of  a  second  pressure  transducer  signal,  or  a  hot  film  or  laser  signal, 
appropriately  located  at  the  outer  edge  of  the  mean  shear  layer.  Proper  location  of  this  reference  Is 
necessary  to  avoid  phase  gradients  associated  with  the  fluctuating  field,  as  well  as  high  frequency  noise 
from  the  flow.  In  linking  the  unsteady  surface  pressure  measurements  to  the  visualized  flow  field,  it  is 
necessary  to  determine,  at  a  given  Instant,  the  phase  shift  between  the  visualized  flow  struc*’ire  and  the 
instantaneous  pressure.  One  possible  technique  for  accomplishing  tills  Is  to  record  simultaneously  the 
flow  structure  and  the  oscilloscope  trace  of  the  reference  signal  on  a  split-screen  video  system. 

Obviously,  there  are  other  methods  of  Indirectly  determining  this  phase  link. 

Aside  from  the  advantages  of  visualization,  use  of  water  as  a  working  medium  allows  use  of  remote 
transducer  locations,  as  opposed  to  the  classical  surface-mounting  of  transducers.  The  characteristic 
frequencies  of  the  pressure  fluctuations  In  water  are  typically  very  low  -  on  the  order  of  I  Hz.  At 
these  low  frequencies,  the  amplitude  and  phase  distortion  arising  from  locating  the  transducer  away  from 
the  surface  are  Insignificant,  provided  the  connecting  line  between  the  surface  tap  and  the  transducer 
Is  properly  designed. 

A  means  of  visually  displaying  the  Instantaneous  flow  structure  and  surface  pressure  has  been 
developed  In  our  CAD  laboratory.  First,  the  flow  Images  are  digitized  and  processed.  Then  each  of  these 
Images  Is,  In  turn,  displayed  with  the  Instantaneous  surface  pressure  field.  This  field  Is  obtained  from 
cross-spectral  analysis  between  the  surface  pressure  at  a  given  location  and  an  appropriate  phase 
reference.  In  the  following,  we  discuss  cases  where  the  raw.  Instantaneous  visualization  Images,  without 
any  processing  or  enhancement,  are  shown  In  conjunction  with  the  Instantaneous  surface  pressure. 

Figure  3  shows  the  Interaction  of  an  Incident  vortex  with  an  elliptical  leading-edge.  In  this  case, 
the  vortex  Is  generated  In  a  mixing-layer  flow  arising  from  the  Inherent  Instability  of  the  shear  layer. 

As  a  consequence,  the  process  of  vortex-edge  Interaction  Is  harmonic.  Full  details  of  this  experiment, 
Including  the  circulation  and  degree  of  concentration  of  vortlclty  of  the  Incident  vortex.are  given  by 
Sohn  and  Rockwell***.  In  the  left  column  of  Figure  3,  the  vertical  hatched  lines  represent  the  amplitude 
envelope,  and  the  closely- spaced  hatched  lines  represent  the  Instantaneous  surface  pressure  fluctuations. 
Considering  the  visualization  In  the  right  column  of  photos  In  Figure  3,  It  Is  evident  that  the  Incident 
vortex  undergoes  substantial  distortion  as  It  encounters  the  leading-edge  region.  Concurrently,  there 
Is  local  flow  separation  from  the  leading  portion  of  the  edge,  giving  rise  to  a  small-scale  secondary 
vortex  that  becomes  nested  within  the  Incident  primary  vortex.  The  third  row  of  pressure  -  visual Izatlon 
shows  that  onset  of  flow  separation  Induces  a  local  negative  pressure;  the  fourth  and  fifth  rows  show  that 
this  negative  pressure  region  moves  downstream  as  the  formation  of  the  secondary  vortex  evolves.  It  Is 
clear,  especially  from  the  fourth  row  of  Figure  3,  that  there  are  very  large  gradients  of  Instantaneous 
pressure  In  the  vicinity  of  the  tip.  In  fact,  maximum  positive  pressure  occurs  at  the  tip  Itself,  and 
a  short  distance  downstream  of  that  location,  there  Is  a  maximum  of  negative  pressure. 

Figure  4  shows  the  Instantaneous  loading  and  corresponding  flow  structure  for  the  case  of  a  fully 
turbulent  boundary  layer  separating  from  an  actively-controlled  tral I Ing-edge.  In  this  case,  the 
displacement  perturbat Ions  of  the  tral I Ing-edge  are  about  two  orders  of  magnitude  smaller  than  Its 
thickness  Nevertheless,  the  surface  pressure  fluctuations  arising  from  the  edge  perturbations  overshado* 
those  of  the  beckground  turbulent  flow.  Since  the  first  harmonic  component  of  these  pressure  fluctuations 
Is  very  small,  we  can  represent  the  fluctuating  surface  pressure  as  approximately  sinusoidal,  and  employ 
the  edge  displacement  as  a  phase  reference  In  relating  the  Instantaneous  pressure  to  the  flow  structure. 

The  top  row  of  pressure-visualization  shows  the  case  of  excitation  at  a  frequency  fe  Just  below  the 
natural  shedding  frequency  f£,  and  the  middle  row  at  a  frequency  just  above  It.  From  the  visualization, 
we  note  that  the  Initially  shed  vortex  switches  from  the  bottom  to  the  top  corner  of  the  tral I Ing-edge 
as  the  excitation  frequency  Increases.  The  bottom  row  of  pressure-visualization  represents  the  case  of 
excitation  at  a  frequency  well  above  the  natural  shedding  frequency.  It  Is  apparent  that  the  Initial 
vortex  Is  still  shed  from  the  upper  surface.  In  fact,  only  when  the  edge  Is  excited  at  a  frequency  near 
that  of  the  natural  shedding  Is  a  switch  In  phasing  of  the  Initial ly-shed  vortex  possible. 

In  the  corresponding  pressure  plots,  the  solid  line  and  the  cross-hatched  regions  represent 
respectively  the  amplitude  envelope  and  the  instantaneous  value  of  the  pressure  amplitude  at  the  forcing 
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Figure  4:  Flow  structure  and  Instantaneous  surface  pressure  Induced  by  oscillating  tral I Ing-edge.  Approach 
boundary  layer  Is  fully  turbulent  and^pressure  field  Is  that  resulting  from  motion  of  edge.  Frequency  fe 
of  excitation  relative  to  frequency  f0  of  shedding  from  stationary  edge  Is:  fe/fo  “  0*87  (top  row);  1.02 
(middle  row);  and  2.0  (bottom  row)  (Staubll  and  Rockwell'*). 


frequency.  Moreover,  within  the  schematic  of  the  edge  are  shown  the  amplitude  envelopes  of  the  fluctuating 
force  Fe  at  the  excitation  frequency  f*  and  the  force  F0  at  the  self-excited  frequency  component  f0.  The 
vertical  bar  represents  the  Instantaneous  location  of  the  force  Fe.  Finally,  the  phase  clock  In  the  upper 
left  of  each  schematic  shows  that  the  force  Fe  lags  the  edge  displacement  ne  at  the  excitation  frequency 
below  the  natural  shedding  frequency  (top  diagram),  while  It  leads  at  a  frequency  above  natural  frequency 
(middle  diagram)  and  again  !«gs  at  high  excitation  frequency  (bottom  diagram).  For  the  conditions  of  the 
top  and  middle  diagrams,  the  amplitude  of  the  fluctuating  pressure  rapidly  Increases  as  the  corner  (B)  of 
the  edge  Is  approached;  it  goes  to  zero  at  the  middle  of  the  vertical  base  (BC).  In  the  bottom  diagram 
this  Is  no  longer  the  case;  due  to  the  high  excitation  frequency,  the  nonclrculatory  ("added-mass") 
contribution  to  the  surface  pressure  dominates  that  associated  with  shedding  of  vortlclty,  except  In  the 
base  region  BC. 

These  types  of  flow  structure  -  pressure  correlations  have  the  potential  for  providing  Insight  Into 
the  source  of  unsteady  lift  and  drag  acting  on  bodies;  they  are  currently  being  Interpreted  In  this 
context. 
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5.  THREE-D I HENS I ONAL  RECONSTRUCTION  FROM  FLOW  IMAGES 

Determination  of  the  Instantaneous  structure  of  time-dependent,  three-dimensional  flows  poses  a 
particular  challenge.  By  use  of  a  timeline  marking  technique,  It  Is  possible  to  reconstruct  the  three- 
dimensional  evolution  from  two-dimensional  Images  using  computer-aided  design  methods.  There  are 
basically  two  types  of  approaches:  a  dual-view  method,  where  the  three-dimensional  structure  Is 
determined  from  two  arbitrary  views  of  the  flow  field;  and  a  single-view  method.  Involving  acquisition 
of  phase-locked  Images  from  a  single  view  of  the  flow  field. 

Smith  and  Paxson'®  have  constructed  a  single  three-dimensional  surface  by  generating  timelines  from 
a  single  wire  oriented  In  the  spanwf se  direction,  and  using  a  two-camera  system  with  two  simultaneous 
views.  Their  study  provides  new  Insight  Into  the  character  of  a  turbulent  boundary  layer.  Ongoren  and 
Rockwell'7  have  generated  a  family  of  three-dimensional  surfaces  representing  the  wake  from  a  three- 
dimensional  body  by  employing  phase-locked  Images  of  timelines  from  a  single  view.  In  this  technique, 
the  body  was  subjected  to  controlled  oscillations.  Its  displacement  served  as  a  phase  reference,  thereby 
allowing  definition  of  phase-locked  Images  at  successive  spanwlse  planes.  In  the  following,  we  demonstrate 
the  dual-view  method  for  obtaining  three-dimensional  surfaces  representing  separated  flows  from  stationary 
and  oscillating  bodies. 

Figure  5  shows  time  sequences  of  end  and  plan  views  of  a  tip  vortex.  In  generating  these  Images,  a 
one  mil,  kinked  platinum  wire  was  employed.  By  pulsing  the  wire  at  a  desired  time,  and  over  a  defined 
Interval,  It  Is  possible  to  reveal  the  principal  features  of  the  flow  structure.  In  the  end  view  of 
Figure  5,  the  pulse  width  Is  relatively  long,  being  Initiated  just  prior  to  the  first  photo  and  terminated 
prior  to  the  last  photo.  The  plan  view  of  Figure  5  shows  that  the  free-stream  Is  moving  In  the  downward 
direction.  The  vortex  core  has  a  significantly  lower  axial  velocity  than  that  of  the  free-stream.  In 
these  end  and  plan  views,  it  Is,  In  principle,  possible  to  estimate  the  unsteady  velocity  field.  In 
addition,  as  depicted  at  the  bottom  of  Figure  5,  three-dimensional  surfaces  can  be  constructed  using  the 
CAD  system.  The  various  views  shown  therein  are  from  arbitrary  perspectives.  Real  time,  three-dimensional 
simulation  of  the  structure  of  this  tip  vortex  can  be  carried  out  either  on  the  Unigraphics  UG-2  terminal 
or  the  Evans  and  Sutherland  PS-300  terminal  In  the  Lehigh  University  CAD  laboratory. 

Another  class  of  three-dimensional  flows  that  can  be  characterized  effectively  by  the  dual  view 
technique  Is  that  generated  by  an  oscillating  delta  wing.  To  Illustrate  the  complexity  of  the  flow,  we 
first  consider  a  plan  view  of  a  delta  wing  oscillating  In  pitch  about  Its  trai I Ing-edge  (Figure  6). 

Since  the  tip  of  the  delta  wing  undergoes  large  excursions  during  an  oscillation  cycle,  only  limited 
information  can  be  obtained  by  examining  timeline  patterns  from  a  fixed  wire.  A  "flying  wire"  technique 
has  been  developed  (Atta  and  Rockwell'®);  the  platinum  wire  Is  aligned  In  the  spanwlse  (z)  direction, 
stretched  between  two  supports  located  well  away  from  the  wing,  and  fixed  to  the  mechanism  that  controls 
the  wing  motion.  Consequently,  the  wire  "files"  at  the  same  angular  velocity  as  the  wing,  and  provides 
an  Indication  of  the  degree  to  which  the  flow  entering  a  plane  parallel  to  the  wing  passes  below,  as 
opposed  to  above,  the  wing  surface.  For  the  photos  of  Figure  6,  the  wire  Is  at  the  tip.  Moreover,  the 
reduced  frequency  K  -  2irfC/U_  ”  k.O,  the  mean  angle  of  attack  a  »  35°  and  fluctuating  angle  of  attack 
Aa  -  10°.  From  quasi -steady  considerations,  one  expects  the  pattern  of  vortex  formation  to  be  the  same 
at  a  given  value  of  angle  of  attack  a.  Irrespective  of  whether  the  tip  of  the  wing  is  moving  in  the  upward 
or  downward  direction.  Comparison  of  photos  C  and  E,  taken  at  nearly  the  same  value  of  a,  shows  that 
the  flow  structure  Is  markedly  dissimilar.  Evidently,  there  Is  strong  hysteresis  of  the  vortex  rollup 
relative  to  the  wing  motion. 

With  bottom  views  of  the  sort  shown  In  Figure  6,  taken  simultaneously  with  corresponding  end  views, 

It  Is  possible  to  reconstruct  the  three-dimensional  character  of  the  flow  on  the  CAD  system.  Figure  7 
shows  three-dimensional  surfaces  obtained  using  the  flying  wire  technique  for  a  similar  experimental 
situation  as  that  of  Figure  7,  l.e.  a  delta  wing  pitching  about  Its  tral 1 Ing-edge.  (Note  that  the  hidden 
lines  therein  appear  as  solid  lines).  The  top  set  of  photos  shows  three  basic  views  at  relatively  high 
reduced  frequency.  In  this  set,  the  bottom  (plan)  view  shows  that  pronounced  rollup  of  the  sheet 
commences  at  about  one-third  chord  and  extends  to  about  two-thirds  chord.  This  rollup  Is  particularly 
evident  tn  the  end  view.  Upstream  of  this  region  of  rollup,  there  appear  bulges  of  opposite  curvature. 

In  the  bottom  sets  of  plots,  plan  and  end  views  are  compared  at  the  same  Instantaneous  angle  of  attack. 
Within  each  set,  from  top  to  bottom,  eases  of  static,  moderate  frequency,  and  high  frequency  are  compared. 
Efforts  are  underway  to  relate  basic  patterns  of  three-dimensional  distortion  to  the  concepts  of  hysterlsis 
and  vortex  breakdown. 

6.  FLOW  STRUCTURE  INTERPRETATION  FROM  IMAGES:  INTERACTIVE  VISUALIZATION  AND  SIMULATION 

In  any  type  of  flow  visualization  technique,  the  visualized  flow  at  any  location  represents  the 
integrated  history  of  the  marked  fluid  elements.  This  history  extends  from  the  upstream  location  where 
the  flow  Is  marked  to  the  downstream  location  where  the  flow  structure  Is  viewed.  This  concept  can  be 
used  In  conjunction  with  simple  theoretical  simulations  of  the  flow  structure  In  order  to  Identify  basic 
classes  of  vorticlty  concentration  and  circulation,  as  shown  by  Lusseyran  and  Rockwell1*.  In  essence,  they 
employ  the  velocity  eigenfunctions  from  linear  stability  theory,  In  conjunction  with  the  mean  flow,  to 
provide  simulation  of  the  vorticlty  field  of  an  unstable  wake.  Since,  In  practice,  we  can  often  estimate 
a  priori  certain  features  of  the  mean  flow,  It  should  be  possible  to  estimate  the  local  circulation  of 
the  vorticlty  field  when  flow  visualization  Is  Interpreted  In  conjunction  with  such  a  theoretical 
simulation. 

Figure  8  shows  simulation  of  the  timeline  patterns  for  a  fixed-  (upper  diagram)  and  moving-  (lever 
diagram)  timeline  marker  (e.g.,  hydrogen  bubble  wire)  for  the  case  of  a  vorticlty  field  having 
successively  Increasing  values  of  dimensionless  circulation  r*  “  I,  2,  3,  and  k.  As  Illustrated  In  the 
middle  schematic  of  Figure  8,  the  vorticlty  field,  defined  by  the  concentric  contours  Is  spatially 
periodic.  Moreover,  these  contours  of  constant  vorticlty  represent  a  neutral  disturbance,  l.e.  the 
vorticlty  Is  neither  amplifying  or  decaying.  The  amplitudes  of  the  vorticlty  contours  are  related  to  the 
overall  circulation;  however,  the  shapes  of  the  contours  are  relatively  Independent  of  the  circulation. 


Figure  7:  Various  views  of  three-dimensional  surface  constructed  from  dual  views  of  oscillating  delta 
wing  using  CAO  techniques.  Original  visualization  data  obtained  from  "flying  wire"  technique  applied  to 
wing  pitching  about  Its  tral I Ing-edge.  Top  set  of  views  corresponds  to  relatively  high  reduced  frequency 
Bottom  sets  of  views  (from  top  to  bottom)  represents  static,  moderate,  and  high  reduced  frequencies.  All 
views  In  this  figure  are  at  the  same  angle  of  attack  (Magness,  Lawson,  and  Rockwell20). 
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In  each  of  the  four  timeline  plots,  Indicated  In  the  upper  part  of  Figure  8,  the  fixed  timeline 
generator  Is  located  at  the  left  side  of  each  plot.  Comparing  these  four  timeline  patterns.  It  Is 
evident  that  the  rate  of  "roll-up"  of  the  timelines  Into  a  "vortex"  appears  to  occur  relatively  slowly 
or  rapidly  depending  upon  the  circulation  r*.  The  dots  In  these  diagrams  represent  locations  of  vortlcity 
extrema  of  the  vortlcity  field  Illustrated  In  the  middle  diagram  of  Figure  8.  The  delay  In  deformation 
of  the  timeline  pattern,  from  an  Initially  vertical  line,  Is  related  quantitatively  to  the  circulation. 

By  knowing  the  distance  from  the  timeline  generator  required  to  turn  an  Initially  vertical  timeltne  to 
an  angle  of,  say  45°,  a  direct  estimate  of  the  circulation  follows. 

Further  interpretation  of  the  flow  structure  follows  from  use  of  a  moving  timeltne  generator  as 
indicated  in  the  bottom  series  of  plots  In  Figure  8.  The  wire  velocity  Uw  In  the  downstream  section  Is  a 
given  fraction  of  the  streamwlse  velocity  U„.  In  all  diagrams,  the  vortlcity  field  Is  the  neutral  one 
given  In  the  middle  plot  of  Figure  8;  the  correspond Ing  circulation  Is  r*  ”  3.  It  Is  apparent  that  the 
apparent  rate  of  roll-up  of  the  timeline  pattern  changes  with  the  value  of  Uw,  even  though  the  vortlcity 
field  has  the  same  circulation  In  all  cases.  This  concept  can  be  used  to  our  advantage  to  Indirectly  * 
determine  the  circulation  In  conjunction  with  the  aforementioned  flow  simulation.  For  each  value  of  T  , 
there  Is  a  unique,  threshold  value  of  U.  for  which  a  timeltne  marker  moves  upstream  of  the  wire.  For  the 
parameters  here,  it  occurs  at  U*  ■  0.8  U«. 


CONCLUDING  REMARKS 

From  the  foregoing,  It  Is  evident  that  there  are  a  wide  variety  of  possibilities  for  quantitatively 
Interpreting  flow  visualization  In  conjunction  with  computer-aided  Image  processing  and  simultaneous 
acquisition  of  pressure,  lift,  etc.  Central  to  this  approach,  however,  Is  a  physically  meaningful  basis 
for  defining  the  Image.  In  fact.  If  portions  of  an  Image  or  portions  of  successive  Images  are  to  be 
correlated  with  one  another,  then  It  Is  essential  to  have  a  means  of  Identifying  these  domains.  In  a 
general  sense,  It  Is  possible  to  generate  an  array  of  streak! Ines,  pathllnes,  or  timelines  in  order  to 
accomplish  this  Identification.  However,  one  must  be  aware  of  the  limitations  Inherent  to  flow 
visualization  methods,  namely  that  the  visualized  flow  structure  may  not  necessarily  be  an  accurate 
representat Ion  of  the  underlying,  unsteady  vortlcity  field.  Insofar  as  possible.  It  Is  desirable  to 
implement  techniques  that  employ  a  time  of  flight  concept.  That  Is,  even  though  the  velocity  or  vortlcity 
may  not  be  directly  evident  from  the  Image,  the  marked  fluid  elements  In  the  Image  should  provide  a  basis 
for  determining  these  central  quantitative  parameters. 
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CONTROL  OF  LEADING-EDGE  VORTICES  ON  A  DELTA  WING1 

by 

C.  Magness,  O.  Robinson,  and  D.  Rockwell 


L  INTRODUCTION 

The  unsteady  flow  structure  of  leading-edge  vortices  on  a  delta  wing  has  been 
investigated  using  new  types  of  experimental  techniques,  in  order  to  provide  insight  into  the 
consequences  of  various  forms  of  active  control.  These  investigations  involve  global  control  of 
the  entire  wing  and  local  control  applied  at  crucial  locations  on  or  adjacent  to  the  wing. 
Transient  control  having  long  and  short  time-scales,  relative  to  the  convective  time-scale 
C/Uoo>  allows  substantial  modification  of  the  unsteady  and  time-mean  flow  structure. 

Global  control  at  long  time- scale  involves  pitching  the  wing  at  rates  an  order  of 
magnitude  lower  than  the  convective  time-scale  C/Uqq,  but  at  large  amplitudes.  The 
functional  form  of  the  pitching  maneuver  exerts  a  predominant  influence  on  the  trajectory  of 
the  feeding  sheet,  the  instantaneous  vortidty  distribution,  and  the  instantaneous  location  of 
vortex  breakdown. 

Global  control  at  short  time-scales  of  the  order  of  the  inherent  frequency  of  the  shear 
layer  separating  from  the  leading-edge  and  the  natural  frequency  of  vortex  breakdown  shows 
that  “resonant”  response  of  the  excited  shear  layer-vortex  breakdown  system  is  attainable. 
The  spectral  content  of  the  induced  disturbance  is  preserved  not  only  across  the  entire  core  of 
the  vortex,  but  also  along  the  axis  of  the  vortex  into  the  region  of  vortex  breakdown.  This 
unsteady  modification  results  in  time-mean  alteration  of  the  axial  and  swirl  velocity  fields  and 
the  location  of  vortex  breakdown. 

Localized  control  at  long  and  short  time-scales  involves  application  of  various  transient 
forms  of  suction  and  blowing  using  small  probes  upstream  and  downstream  of  the  location  of 
vortex  breakdown,  as  well  as  distributed  suction  and  blowing  along  the  leading-edge  of  the 
wing  applied  in  a  direction  tangential  to  the  feeding  sheet.  These  local  control  techniques  can 
result  in  substantial  alteration  of  the  location  of  vortex  breakdown;  in  some  cases,  it  is  possible 
to  accomplish  this  without  net  mass  addition  to  the  flow  field. 


II.  EXPERIMENTAL  TECHNIQUES 

The  unsteady  flow  structure  from  the  leading-edge  of  a  delta  wing  subjected  to  various 
forms  of  active  control  has  been  characterized  using  new  types  of  laser-diagnostic  systems  and 
image-processing  techniques.  These  methods  are  integrated  with  active  control  systems,  driven 
by  central  microcomputers.  Using  these  approaches,  it  is  possible  to  impose  active  control  of 
arbitrary  functional  form  and  examine  the  response  of  the  instantaneous  flow  structure.  The 
two-  and  three-dimensional  flow  structure  is  interpreted  with  the  aid  of  newly-released  graphics 
supercomputers . 


III.  GLOBAL  CONTROL  AT  LONG  TIME  SCALES 

The  concept  of  a  phase  shift  between  the  unsteady  motion  of  the  wing  and  the 
development  of  the  leading-edge  vortex  is  well  known.  In  qualitative  visualization  studies, 
Lambourne  et  al.  (1969),  Gad*el-Hak  and  Ho  (1985a, b,  1986),  and  Atta  and  Rockwell  (1989) 
reveal  various  features  of  the  visualized  cross-section  of  the  vortex  during  its  unsteady 
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development.  There  also  occurs  a  phase  shift  of  the  location  of  vortex  breakdown  relative  to 
the  wing  motion;  it  has  been  characterized  from  various  perspectives  by  Woffelt  (1986), 
Rockwell  et  al.  (1987),  Atta  and  Rockwell  (1987),  Reynolds  and  Abtahi  (1987),  Gilliam, 
Robinson,  Walker,  Wisser  (1987),  and  Lemay,  Batill,  and  Nelson  (1988). 

The  following  unresolved  issues  are  the  focus  of  this  investigation:  the  effect  of 
arbitrary  forms  of  pitching  maneuver  on  the  instantaneous  structure  of  the  leading-edge  vortex 
including  trajectories  of  feeding  sheets  and  distributions  of  vorticity;  the  influence  of  vortex 
breakdown  over  a  portion  of  the  cross-section  of  the  vortex;  and  the  response  of  the  axial 
location  of  vortex  breakdown  in  relation  to  all  of  these  features. 

Concerning  the  nature  of  the  instantaneous  structure  of  the  leading-edge  vortices, 
obtained  from  particle  tracking  techniques,  the  following  represent  the  major  findings: 

(i)  For  locations  upsteam  of  vortex  breakdown,  the  shape,  degree  of  concentration, 
and  the  location  of  the  maxima  of  the  instantaneous  vorticity  distribution  across 
the  vortex  core  are  quite  different  for  the  up-  and  downstrokes  of  the  continuous 
pitch-up-down  maneuvers  of  the  wing.  This  finding  emphasizes  the  importance  of 
accounting  for  the  instantaneous  cross-sectional  structure  of  the  vortex,  and  not 
simply  the  instantaneous  location  of  vortex  breakdown,  in  determining  the  overall 
loading  on  the  wing. 

(ii)  Comparison  of  the  vorticity  distribution  of  the  leading-edge  vortex  with  the 
trajectory  of  the  feeding  sheet  from  the  edge  of  the  wing  shows  the  relationship 
between  the  possible  trajectories  of  the  feeding  sheet  and  the  corresponding 
vorticity  field.  A  major  factor  is  the  occurrence  or  non-occurrence  of  vortex 
breakdown  within  the  core  of  the  vortex. 

Figure  1  shows  an  excerpt  from  the  current  investigation.  Contours  of  constant 
vorticity  were  obtained  by  direct  particle  tracking  and  image  processing  techniques.  The 
experimental  parameters  correspond  to  a  pitching  motion  of  15*  <  a  <  40*  for  a  continuous 
pitch-up-down  motion  at  a  pitching  rate  aC/2Uoo  =  0.15.  The  surface  of  the  wing  is 
indicated  by  the  bold  horizontal  line.  The  contours  of  constant  vorticity  on  the  left  side 
correspond  to  the  pitch-up  portion  of  the  maneuver,  and  those  on  the  right  side  to  the  pitch- 
down  portion.  The  differences  in  elevation,  orientation,  and  scale  of  the  vorticity  distributions 
are  evident.  They  are  dependent  upon  the  history  of  the  wing  motion  and  appear  to  be  most 
pronounced  at  the  smallest  angle  of  attack  a  =  20*. 

The  importance  of  accounting  for  vortex  breakdown  within  the  core  of  the  vortex  is 
illustrated  in  Figure  2.  Instantaneous  positions  of  the  feeding  sheets  and  contours  of  constant 
vorticity  axe  shown  for  the  same  parameters  as  in  Figure  2,  but  at  a  =  40*  for  two  different 
types  of  maneuvers.  The  shaded  black  region  represents  the  extent  of  breakdown  within  the 
vortex.  For  the  simple  pitch-down  motion,  or  =  40*  represents  the  static  condition 
immediately  preceding  the  onset  of  the  maneuver,  while  the  pitch-up-down  case  at  a  =  40* 
includes  the  integrated  history  of  the  upstroke  portion  of  the  maneuver.  It  is  evident  that  the 
positions  of  the  feeding  sheet  and  the  contours  of  vorticity  are  substantially  different  for  these 
two  cases. 

The  structure  of  the  leading-edge  vortex  at  a  given  cross-section  must,  of  course,  be 
considered  in  conjunction  with  the  axial  movements  of  the  location  of  vortex  breakdown. 
Magness,  Robinson,  and  Rockwell  (1989)  preliminarily  addressed  the  effect  of  the  type  of 
pitching  maneuver  on  the  general  response  of  the  vortex  breakdown  location  as  a  function  of 
angle  of  attack.  Recent  studies  have  focussed  on  the  vortex  response  to  different  classes  of 
maneuver,  and  the  detailed  structure  of  the  leading-edge  vortices.  Regarding  the  response  of 
the  location  of  the  vortex  breakdown,  the  major  findings  are: 
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(i)  Continuous  pitch  up-down  motions  of  the  wing  can  preclude  occurrence  of  vortex 
breakdown  on  the  upstream  portion  of  the  wing,  relative  to  that  occurring  for 
simple  pitch-up  and  pitch-down  motion  where  relaxation  processes  having  long 
time  scales  are  allowed  to  occur. 

(ii)  For  the  continuous  pitch-up-down  maneuver  of  the  wing,  the  consequence  of  not 
allowing  the  vortex  breakdown  to  relax  to  its  equilibrium  state  is  to  produce 
upstream  movements  of  the  vortex  breakdown  location  towards  the  apex  for 
initial  decreases  in  angle  of  attack  or. 

(iii)  Combinations  of  simple  ramp-type  motions  to  form  a  hybrid  pitching  motion 
produce  overshoots  of  the  static  characteristic  of  vortex  breakdown  location  versus 
angle  of  attack,  beyond  that  attainable  with  any  of  the  simple  ramp  motions 
alone. 

Figure  3  shows  plots  characterizing  the  first  two  of  these  three  principal  findings  at  two 
extreme  values  of  reduced  frequency.  This  sort  of  characterization  of  the  breakdown  location 
serves  as  a  basis  for  detailed  investigations  of  the  flow  structure  of  the  leading-edge  vortices. 


IV.  GLOBAL  CONTROL  AT  SHORT  TIME  SCALES 

Perturbation  of  a  delta  wing  in  the  pitching  mode  at  sufficiently  high  frequency  and 
very  low  amplitude  allows  control  of  the  detailed  flow  structure  of  the  leading-edge  vortex.  In 
essence,  the  vortex  development  and  breakdown  on  a  delta  wing  involves  two  classes  of 
characteristic  frequencies:  the  inherent  instability  frequency  of  the  shear  layer  from  the 
leading-edge;  and  the  frequency  at  which  vortex  breakdown  occurs.  The  major  issues  here  are: 
the  structure  of  the  perturbed  feeding  sheet;  the  nature  of  the  perturbed  onset  of  vortex 
breakdown;  and  the  corresponding  alteration  of  the  time-mean  vortex  flow. 

Simple  considerations  of  hydrodynamic  instability  show  that  the  processes  of 
disturbance  amplification  in  the  shear  layer  and  in  the  vortex  core  during  the  breakdown 
process  are  receptive  to  a  wide  range  of  excitation  frequencies.  As  a  consequence,  it  is  possible 
to  attain  “resonant”  excitation,  leading  to  large  alteration  of  the  separating  shear  layer  from 
the  edge  of  the  wing  and  the  breakdown  of  the  vortex  core.  The  preliminary  phase  of  this 
investigation  was  reported  by  Rockwell  et  al.  (1987).  This  work  is  described  in  its  completed 
form  by  Kuo,  Magness,  and  Rockwell  (1989). 

The  principal  findings  of  this  investigation  are,  in  short: 

(i)  Small  amplitude  perturbations  of  the  leading-edge  lead  to  substantial  alteration  of 
the  structure  of  the  shear  layer  separating  from  it  without  occurrence  of  the 
classical  mechanism  of  small-scale  vortex  coalescence. 

(ii)  The  spectral  content  of  the  disturbance  induced  in  the  shear  layer  separating  from 
the  leading-edge  is  preserved  not  only  across  the  core  of  the  vortex,  but  also  along 
the  streamwise  extent  of  the  core  into  the  region  of  vortex  breakdown. 

(iii)  Substantial  alteration  of  the  time-mean  characteristics  of  the  leading-edge  vortex 
include  changes  in  the  axial  and  swirl  velocity  fields  and  modification  of  the 
location  of  vortex  breakdown. 

Selected  excerpts  describing  certain  of  the  foregoing  phenomena  are  given  in  Figures  4 
through  6.  Figure  4  shows  the  visualization  obtained  by  locating  a  vertical  hydrogen  bubble 
wire  along  the  leading-edge  of  the  wing.  The  laser  sheet  that  illuminated  the  marker  bubbles 
was  translated  to  the  downstream  locations  x/C  indicated  in  the  photos.  Excitation  frequency 
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fe  is  normalized  with  respect  to  the  inherent  instability  frequency  f-  of  the  shear  layer 
separating  from  the  leading-edge.  Large-scale  vortical  structures  are  induced  over  the  cross- 
section  of  the  vortex  in  the  presence  of  excitation  at  the  inherent  instability  frequency  of  the 
feeding  sheet.  No  small-scale  vortex  coalescence  occurs. 

Figure  5  shows  spectra  of  the  stream  wise  component  u,  i.e.  SQ,  taken  at  various 
locations  upstream  and  downstream  of  the  onset  of  vortex  breakdown.  The  edge  excitation 
frequency  is  normalized  by  the  inherent  vortex  breakdown  frequency  f^.  For  excitation  at 
the  first  harmonic  of  the  vortex  breakdown  frequency,  i.e.  at  ^/f,  =  2,  the  spectral  content 
shows  predominance  of  the  excitation  frequency  and  its  associated  higher  harmonics  in  regions 
before  and  after  occurrence  of  vortex  breakdown.  In  this  case,  the  higher  harmonic  content 
persists  well  downstream  of  the  onset  of  breakdown.  For  excitation  at  f^/f,  =  1»  there  also 
occur  a  large  number  of  higher  harmonics  due  to  the  strong  nonlinearity  of  the  shear  layer 
response.  This  spectral  content  is  maintained  over  the  entire  cross-section  of  the  vortex  core 
prior  to  the  occurence  of  breakdown,  emphasizing  the  nonoccurence  of  vortex-vortex 
interactions  (i.e.  coalescence)  in  the  shear  layer  as  it  is  wrapped  inwards  toward  the  center  of 
the  core.  Downstream  of  vortex  breakdown,  the  predominant  excitation  peak  at  fe/f^  —  1 
persists,  but  the  coherent  higher  harmonic  components  are  attenuated. 

Figure  6  shows  contours  of  constant  mean  axial  velocity  0  and  constant  fluctuating 
velocity  u  over  the  entire  cross-section  of  the  leading-edge  vortex  at  values  of  excitation 
frequency  f^,  relative  to  the  inherent  vortex  breakdown  frequency  f^,  i.e.  f^/f^  =  1  (left 
column)  and  2  (right  column).  The  effect  of  the  matched  excitation  at  f^/f^  =  1  is  to  induce 
large  amplitude  fluctuations  in  the  separating  shear  layer  surrounding  the  core  of  the  vortex, 
located  at  the  peak  of  the  contours  of  constant  Q  and  designated  by  the  symbol  -{-.  At 
f«/f^  =  2,  the  location  of  the  core  of  the  vortex  moves  downward  towards  the  surface  of  the 
wing  and  outward  towards  the  leading-edge.  The  maximum  amplitude  of  the  fluctuation  u  is 
coincident  with  the  location  of  the  core  of  the  vortex.  This  coincidence  of  the  maxima  of  H 
and  u  corresponds  to  the  early  onset  of  vortex  breakdown  at  the  higher  excitation  frequency 

Vfb  =  2. 

V.  LOCAL  CONTROL  AT  MODERATE  AND  LONG  TIME  SCALES 

Local  control  involves  localized  injection  or  suction  of  the  flow  at  defined  locations  in 
the  flow  field  and/or  the  surface  of  the  wing.  In  a  practical  sense,  this  can  be  achieved  by  use 
of  small  probes,  whose  tips  are  located  at  crucial  locations  in  the  vortex  core,  or  slits  along  the 
leading-edge  of  the  wing.  In  essence,  these  techniques  simulate  localized  point  sources/sinks  or 
distributions  of  them.  The  major  issues  here  are:  determination  of  the  most  sensitive  location 
of  the  applied  control;  and  optimization  of  the  functional  form  of  the  unsteady  control  in  the 
form  of  blowing/suction. 

For  the  case  of  localized  blowing  along  the  leading-edge  of  the  wing,  Wisser,  Iwanski, 
Nielson,  and  Ng  (1988)  most  recently  have  revealed  an  increase  in  length  of  the  vortex  core 
prior  to  breakdown  and  an  increase  in  lift  acting  on  the  wing.  Not  until  this  past  year  has  the 
case  of  localized  suction  been  explored;  such  simulations  of  a  localized  sink  are  described  by 
Parm  cter  and  Rockwell  (1989).  Location  of  a  probe  in  the  region  downstream  of  vortex 
breakdown  allows  efficient  restabilization  of  the  vortex  core.  Among  the  principal  findings  are: 

(i)  Locations  of  the  simulated  point  sink  downstream  of  the  occurrence  of  vortex 
breakdown  produces  stabilization  of  the  core;  such  stabilization  is  attainable  at 
relatively  low  values  of  dimensionless  suction  coefficient  C^.  The  transient 
response  time  of  the  stabilization  process  due  to  an  imposed  transient  (unsteady 
sink  flow)  scales  as  the  magnitude  of  the  imposed  transient  suction. 

(ii)  Hysteresis  effects  occur  due  to  relaxation  of  the  vortex  breakdown  (on  a 
stationary  wing)  after  abrupt  onset  or  cessation  of  suction.  These  hysteresis 
effects  simulate  those  on  a  pitching  delta  wing. 
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Localized  control  involving  simulations  of  distributed  sources/sinks  in  the  form  of  a 
blowing/suction  slit  along  the  leading-edge  have  received  little  attention  except  for  the  steady 
blowing  experiments  of  Wood  and  Roberta  (1987),  Wood,  Roberts  and  Lee  (1987)  and 
Roberts,  Hesselink,  Kroo,  and  Woods  (1987),  and  the  (high  frequency)  sinusoidal 
perturbations  employed  in  the  investigation  of  Gad-el-Hak  and  Blackwelder  (1987).  The 
consequences  of  this  type  of  control  on  the  structure  of  the  large-scale  vortex  have  remained 
unexplored.  Moreover,  the  possible  modification  of  the  nature  and  location  of  onset  of  vortex 
breakdown  has  not  been  pursued.  Important  considerations  in  our  recent  investigations 
include  not  only  the  case  of  steady  blowing,  but  also  the  corresponding  case  of  steady  suction 
and,  most  significantly,  the  case  of  cyclic  blowing  and  suction.  The  major  findings  of  this 
investigation  are: 

(i)  Both  steady  suction  and  steady  blowing  are  effective  at  low  values  of  C^,  i.e.  both 
result  in  lengthening  of  the  vortex  core  prior  to  the  onset  of  breakdown. 

(ii)  The  most  effective  and  robust  control  involves  cyclic  suction  and  blowing  at  an 
appropriate  frequency.  This  approach  involves  no  net  mass  addition  to  or  from 
the  flow. 

The  use  of  cyclic  blowing  and  suction  applied  tangentially  in  the  form  of  a  jet  Vj(t)  at 
the  rounded  leading-edge  is  represented  in  Figure  7;  it  is  compared  with  the  case  of  no 
blowing/suction,  i.e.  V-(t)  =  0.  (These  data  were  acquired  by  Professor  W.  Gu,  a  member  of 
our  research  group.)  Comparison  of  these  velocity  fields  of  Figure  7  suggests  that  application 
of  the  control  results  in  restabilization  of  the  vortex  from  a  stalled  condition  to  a  well-defined, 
large-scale  vortical  structure  and  downward  deflection  of  the  separation  streamline  from  its 
approximately  horizontal  position.  These  trends  are  associated  with  downstream  movement  of 
the  location  of  vortex  breakdov  n. 
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Figure  1:  Instantaneous  contours  of  constant  vorticity  at  midchord  for  continuous  pitch-up* 
down  maneuver  of  delta  wing.  Sweep  angle  =  75*;  pitch  rate  6rC/2U  =  0.15;  pitching  axis  at 
midchord. 
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Figure  2:  Instantaneous  contours  of  constant  vorticity  and  positions  of  feeding  sheet  at 
a  =  40*  for  continuous  pitch-up-down  and  pitch-down  maneuvers.  Sweep  angle  =  75*;  pitch 
rate  aC/2U  =  0.5;  pitching  axis  at  midchord. 
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Figure  3:  Instantaneous  location  of  vortex  breakdown  as  function  of  angle  of  attack  for  three 
basic  types  of  delta  wing  maneuvers  and  two  extreme  values  of  dimensionless  pitching  rate. 
Sweep  angle  of  delta  wing  ss  75*. 


Figure  4:  Visualization  of  flow  structure  of  separating  shear  layer  at  several  cross-sections 
along  leading-edge  of  delta  wing.  Ratio  of  excitation  frequency  fe  to  inherent  instability 
frequency  fj  of  separating  shear  layer  is:  0  (left  column);  0.5  (middle  column);  and  1.0  (right 
column).  Angle  of  attack  o  s  5  +  Oq  sin  2rfet;  5  s  20*,  Oq  =s  1*. 
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Figure  5:  Overview  of  evolution  of  spectral  content  of  vortex  core  upstream  and  downst-»am 
of  vortex  breakdown.  Ratio'  of  excitation  frequency  f^  to  vortex  breakdown  frequency  has 
values  indicated,  fj/f^  =  2. 
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FisTi-e  6:  Distribution  of  mean  U  component  of  axial  velocity  over  cross-section  of  vortex 
upstream  of  occurrence  of  vortex  breakdown.  Data  acquired  at  reference  station  x/C  =  0.3S 
upstream  of  vortex  breakdown.  Ratio  of  excitation  frequency  f*  to  frequency  f.  of  inherent 
vortexbreakdown  is  =  1  (left  column)  and  f^/fb  =  2  (right  column).  Ratio  of  inherent 
instability  frequency  f.  of  separating  shear  layer  to  freouenev  f,  i.»  f.  /f.  —  o 
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Figure  7:  Comparison  of  velocity  fields  over  cross-section  of  half  delta  wing  at  o  =  45*  for 
cases  with  and  without  sinusoidal  suction/blowing  applied  tangentially  at  rounded  leading- 
edge.  Laser  sheet  defining  cross-section  of  visualized  vortex  located  at  approximately  one-third 


